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Introduction 


This  report  summarizes  the  works  done  for  the  second  year  of  this  3  year  grant,  between  May  1,  2004  and 
April  30,  2005. 

1.  Scope  of  work 

The  scope  of  the  work  includes  software  and  instrument  development  outline  below. 

Task  1.  To  develop  instruments  to  measure  the  optical  properties  (1-6  months) 

a.  Construct  a  motorized  detector  probe  (1-3  months). 

b.  Verify  the  feasibility  software  of  the  detector  probe  in  phantom  (4-6  months). 

c.  Develop  a  16-channel  in-vivo  light  dosimetry  system  capable  of  interfacing  with  the 
motorized  detector  probe  (1-6  months). 

Task  2.  To  develop  software  for  light  dose  calculation  and  optimization,  and  accumulate 
clinical  data  (6-18  months) 

a.  Develop  a  kernel  based  dose  calculation  algorithm  and  compare  with  optical  phantom 
measurements  (6-12  months). 

b.  Verify  the  accuracy  of  the  calculation  of  light  dosimetry  in  patients  (6  -  24  months). 

c.  Develop  the  software  to  optimize  the  weights  of  PDT  light  sources  (1-18  months). 

d.  Improve  in  situ  fluorescence  measurements  of  photosensitizer  (1-18  months). 

Task  3.  To  develop  instruments  and  software  for  integrated  real-time  in-vivo  light  monitoring 

and  optimized  light  delivery  (18-36  months) 

a.  Develop  an  8-way  beamsplitter  system  for  light  delivery  (18-24  months). 

b.  Interface  the  in-vivo  light  dosimetry  system  with  the  calculation  engine  (18-24 
months). 

c.  Interface  the  fluorescence  and  absorption  measurement  to  determine  drug 
concentration  and  tissue  oxygenation  with  the  PDT  dose  calculation  engine  (18-24 
months). 

d.  Refine  the  integrated  system  for  real-time  light  monitoring  and  optimization  in 
patients  (24  -  36  months). 

For  the  scientific  research,  the  scientific  objectives  are  outlined  below  as  specific  aims: 

1:  Characterize  the  optical  properties  (scattering  and  absorption  coefficients)  of  the  prostate  gland  in  vivo : 
2:  Add  the  ability  to  calculate  light  fluence  rate  within  the  prostate  volume  during  PDT  to  the  current  light 
dosimetry  system  that  measures  light  fluence  at  selected  points: 

3:  Measure  photosensitizer  (MLu)  drug  concentration  and  tissue  oxygenation  measured  with  optical 
fluorescence  and  absorption  spectra  and  incorporate  these  values  in  a  determination  of  PDT  dose: 

4:  Develop  optimization  methods  for  placement  of  the  light  source  to  achieve  uniform  PDT  dose  in  the 
prostate  gland: 

2.  Summary  of  accomplishment 

For  the  instrument  and  software  development  works,  we  have  finished  work  on  Task  2  this  year. 
Combined  work  done  in  the  last  two  years,  we  have  completed  Task  1  and  Task  2.  We  are  making  good 
progress  for  Task  3  a  -  c. 

For  the  scientific  research,  we  have  completed  specific  aims  1  -3. 
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Development  of  instrument  to  measure  the  distribution  of  optical  properties  and  light  fluence  rate 


Following  the  successful  completion  of  the  characterization  of  optical  properties  in  12  patients  last  year, 
we  have  published  the  results  (for  14  patients)  in  a  peer-reviewed  journal  (Ref.  1).  We  have  improved  the 
motorized  detector  probe  to  characterize  the  optical  properties  at  different  depths  in  a  catheter  (see 
Section  A  below).  This  new  design  along  with  improved  software  (Section  B)  allows  us  to  measure  the 
distribution  of  optical  properties  in  prostate  (Ref.  2).  We  have  now  characterized  the  accuracy  of 
determining  the  optical  properties  using  the  instrument  (Ref.  3). 

A.  Motorized  detector  probe 

The  motorized  detector  probe  system  has  been  expanded  from  the  single-channel,  single-motor  device 
reported  last  year  to  a  two-motor  system  capable  of  controlling  up  to  5  detectors  simultaneously.  This 
device  allows  one  light  source  and  up  to  four  detectors  to  be  used  to  make  optical  property  measurements 
in  four  regions  of  a  single  prostate  quadrant  at  the  same  time.  The  device  is  shown  in  figure  1.  An 
additional  enhancement  of  the  device  is  the  addition  of  plastic  extensions  onto  the  catheters  used  for 
detection.  These  extensions  are  held  by  a  clamp  mounted  on  the  front  of  the  motor  assembly,  and  allow 
the  optical  fibers  to  be  translated  along  curved  paths  without  binding  or  interference,  greatly  improving 
the  accuracy  of  positioning. 


Figure  1:  Photograph  of  the  multichannel  step  motor  assembly  used  to  move  multiple  detectors  and  a  single 
source  simultaneously. 


B.  Software 

The  introduction  of  the  two-motor  positioning  system  has  been  accompanied  by  an  enhancement  in  our 
dosimetry  software,  allowing  us  to  use  the  motorized  positioner  to  characterize  the  fluence  rate 
distribution  during  treatment.  Previously,  the  maximum  value  in  each  quadrant  (to  which  the  prescription 
is  specified)  was  found  by  manual  positioning  of  the  fiber.  In  the  new  system,  the  detection  fiber  in  the 
quadrant  being  treated  is  controlled  by  the  motorized  positioner.  The  detector  is  scanned  along  its 
catheter  under  computer  control  without  interruption  of  dosimetry  measurements.  The  resulting  fluence 
rate  profile  is  stored  for  future  analysis,  and  the  motor  automatically  moves  the  detector  to  the  position  of 
maximum  fluence  rate.  This  not  only  greatly  improves  the  accuracy  and  speed  with  which  the  point  of 
maximum  fluence  rate  can  be  found,  but  also  gives  a  quantitative  measure  of  the  fluence  rate  profile 
within  the  quadrant.  Because  this  measurement  requires  only  one  of  the  two  motors,  the  other  can  be 
used  to  make  profile  measurements  in  adjacent  quadrants.  This  allows  characterization  of  the  diffusion  of 
light  from  one  quadrant  into  another  and  also  ensures  that  the  dosimetric  measurements  will  accurately 
reflect  the  total  dose  given  to  each  quadrant,  including  contributions  from  adjacent  quadrants. 
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C.  Instrument  accuracy 

The  instrument  accuracy  has  been  characterized  in  optical  phantom.  We  find  that  |_ia  and  ps’  can  be 
determined  by  this  method  with  a  standard  (maximum)  deviation  of  8%  (15%)  and  18%  (32%)  for  pa  and 
ps’,  respectively  (see  Abstract  and  Fig  9  of  Appendix  3,  Ref.3). 

Development  of  light  dose  calculation  and  optimization  and  accumulate  clinical  data 

We  have  developed  a  kernel  model,  which  allows  predict  light  fluence  rates  in  human  prostates  during 
PDT.  We  have  examined  the  model  by  comparing  the  calculations  with  measurements  in  14  patients. 

A.  Kernel  model 


A  kernel  model  for  light  fluence  rate  calculation  is  developed  for  each  linear  light  source  used  in  a 
treatment,  which  is  based  on  the  solution  of  the  diffusion  equation  for  a  point  source. 


i  =  z 
1  =  1 
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4wan 
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where  s  is  the  source  strength  of  each  segment  along  the  linear  light  source,  Ax,  is  the  division  along  the 


direction  parallel  to  the  orientation  of  the  linear  source,  /ieff  is  the  effective  attenuation  coefficient,  f.ia  is 

the  absorption  coefficient,  and  r,  is  the  distance  between  the  source  segment  and  the  detector.  The 
summation  in  Eq.  (1)  is  over  all  the  segments  composing  the  linear  source.  The  fluence  rate  at  a  point  in  a 
prostate  is  the  summation  of  the  fluence  rates  generated  by  each  linear  source.  In  a  calculation,  ultrasound 
images  of  a  prostate  are  imported  and  the  calculated  light  fluence  rate  distribution  is  superimposed  on  the 
contours  of  the  prostate. 


B.  Application  and  examination  of  the  model 


B.l.  Application  in  treatment  planning 


NM,  #16,  z  =  0.5  cm,  100  J/cm' 


350  r 


RLQ,  LLQ,  Peak  =  327.8  mW/crr/ 

f  \ 

f  \ 

Urethra,  Peak  =>256.6  mW/cm2 


2 

Z  (cm) 

(a)  (b) 

Figure  2:  Calculated  light  fluence  rate  distribution  superposed  on  an  ultrasound  image  of  prostate  cross- 
section  at  depth  of  0.5  cm,  which  was  used  in  PDT  planning.  The  outer  white  line  indicates  the  predicted 
100%  isodose  line  and  the  inner  line  is  the  contour  of  the  prostate.  The  “o”  symbols  represent  the  source 
positions  and  the  predicted  light  fluence  rates  at  “x”  symbol  positions  are  displayed  next  to  the  symbols,  (b) 
Light  fluence  rates  along  the  z  direction  at  the  detector  positions  and  in  the  urethra. 

With  the  kernel  model,  we  have  been  able  to  calculate  light  fluence  rates  in  a  three-dimensional  prostate 
volume.  And  we  have  applied  the  calculation  in  current  PDT  treatment  planning. 
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Figure  2(a)  shows  predicted  light  fluence  rate  distribution  (the  outer  thin  line)  superposed  on  an 
ultrasound  image  of  a  prostate  cross-section,  which  is  a  calculation  result  used  in  a  treatment  plan.  The 
outer  white  line  shows  the  isodose  line  of  100  J/cm2  and  the  inner  line  shows  the  contour  of  the  prostate. 
The  “o”  symbols  represent  the  source  positions  and  the  “x”  symbols  represent  the  detector  positions,  next 
to  which  the  corresponding  light  fluence  rates  are  displayed.  Figure  2(b)  shows  light  fluence  rates  in  the  z 
direction  at  the  detector  positions  and  in  the  urethral  as  well.  We  actually  can  calculate  light  fluence  rates 
anywhere  in  the  prostate.  The  treatment  time  is  determined  using  the  predicted  light  fluence  rates.  The 
calculation  is  fast,  taking  only  a  few  seconds.  We  are  planning  to  apply  the  calculation  in  real-time 
dosimetry. 


Patient  number  Patient  number 


Figure  3:  Comparison  between  light  fluence  rates  calculated  using  measured  optical  properties  and 
estimated  optical  properties  with  actual  source  strengths. 

B.2.  Examination  of  the  model:  comparison  of  calculation  and  measurement 

To  examine  the  model  and  check  the  accuracy  of  the  calculation,  we  have  made  comparisons  between  the 
calculations  and  the  measurements  in  14  patients.  Calculations  were  done  under  different  conditions,  e.g., 
different  source  strengths  and  different  optical  properties.  Overall,  the  calculations  were  categorized  into 
two  kinds:  (1)  with  assumed  parameters,  and  (2)  with  actual  parameters  used  in  each  treatment.  The 
former  was  like  what  we  have  been  doing  in  treatment  planning.  In  that  case,  we  assumed  that  the  optical 
properties  were  homogeneous  in  a  prostate  and  we  used  the  mean  optical  properties  for  the  whole  prostate 
in  the  calculation,  which  was  the  mean  of  all  the  optical  properties  measured  in  the  patients  treated 
previously,  i.e.,  pa  =  0.3  cm"1  and  pf=  14  cm'1.  Also,  we  assumed  the  source  strength  to  be  150 

mW/cm,  which  was  the  source  strength  used  in  most  of  the  treatments  except  in  the  case  of  bleeding.  In 
the  latter  case,  we  used  actual  source  strengths  and  incorporated  heterogeneity  of  optical  properties  of  a 
prostate,  i.e.,  we  used  optical  properties  measured  in  different  regions  in  a  prostate.  We  had  measured 
optical  properties  in  most  of  the  patients’  prostates. 

During  each  treatment,  we  measured  light  fluence  rates  in  four  regions  in  a  prostate,  which  we  call 
quadrants  (RUQ,  LUQ,  RLQ,  and  LLQ).  We  compared  our  calculation  results  with  these  measured  data 
in  each  quadrant.  A  percentage  error  is  defined  as  the  percentage  difference  between  the  calculated  and 
the  measured  light  fluence  rates,  i.e.,  (Calculation-Measurement)/Measurement.  Flere  we  present  some  of 
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the  results.  Figure  3  shows  comparison  between  light  fluence  rates  calculated  using  measured  optical 
properties  and  the  mean  optical  property,  with  actual  source  strengths.  The  maximum  error  is  over  600% 
(LUQ,  patient  3).  This  large  error  may  be  due  to  the  inaccuracy  of  the  light  fluence  measurement  rather 
than  the  calculation,  because  the  measured  light  fluence  rate  in  this  case  was  unexpectedly  low.  It  is 
possible  that  this  anomalously  low  reading  was  caused  by  bleeding  around  the  catheter.  Since  the 
measurements  of  optical  properties  and  light  fluence  rates  in  patients  11-17  were  well  controlled,  the 
measurement  data  in  these  cases  are  more  reliable. 

The  histograms  of  the  comparison  results  for  patients  1 1-17  are  shown  in  Fig.  4,  which  are  the  statistics  of 
all  the  comparisons,  i.e.,  the  comparisons  of  the  calculations  in  all  the  quadrants  for  these  patients.  In  the 
figure,  we  used  the  absolute  values  of  the  percentage  errors.  Fig.4a  shows  the  results  of  the  calculations 
using  actual  source  strengths  and  measured  optical  properties.  The  mean  of  the  absolute  percentage  errors 
is  59.5%.  Fig.  3b  shows  the  results  of  the  calculations  using  actual  source  strengths  and  the  mean  optical 
property,  where  the  mean  of  the  absolute  percentage  errors  is  70.0%,  which  is  larger  than  that  in  Fig.  4a. 
This  indicates  that  accurate  optical  properties  are  needed  to  improve  the  accuracy  of  the  calculation.  Fig. 
4c  shows  the  results  of  the  calculations  using  a  source  strength  of  150  mW/cm  and  the  mean  optical 
properties,  analogous  to  the  calculations  used  in  treatment  planning.  The  mean  of  the  errors  is  88.7%,  not 
significantly  worse  than  the  59.5%  error  found  when  using  actual  source  strengths  and  measured  optical 
properties. 


Statistics  of  the  comparison  results,  patients  11-17  Statistics  of  the  comparison  results,  patients  11-17 


Absolute  percentage  error  (%)  Absolute  percentage  error  (%) 


Statistics  of  the  comparison  results,  patients  11-17 


Absolute  percentage  error  (%  j 


(a)  (b)  (c) 

Figure  4:  Histograms  of  the  comparison  results  in  patients  11-17  for  calculations  using:  (a)  actual  source 
strengths  and  measured  optical  properties;  (b)  actual  source  strengths  and  the  mean  optical  properties  (pa  = 
0.3  cm'1  and  ps’=  14  cm'1);  and  (c)  assumed  source  strength  of  150  mW/cm  and  the  mean  optical 
properties. 


C.  Conclusions 


We  have  developed  and  applied  a  kernel  model  in  human  prostate  PDT  treatment.  With  the  model,  we  are 
able  to  calculate  light  fluence  rates  inside  a  human  prostate  volume.  We  are  applying  the  calculation  in 
treatment  planning  and  real-time  treatment  guiding.  We  have  compared  the  calculations  with  the 
measurements  in  14  patients.  The  comparisons  with  those  measurements  which  had  more  reliable  data 
show  that,  the  standard  uncertainty  is  59.5%  when  actual  source  strengths  and  measured  optical  properties 
were  used  in  the  calculations,  and  the  standard  uncertainty  is  88.7%  when  the  assumed  source  strength  of 
150  mW/cm  and  the  mean  optical  properties  were  used  in  the  calculations. 

Characterization  of  absorption  spectra  in  vivo 

The  absorption  spectroscopy  methods  described  previously  (Ref.  4)  have  been  expanded  by  the  use  of  the 
multichannel  motorized  positioning  system  described  above.  As  before,  the  diffuse  fluence  rate  spectrum 
is  measured  at  a  number  of  distances  from  a  white  light  source  of  known  intensity.  The  positions  of  the 
source  and  detector  are  controlled  independently  by  the  multichannel  motorized  positioner.  The  source 
position  is  fixed,  and  measurements  are  made  at  different  detector  positions.  The  source  is  then  moved  by 
a  predetermined  distance,  and  the  process  is  repeated.  By  repeating  this  process  multiple  times,  we  can 
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map  the  distribution  of  absorption  and  scattering  spectra  within  the  prostate.  Each  absorption  spectrum 
can  then  be  analyzed  using  a  singular  value  decomposition  (SVD)  algorithm  to  determine  the 
concentrations  of  oxyhemoglobin  (Hb02),  deoxyhemoglobin  (Hb),  and  MLu,  as  shown  in  Fig.  5a  (Ref.  4). 


Figure  5:  (a)  Components  of  absorption  spectra  acquired  prior  to  PDT  in  vivo.  The  measured 
data  are  labeled  as  symbols  ('o’)  and  the  components  are  labeled  as  lines.  The  components 
include:  MLu,  water  (H20),  deoxyhemoglobin  (Hb),  and  hemoglobin  (Hb02).  The  ‘Fourier’ 
component  is  a  Fourier  series  designed  to  account  for  unknown  absorbers  or  fluorophores.  (b) 

The  corresponding  in  vivo  distribution  of  (a)  St02,  blood  volume,  and  MLu  concentration 
determined  using  the  absorption  spectra. 

Combining  these  measurements  allows  us  to  construct  a  map  of  the  distribution  of  the  sensitizer,  the  total 
hemoglobin  concentration,  and  the  hemoglobin  oxygen  saturation  (St02)  (Ref  2).  One  such  distribution, 
obtained  in  vivo  from  a  human  prostate,  is  shown  in  figure  4(b).  This  measurement  demonstrates  that  the 
heterogeneity  in  optical  properties  observed  in  our  single-wavelength  measurements  is  related  to  the 
heterogeneity  in  blood  and  sensitizer  distribution. 

Fluorescence  spectroscopy  in  vivo 

* 

The  local  concentration  of  sensitizer  can  also  be  measured  via  fluorescence  spectroscopy  (Ref.  5).  We 
have  developed  a  measurement  technique  to  characterize  the  fluorescence  of  tissue  using  a  side-firing 
fiber  and  the  multichannel  motorized  positioner  described  above.  A  schematic  of  this  device  is  shown  in 
Fig.  6.  The  fluorescence  system  consists  of  an  optical  fiber  terminated  in  a  beveled  end  that  emits  and 
collects  light  at  a  right  angle  to  its  axis.  This  sidefiring  fiber  is  connected  to  a  dichroic  beamsplitter  that 
delivers  465-nm  excitation  light  provided  by  a  light-emitting  diode  (LED)  to  the  fiber.  The  fluorescence 
excited  by  this  light  is  then  collected  by  the  fiber  and,  by  virtue  of  its  longer  wavelength,  passes  through 
the  dichroic  beamsplitter  to  the  spectrograph  and  CCD. 

The  spectra  captured  by  the  CCD  are  analyzed  using  the  same  SVD  algorithm  used  to  analyze  the 
absorption  spectra  described  above.  In  this  case,  the  basis  spectra  used  are  the  inherent  background 
fluorescence  of  the  catheter  and  optical  fiber  and  the  fluorescence  emission  spectrum  of  MLu.  The  MLu 
fluorescence  is  corrected  for  variations  in  lamp  intensity  by  dividing  by  the  inherent  background 
fluorescence  and  then  scaled  by  an  empirically  determined  value.  By  moving  the  sidefiring  fiber  along 
its  catheter,  we  can  sample  the  distribution  of  MLu  concentration  within  the  prostate.  The  distributions  in 
all  four  quadrants  of  a  typical  prostate  are  shown  in  figure  7.  Profiles  acquired  before  and  after  irradiation 
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are  shown.  In  general,  the  shape  of  the  distribution  dose  not  appears  to  change  with  irradiation,  indicating 
that  there  is  little  redistribution  of  MLu  during  treatment. 


465  nm  LED 


Figure  6:  Fluorescence  spectroscopy  setup.  The  computer  that  acquires  and  stores  fluorescence  spectra  also 
controls  the  position  of  the  detection  fiber  via  a  step-motor  positioner  (not  shown).  The  dichroic  allows  a 
single  fiber  to  deliver  excitation  light  and  collect  emitted  fluorescence. 

LUQ  RUQ 


Figure  7:  MLu  fluorescence  profiles  acquired  in  the  four  quadrants  of  a  single  prostate  before  (dashed 
lines)  and  after  (solid  lines)  PDT  treatment.  Frames  a  through  d  depict  the  results  from  the  right  upper, 
left  upper,  left  lower,  and  right  lower  quadrants,  respectively,  as  indicated  in  the  panel  titles. 


We  have  described  above  three  independent  methods  of  obtaining  MLu  profiles  in  vivo ,  namely 
fluorescence  spectroscopy,  absorption  spectroscopy,  and  single-wavelength  absorption  measurements.  In 
one  patient,  we  have  obtained  all  three  measurements  in  a  single  quadrant  before  and  after  PDT.  The 
distributions  obtained  from  this  patient  are  shown  in  figure  8.  The  agreement  in  shape  and  amplitude 
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among  the  three  methods  confirms  that  all  three  are  sensitive  to  the  same  distribution.  The  single¬ 
wavelength  and  spectroscopic  absorption  measurements  therefore  provide  verification  of  the  accuracy  of 
the  MLu  concentration  measured  by  fluorescence  spectroscopy. 


Figure  8:  MLu  fluorescence  (filled  circles),  absorption  spectroscopy  (triangles),  and  732-nm  absorption 
(squares)  profiles  as  functions  of  position  within  the  right  upper  quadrant  of  patient  #13.  The  profiles 
measured  before  (panel  a)  and  after  (panel  b)  PDT  treatment  are  similar  in  shape,  but  indicate  some 
photobleaching,  especially  in  the  center  of  the  prostate. 


Finally,  we  have,  in  addition  to  the  three  independent  optical  measurements  described  above,  an 
independent  measurement  of  tissue  MLu  concentration  based  on  an  ex  vivo  fluorescence  analysis  of  tissue 
biopsies  taken  before  and  after  PDT.  As  these  biopsies  are  taken  in  only  one  position  in  the  prostate,  we 
cannot  obtain  a  distribution  of  MLu  from  them,  however  we  can  compare  the  ex  vivo  measurement  with 
the  means  of  the  values  obtained  from  the  three  other  methods.  The  results  of  this  comparison  for  five 
patients  are  shown  in  figure  9.  The  four  measurements  agree  remarkably  well  considering  the 
heterogeneity  obvious  in  the  distributions  shown  in  figure  8  and  the  variations  among  prostates. 
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Figure  9:  MLu  concentration  in  5  human  prostates  measured  by  fluorescence  spectroscopy  (circles),  ex  vivo 
fluoroscopy  (squares),  absorption  spectroscopy  (diamonds),  and  732  nm  absorption  (triangles).  Both  pre- 
PDT  (open  symbols)  and  post-PDT  data  (filled  symbols)  are  plotted.  The  dotted  lines  indicate  the  mean  of 
all  measurements  in  each  patient,  and  the  shaded  areas  include  one  standard  deviation. 

Develop  optimization  methods  for  weights  and  placements  of  the  light  sources  to  achieve  uniform 
PDT  dose  in  the  prostate  gland 

To  deliver  uniform  photodynamic  therapy  (PDT)  dose  to  the  prostate  gland,  it  is  necessary  to  develop 
algorithms  that  optimize  the  locations  and  strengths  (power  x  exposure  times,  in  units  of  Joule)  of  the 
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light  sources.  Since  tissue  optical  properties  may  change  with  time,  rapid  (almost  real-time)  optimization 
is  desirable.  We  evaluated  the  use  of  the  Cimmino  feasibility  algorithms  for  this  purpose  because  it  is 
fast,  linear,  and  always  converges  reliably.  PDT  for  the  prostate  is  performed  with  cylindrical  diffusing 
fibers  (CDF)  of  various  lengths  inserted  perpendicular  to  a  base  plate  to  obtain  longitudinal  coverage  by  a 
matrix  of  parallel  catheters.  The  standard  plan  for  the  protocol  uses  equal  unit-length  source  strength  and 
equal  spaced  (1-cm)  loading.  We  have  developed  several  search  procedures  to  aid  the  user  in  choosing 
the  positions,  lengths,  and  intensities  of  the  CDFs.  The  Cimmino  algorithm  can  be  used  in  these 
procedures  to  determine  the  strengths  of  the  light  sources  at  each  step  of  the  iterative  selection  process. 
Maximum  and  minimum  bounds  for  allowed  doses  to  the  sample  points  in  four  volumes  (prostate, 
urethra,  rectum,  and  background)  constrained  the  solutions  for  the  strengths  of  the  linear  light  sources. 
Uniform  optical  properties  were  assumed.  Compared  with  the  standard  plan,  our  algorithm  was  shown  to 
quickly  produce  a  plan  in  real  time  that  covered  the  planned  target  volume  (PTV)  and  spared  the  urethra 
and  rectum. 


We  have  studied  three  kinds  of  problems.  (1)  Cimmino  1:  Given  cylindrical  diffusing  fibers  (CDFs)  with 
every  quantity  specified,  namely,  the  number  of  CDFs,  the  template  slots,  the  source  lengths,  and  the 
retractions  of  the  CDF  into  the  slot,  find  the  source  strengths  (emitted  power  multiplied  by  duration  of 
illumination  [J])  to  satisfy  the  prescribed  dose  constraints.  (2)  Cimmino  2  (or  3):  Given  only  the  number 
of  CDFs  and  the  allowed  set  of  template  slots,  find  the  particular  source  slots,  source  parameters,  and 
source  strengths  that  are  optimal.  (3)  Cimmino  4:  Given  the  number  and  locations  of  the  source  slots,  find 
the  optimal  source  parameters  (source  lengths  and  retractions)  and  source  strengths.  The  only  difference 
between  Cimmino  2  and  3  is  the  different  constraints  for  the  rectum. 


Fig.  10.  DVH  comparison  of  the  manual  standard  plan  vs.  the 
Cimmino-based  search  results  for  optical  properties  /ja  =  0.3 
cm'1,  fis'  =14  cm'1  for  (a)  prostate  (b)  urethra  (c)  rectum. 
The  standard  plan  uses  uniform  1-cm  source  loading  with 
uniform-strength.  The  optimized  results  are:  Cimmino  1  uses 
the  same  fixed  source  positions  and  source  parameters  as  the 
standard  plan  but  Cimmino  optimized  weights;  Cimmino  2 
uses  Cimmino  optimized  source  lengths,  loading,  and  template 
for  12  linear  sources  with  upper  constraints  of  300%  for 
rectum  and  urethra;  Cimmino  3  is  the  same  as  Cimmino  2  but 
with  upper  constraints  of  200%  for  rectum  and  urethra; 
Cimmino  4  uses  Cimmino  optimized  source  lengths,  loading 
and  template  for  51  linear  sources.  See  Table  4  for  the  upper 
dose  bounds  used  for  each  Cimmino-search  algorithm 


We  have  compared  the  dose-volume  histograms  (DVHs)  obtained  in  the  standard  plan,  i.e.,  the  case  of 
uniform  source  strength  with  all  source  slots  and  parameters  specified,  which  is  the  present  practice  in  the 
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clinic,  with  those  obtained  in  Cimmino  1,  Cimmino  2,  Cimmino  3,  and  Cimmino  4.  Figure  10  shows  the 
DVHs  for  the  average  optical  properties  (pia  =  0.3  cm'1  and  //s,  ’  =  14  cm'1),  in  which  dosimetric 
improvement  is  observed  in  the  above  order.  Cimmino  1,  which  optimizes  the  source  strength  only,  did 
not  provide  significant  improvement  over  the  standard  plan.  Cimmino  4  which  had  sufficient  number 
(51)  of  CDFs  to  cover  the  entire  prostate  gave  the  best  DVH,  although  it  is  clinically  impractical  to  use  so 
many  sources.  The  DVH  for  Cimmino  4  represents  the  best  possible  mathematical  solution.  Cimmino  2 
and  3,  which  optimized  the  source  locations,  showed  substantial  improvements  over  the  standard  plan,  not 
only  for  the  coverage  of  the  target  (prostate)  but  also  for  the  dose  reduction  of  the  urethra  and  rectum. 
Cimmino  3,  which  used  200%  upper  constraint  for  rectum,  gave  better  results  than  Cimmino  2,  which 
used  300%  upper  constraint  for  rectum. 

Our  study  has  shown  that  significant  sparing  of  the  urethra  and  rectum  could  be  implemented  with  the 
optimization  algorithm.  The  combined  selection  of  positions,  lengths,  and  strengths  of  interstitial  light 
sources  could  improve  the  outcome.  Compared  with  the  standard  plan,  our  algorithm  has  been  shown  to 
quickly  produce  a  plan  in  real  time  that  covered  the  planned  target  volume  (PTV)  and  spared  the  urethra 
and  rectum.  The  clinical  implication  of  our  findings  is  that  changes  in  the  treatment  plan  may  be  made  in 
the  operating  room  before  and  during  PDT  to  optimize  light  delivery.  The  calculation  time  is  always  less 
than  300  s  (see  Table  2  in  Appendix  5,  Ref.  6). 

Key  Research  Accomplishments 

1.  Completion  of  comparison  between  measurement  and  calculation  in  prostate  in-vivo 
assuming  homogeneous  optical  properties. 

2.  Development  of  instrument  for  fluorescence  spectra  interstitially. 

3.  Completion  of  the  optimization  algorithm  for  light  source  weights,  locations,  and 
lengths. 

4.  Development  of  instrument  for  measuring  distribution  of  optical  properties,  MLu 
concentration,  and  St02  in  vivo  and  characterization  of  measurement  uncertainty  of 
the  instrument. 

Reportable  outcomes 

1.  Zhu  TC,  Dimofte  A,  Finlay  JC,  Stripp  D,  Busch  T,  Miles  J,  Whittington  R,  Malkowicz  SB, 
Tochner  Z,  Glatstein  E,  Hahn  SM,  “Optical  properties  of  Human  Prostate  at  732nm  Measured  in 
vivo  during  motexafin  lutetium-mediated  photodynamic  therapy,”  Photochem  photobiol.  81 :  96- 
105(2005). 

2.  Zhu  TC,  Finlay  JF,  Hahn  SM,  “Determination  of  the  distribution  of  light,  optical  properties, 
drug  concentration,  and  tissue  oxygenation  in-vivo  in  human  prostate  during  motexafin 
lutetium-mediated  photodynamic  therapy”.  J.  Photochem  Photobiol  B.  79:231-241  (2005) 

3.  Dimofte  A,  Finlay  JC,  Zhu  TC,  “A  method  for  determination  of  the  absorption  and  scattering 
properties  interstitially  in  turbid  media,”  Phys  Med  Biol  50  2291-2311  (2005). 

4.  Finlay  JC,  Zhu  TC,  Dimofte  A,  Stripp  D,  Malkowicz  SB,  Whittington  R,  Miles  J,  Glatstein 
E,  Hahn  SM,  “In  vivo  measurement  of  fluorescence  emission  in  the  human  prostate  during 
photodynamic  therapy”,  Proc  SPIE  5689:299-310  (2005) 

5.  Altschuler  MD,  Zhu  TC,  Jun  Li,  Hahn  SM,  “Optimization  of  light  sources  for  prostate 
photodynamic  therapy,”  Proc.  SPIE  5689:186-197  (2005). 
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6.  Zhu  TC,  Dimofte  A,  Finlay  JC,  Glatstein  E,  Hahn  SM,  “Detector  Calibration  factor  for 

interstitial  in-vivo  light  dosimetry  using  isotropic  detectors  with  scattering  tip,”  Proc.  SPIE 

5689:174-185  (2005). 

Conclusions: 

We  have  completed  specific  aim  2  of  the  project,  i.e.,  light  dose  calculation  and  optimization, 
characterization  of  photosensitizer  drug  concentration,  and  accumulate  clinical  data  in  vivo  in 
human  prostate  at  732  nm.  We  have  improved  the  motorized  probe  and  associate  software  to 
measure  the  distribution  of  optical  properties,  drug  concentration,  and  StC>2  interstitially.  We 
have  developed  a  fluorescence  technique  to  measure  in-vivo  phototosensitizer  distribution  interstitially  in 
5  patients  and  have  verified  its  accuracy  using  the  absorption  spectra.  We  have  also  developed  software 
to  automatically  optimize  the  light  source  weights,  lengths,  strengths  in  near  real-time  to  improve  the  light 
fluence  rate  distribution  in  prostate  and  spare  dose  to  critical  structures.  We  have  compared  the 
measurement  and  calculation  of  in-vivo  light  fluence  rate  in  14  patients  assuming  homogeneous  optical 
properties  in  prostate.  They  agreed  within  a  standard  uncertainty  of  59.5%.  Our  measurements 
demonstrate  heterogeneity  in  both  the  optical  properties  and  the  photosensitizer  drug  distribution.  Thus, 
we  concluded  that  it  is  necessary  to  determine  the  3D  distribution  of  optical  properties  and  drug 
concentration  in  order  to  optimize  the  PDT  dose  distribution  for  each  individual  patient. 
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ABSTRACT 

Characterization  of  the  tissue  light  penetration  in  prostate 
photodynamic  therapy  (PDT)  is  important  to  plan  the 
arrangement  and  weighting  of  light  sources  so  that  sufficient 
light  fluence  is  delivered  to  the  treatment  volume.  The  optical 
properties  (absorption  [pa],  transport  scattering  [p^  and 
effective  attenuation  [pefr]  coefficients)  of  13  patients  with 
locally  recurrent  prostate  cancer  were  measured  in  situ  using 
interstitial  isotropic  detectors.  Measurements  were  made  at 
732  nm  before  and  after  motexafin  lutetium  (MLu)-mediated 
PDT  in  four  quadrants.  Optical  properties  were  derived  by 
applying  the  diffusion  theory  to  the  fluence  rates  measured  at 
several  distances  (0.5-5  cm)  from  a  point  source.  pa  and  p' 
varied  between  0.07  and  1.62  cm-1  (mean  0.37  ±  0.24  cm-1) 
and  1.1  and  44  cm-1  (mean  14  ±  11  cm-1),  respectively. 
pa  was  proportional  to  the  concentration  of  MLu  measured  by 
an  ex  vivo  fluorescence  assay.  We  have  observed,  on  average, 
a  reduction  of  the  MLu  concentration  after  PDT,  presumably 
due  to  the  PDT  consumption  of  MLu.  pe(T  varied  between  0.91 
and  6.7  cm-1  (mean  2.9  ±  0.7  cm-1),  corresponding  to  an 
optical  penetration  depth  (8  =  l/pc(T)  of  0.1-1.1  cm  (mean  0.4  ± 
0.1  cm).  The  mean  penetration  depth  at  732  nm  in  human 
prostate  is  at  least  two  times  smaller  than  that  found  in  normal 
canine  prostates,  which  can  be  explained  by  a  four  times 
increase  of  the  mean  value  Of  p'  in  human  prostates.  The  mean 
light  fluence  rate  per  unit  source  strength  at  0.5  cm  from 
a  point  source  was  1.5  ±  1.1  cm-2,  excluding  situations  when 
bleeding  occurs.  The  total  number  of  measurements  was  N  = 
121  for  all  mean  quantities  listed  above.  This  study  showed 
significant  inter-  and  intraprostatic  differences  in  the  optical 
properties,  suggesting  that  a  real-time  dosimetry  measurement 
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and  feedback  system  for  monitoring  light  fluences  during 
treatment  should  be  considered  for  future  PDT  studies. 


INTRODUCTION 

Photodynamic  therapy  (PDT)  is  a  treatment  modality  using  light  of 
an  appropriate  wavelength  in  the  presence  of  oxygen  to  activate 
a  photosensitizing  drug,  which  then  causes  localized  cell  death  or 
tissue  necrosis.  Using  a  surface  illumination  technique,  PDT  has 
been  used  to  treat  superficial  tumors,  including  those  of  the  skin, 
lung,  esophagus  and  bladder  (1).  This  technique  is,  however, 
inadequate  when  applied  to  large  bulky  tumors  or  solid  organs 
because  of  limited  light  penetration  into  tissue.  Interstitial  light 
delivery,  wherein  optical  fibers  are  placed  directly  into  bulky  tumors 
or  organs,  is  a  more  rational  approach  in  these  cases. 

The  prostate  gland  is  an  organ  that  appears  to  be  a  good  target  for 
interstitial  PDT.  Tumors  of  the  prostate  are  often  confined  to  the 
prostate  itself,  and  brachytherapy  techniques  used  for  the  placement 
of  radioactive  seed  implants  can  be  adapted  for  the  placement  of 
interstitial  optical  fibers  (2).  Several  preclinical  studies  have  evalu¬ 
ated  the  feasibility  of  delivering  PDT  to  the  prostate  via  this 
interstitial  approach  (3-8).  A  trial  of  interstitial  prostate  PDT  in 
humans  has  been  reported  by  Nathan  et  al.  (9),  who  treated  14  men 
with  locally  recurrent  prostate  cancer  using  m<?,so-tetrahydroxy- 
phenyl  chlorin-mediated  interstitial  PDT.  The  light  treatment  was 
directed  toward  regions  from  which  biopsies  showed  cancer  or 
which  were  suspicious  on  imaging  studies.  Because  prostate  cancer 
is  a  multifocal  disease,  our  PDT  protocol  is  designed  to  ablate  the 
prostate  gland  completely.  The  development  of  a  light  delivery 
technique  appropriate  to  this  goal  has  necessitated  an  improved 
understanding  of  light  dosimetry,  which  is  critical  in  planning  the 
configuration  of  multiple  fibers  within  the  organ  or  tumor. 

Several  investigators  have  attempted  to  characterize  the  optical 
properties  of  prostate  tissue  in  animals  (10-12)  and  in  humans  (13- 
15)  to  predict  light  dosimetry  more  reliably.  Using  diffusion  theory 
for  a  point  source,  the  absorption  (pa)  and  transport  scattering  (p's) 
coefficients  of  a  particular  tissue  can  be  determined  yielding  the 
effective  attenuation  coefficient  (peff)>  which  provides  a  measure  of 
light  penetration  in  that  tissue  (11).  This  measurement  is  a  critical 
factor  in  planning  interstitial  light  source  placement.  The  optical 
properties  of  prostate  tissue  have  been  most  extensively  evaluated  at 
a  wavelength  of  630  nm  because  it  is  used  to  activate  the  most  widely 
used  photosensitizer,  Photofrin.  However,  it  has  been  suggested  that 
a  wavelength  of  630  nm  is  suboptimal  in  achieving  uniform  and 
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Figure  1.  (a)  Transrectal  ultrasound  image  of  a  human  prostate,  showing  the 
position  of  source  (closed  circle)  and  detector  (*)  fibers.  The  source  positions 
labeled  1 ,  2,  3  and  4  were  used  for  optical  properties  measurement  for  the 
RUQ,  LUQ,  RLQ  and  LLQ,  respectively,  with  the  detector  placed  in 
a  separate  catheter  in  each  quadrant.  The  open  circle  on  the  right  upper 
quadrant  is  for  a  linear  source  that  passed  through  the  position  but  is  too  short 
to  have  active  light  component  in  the  cross-section  plane.  The  grid  on  the 
template  (“+”)  is  0.5  cm  apart,  (b)  Schematic  of  the  measurement  geometry, 
illustrating  the  coordinates  used  to  determine  the  source-detector  distance. 


adequate  penetration  due  to  the  absorption  of  light  by  naturally 
occurring  chromophores,  such  as  hemoglobin,  melanins,  cyto¬ 
chromes  and  flavoproteins,  at  that  wavelength  (10,15).  In  addition, 
changing  concentrations  of  hemoglobin  (in  particular  deoxyhemo¬ 
globin)  during  PDT  may  also  alter  the  optical  properties  of  tissues 
(11).  As  a  result,  Amfield  et  al.  (10)  investigated  the  use  of  longer 
wavelengths  of  light  beyond  the  absorption  bands  of  these  naturally 
occurring  molecules.  A  comparison  between  630  and  789  nm  light 
revealed  a  consistent  and  substantial  increase  in  light  penetration 
with  the  789  nm  light. 

Based  on  a  feasibility  and  toxicity  study  in  a  canine  prostate  model 
(7),  we  have  started  a  Phase  I  study  of  motexafin  lutetium  (MLu)- 
mediated  PDT  for  prostate  cancer.  MLu  is  a  water-soluble,  second- 
generation  synthetic  photoactive  dmg  that  has  a  Q-band  absorption 
peak  at  732  nm  (16,17).  In  this  study  we  evaluated  the  optical  pro¬ 
perties  of  the  human  prostate  by  illuminating  it  with  a  point  source 
emitting  732  nm  light.  Differences  in  optical  properties  were  com- 


Table  1.  Distribution  of  patients  for  measurements  of  optical  properties. 
The  number  in  parentheses  is  the  number  of  valid  measurements  that  is  not 
affected  by  bleeding 


MLu  Dmg-  Light 

Patient  dose  treatment  fluence 
number  (mg/kg)  interval  (h)  (J/cm2) 

No. 

measured 
before  PDT 

No. 

measured 
after  PDT 

Optical 
properties 
along  a 
catheter? 

2 

0.5 

24 

25 

4(1) 

4(1) 

No 

3 

0.5 

24 

25 

4(1) 

1  (1) 

No 

4 

1 

24 

25 

4(3) 

1  U) 

No 

5 

1 

24 

25 

4(1) 

1  (1) 

No 

6 

1 

24 

25 

3(2) 

3(2) 

No 

7 

1 

6 

25 

4(0) 

4(2) 

No 

8 

2 

6 

25 

4(0) 

4(0) 

No 

9* 

2 

6 

50 

7(6) 

7(7) 

No 

10* 

2 

6 

100 

8(6) 

7(5) 

No 

11* 

2 

3 

25 

8(7) 

7(5) 

No 

12* 

2 

3 

50 

21  (20) 

9(7) 

Yes 

13* 

2 

3 

100 

28  (20) 

20  (20) 

Yes 

14* 

2 

3 

100 

29  (2) 

29  (0) 

Yes 

*Measurements  made  with  a  motorized  probe. 


pared  before  and  after  PDT  in  which  the  entire  prostate  was  treated. 
Intraprostatic  as  well  as  interprostatic  differences  were  evaluated. 

MATERIALS  AND  METHODS 

Patient  selection,  surgical  and  PDT  procedure.  A  Phase  1  clinical  trial  of 
MLu-mediated  PDT  in  patients  with  locally  recurrent  prostate  carcinoma 
was  initiated  at  the  University  of  Pennsylvania.  The  protocol  was  approved 
by  the  Institutional  Review  board  of  the  University  of  Pennsylvania,  the 
Clinical  Trials  and  Scientific  Monitoring  Committee  of  the  University  of 
Pennsylvania  Cancer  Center  and  the  Cancer  Therapy  Evaluation  Program  of 
the  National  Cancer  Institute.  A  total  of  14  patients  were  treated,  of  which 
13  patients  have  undergone  measurement  of  optical  properties.  Each  patient 
who  signed  the  informed  consent  document  underwent  an  evaluation, 
which  included  magnetic  resonance  imaging  of  the  prostate,  bone  scan, 
laboratory  studies  including  prostatic  specific  antigen  and  a  urological 
evaluation.  Approximately  2  weeks  before  the  scheduled  treatment, 
a  transrectal  ultrasound  (TRUS)  was  performed  for  treatment  planning. 
An  urologist  drew  the  target  volume  (the  prostate)  on  each  slice  of  the 
ultrasound  images.  These  images  were  spaced  0.5  cm  apart  and  were 
scanned  with  the  same  ultrasound  unit  used  for  treatment.  A  built-in 
template  with  a  0.5  cm  grid  projected  the  locations  of  possible  light  sources 
relative  to  the  prostate.  A  treatment  plan  was  then  prepared  to  determine  the 
location  and  length  of  light  sources.  Cylindrical  diffusing  fibers  (CDF)  with 
active  lengths  of  1 ,  2,  3, 4  and  5  cm  were  used  as  light  sources.  The  sources 
were  spaced  1  cm  apart,  and  the  light  power  per  unit  length  was  less  than  or 
equal  to  150  mW/cm2  for  all  optica]  fibers.  The  length  of  the  CDF  at 
a  particular  position  within  the  prostate  was  selected  to  cover  the  full  length 
of  the  prostate  (see  Fig.  la).  The  final  plan  often  required  that  the  prostate 
be  divided  into  four  quadrants.  Four  isotropic  detectors  were  used,  each 
placed  in  the  center  of  one  quadrant.  A  fifth  isotropic  detector  was  placed  in 
a  urethral  catheter  to  monitor  the  light  fluence  in  the  urethra. 

MLu  was  administered  intravenously  to  patients  at  various  times  before 
light  delivery  (depending  on  the  dose  level).  In  the  initial  six  patients,  MLu 
was  delivered  24  h  before  light  delivery.  The  drug-treatment  interval  and 
dose  of  MLu  were  subsequently  changed  according  to  a  predetermined 
scheme  dictated  by  the  protocol.  The  drug-treatment  interval  and  the 
amount  of  MLu  for  the  latest  patients  were  3  h  and  2  mg/kg,  respectively 
(see  Table  1). 

The  patients  were  anesthetized  in  the  operating  room  with  general 
anesthesia  to  minimize  patient  movement  during  the  procedure.  TRUS- 
guided  biopsies  for  MLu  measurements  were  obtained  before  light  delivery. 
The  ultrasound  unit  was  used  to  guide  needle  placement  in  the  operating 
room.  A  template  was  attached  to  the  ultrasound  unit  and  was  matched  to 
the  same  5  mm  grid  used  for  treatment  planning.  Four  detector  catheters 
(one  for  each  quadrant)  were  inserted  into  the  prostate.  These  detectors  were 
kept  in  place  during  the  entire  procedure  of  PDT  treatment.  Four  additional 
preplanned  treatment  catheters  for  light  sources  were  then  inserted  0.5  or 
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0.7  cm  away  from  the  detector  catheters.  These  source  catheters  were  used 
for  light  delivery  and  optical  properties  measurements.  The  optical 
properties  of  the  prostate  in  each  quadrant  were  then  measured  before 
light  delivery  using  these  existing  detectors  and  a  point  source  inserted  into 
one  source  catheter.  The  time  for  optical  property  measurements  was 
approximately  5  min.  The  PDT  treatment  was  then  performed  one  quadrant 
at  a  time  by  inserting  the  CDF  into  the  source  catheters.  A  15  W  diode  laser, 
model  730  (Diomed,  Ltd.,  Cambridge,  UK)  was  used  as  the  732  nm  light 
source.  The  treatment  time  for  each  quadrant  was  dependent  on  the  detector 
reading  in  that  quadrant.  The  light  was  delivered  to  ensure  that  the 
cumulative  fluence  reading  of  all  four  isotropic  detectors  reached  the  same 
prescription  fluence.  The  total  treatment  time  varied  but  was  always  less 
than  1  h.  After  light  delivery,  the  optical  properties  of  the  prostate  in  all  four 
quadrants  were  measured  again  using  a  point  source  (for  another  5  min). 
The  light  sources  and  detectors  were  then  removed,  and  posttreatment 
biopsies  were  performed. 

Diffusion  theory  and  optimization  algorithm,  p'  and  pa  coefficients 
characterize  the  scattering  and  absorption  properties  of  tissues.  With  the 
diffusion  approximation,  the  light  fluence  rate  (|i  at  a  distance  r  from  a  point 
source  of  source  strength,  S,  can  be  expressed  as  (18) 

9  4tc /-p.  An r  K) 

where  S  is  the  source  strength  of  the  point  source,  <(>(/•)  is  the  fluence  rate  at 
position  r,  Petr  =  is  the  Fetr  in  tissues  and  is  applicable  for  a  wider 

range  of  pa  and  pa  combinations  than  the  traditional  definition  of  pcff  = 
i/3pa(p'  +  pa)  (19),  although  it  is  obtained  as  an  approximation  of  the 
latter  expression,  r  =  y/x1  +  h2,  where  x  and  h  are  the  parallel  and 
perpendicular  distances  from  the  center  of  the  point  source  (see  Fig.  1).  The 
two  free  parameters  (pa  and  p')  are  inherently  separable  in  that  for  a  point 
source  with  a  given  source  strength,  the  magnitude  of  the  fluence  rate  near 
the  light  source  (x  =  0)  is  determined  by  p'  only  and  the  slope  of  the  spatial 
decay  of  the  light  fluence  rate  is  determined  by  peft  only. 

In  theory,  measurements  of  <J>  at  two  different  distances  r  from  the  point 
source  are  sufficient  to  determine  both  pa  and  p' .  In  reality,  our  measurements 
contain  at  least  four  different  distances  and  up  to  800  distances  for  latter 
measurements.  Measurements  at  multiple  sites  allow  evaluation  of  the 
variation  of  these  optical  characteristics  within  the  prostate  volume.  Because 
Eq.  (1)  is  a  nonlinear  equation  of  two  free  parameters  pa  and  pa,  we  used 
a  differential  evolution  algorithm  developed  by  Stom  and  Price  (20).  This 
algorithm  is  simple  and  robust,  and  converges  faster  than  adaptive  simulated 
annealing  or  the  annealed  Nelder  and  Mead  approach  (20).  We  modified  the 
algorithm  to  require  that  all  free  parameters  (pa,  p's  and  h )  are  positive  (21). 
The  deviations  between  measurement  and  fit  are  represented  by  standard 

deviation,  y  ((measurement,  —  fit,)/fit,)2/(n  —  1),  where  n  is  the 

number  of  data  points  needed  in  the  fit. 

Calibration  of  the  isotropic  detector  and  phantom  verification.  Isotropic 
detectors  described  by  Marijnissen  and  Star  (22)  were  used  for  the  study. 
Each  isotropic  detector  was  made  of  an  optical  fiber  with  a  spherical  tip 
made  of  Ti02  (a  scattering  material).  The  isotropic  detectors  were  made  by 
Rare  Earth  Medical  (now  CardioFocus,  Norton,  MA)  and  have  an  isotropy 
of  better  than  ±30%  from  any  direction  except  for  angles  within  30°  of  the 
optical  fiber  attachment  point.  Each  detector  fiber  was  connected  to 
a  photodiode  via  an  SMA  connector.  The  measured  photovoltage  (V)  from 
the  isotropic  detector  was  converted  to  light  fluence  rate  using: 

<)>  =  a(V  -  b),  (2) 

where  a  was  the  conversion  factor  and  b  characterized  the  leakage  of  the 
photodiode.  This  calibration  was  performed  under  collimated  732  nm  laser 
light  in  air.  When  the  isotropic  detector  was  used  to  measure  fluence  rate  in 
tissue  (or  in  a  liquid  optical  phantom),  a  tissue  (or  water)  correction  factor 
of  2.0  was  used.  This  value  was  determined  by  measuring  the  response  of 
the  isotropic  detectors  in  and  out  of  water  medium  for  the  same  incident 
irradiance  (22).  The  calibration  of  the  isotropic  detectors  was  checked  to  be 
accurate  to  within  5%  using  an  integrating  sphere  before  each  individual 
measurement. 

A  homemade,  12-channel  light  dosimetry  system  was  used  for  all  in  situ 
measurements.  Five  different  isotropic  detectors  were  used.  The  conversion 
factors  were  a  =  129-161  and  a  =  63-72  mW/cm2/V  for  the  four  isotropic 
detectors  with  0.5  mm  scattering  tips  and  the  isotropic  detector  with  1  mm 
scattering  tip  used  in  the  urethra,  respectively,  b  =  0.020  V  for  all  isotropic 


detectors.  The  response  was  found  to  be  linear  over  the  light  fluence  rate 
range  (0-1200  mW/cm2)  to  within  ±5%.  The  minimum  resolution  of  the 
dosimetry  system  was  better  than  0.1  mW/cm2  after  the  introduction  of  an 
additional  software-controllable  analog  gain  to  adjust  the  dynamic  range  of 
the  light  dosimetry  system. 

To  ensure  that  the  optimization  technique  was  robust  and  to  estimate  the 
accuracy  of  the  optimization  algorithm  to  extrapolate  optical  properties,  we 
measured  the  optical  properties  of  an  optical  phantom  with  known  pa  and  pa 
using  the  same  methods  as  those  used  in  human  patients.  The  optical 
phantom  was  composed  of  pure  scattering  medium  (liposyn  III,  30%  Abbott 
Lab,  North  Chicago,  IL)  and  pure  absorbing  medium  (green  ink).  This  type 
of  phantom  has  been  described  by  Flock  et  al.  (23).  The  concentrations  of 
liposyn  used  were  0.25%  and  0.5%  by  volume  in  water.  Three  green  ink 
concentrations  were  used:  0%,  0.014%  and  0.07%.  The  optical  properties 
(p'  and  pj  of  the  liquid  phantom  were  independently  characterized  by 
analyzing  the  decay  of  intensity  along  the  axis  of  a  broad  collimated  light 
beam  illuminating  the  phantom  from  above  (12).  They  were  p'  =  3.6  and 
7.2  cm-1  and  pa  =  0.02,  0.1  and  0.5  cm-1.  The  phantom  was  placed  in 
a  plastic  container,  which  was  painted  black  and  was  large  enough  (18.2  X 
14.6  X  7.7  cm3)  to  protect  the  detector  from  boundary  effects.  A  point  light 
source  with  an  outside  diameter  of  500  pm  (Pioneer  Optics,  Windsor  Locks, 
CT)  was  used  for  light  delivery,  along  with  a  polyurethane  template  and 
plastic  catheters.  Fluence  rates  were  measured  at  various  positions  along  the 
length  of  detector  catheters  placed  at  fixed  distances  (h  —  0.5-0.7  cm)  from 
the  source  catheter. 

A  second  set  of  phantoms  was  designed  to  verify  that  the  pa  of  the  phantom 
varies  linearly  with  MLu  concentration.  These  phantoms  were  constructed 
with  a  lipid  content  of  0.5%.  (Note  a  different  batch  of  Liposyn  was  used.  As 
a  result,  |4  =  4  cm-1  for  0.5%  lipid  content.)  The  measurements  and  data 
analysis  proceeded  as  above  using  732  nm  laser  and  the  same  setup. 

The  transmission  of  the  light  from  the  sources  through  the  transparent 
catheters  was  measured  in  a  6  inch  diameter  integrating  sphere  (IS60, 
LabSphere,  Inc.,  North  Sutton,  NH)  as  a  ratio  of  the  detector  reading  with 
and  without  the  catheter.  The  result  was  1 .0  for  both  point  and  linear  sources. 
The  transmission  of  the  catheter  for  an  isotropic  detector  in  air  was  between 
0.95  and  1  depending  on  the  whether  the  light  source  was  farther  away  or 
very  close  to  the  detector.  The  latter  was  expected  in  tissue  because  the  light 
was  coming  from  all  directions  and  was  very  close  to  the  detector  in  the 
catheter.  Based  on  these  results  no  transmission  correction  was  necessary  for 
the  transparent  catheters  for  either  the  light  source  or  the  isotropic  detectors. 

PDT  and  in  vivo  light  measurements.  A  total  of  14  patients  were  treated. 
MLu  (Pharmacyclics,  Inc.,  Sunnyvale,  CA)  injection  was  administered 
intravenously  (0.5—2  mg/kg)  3-24  h  before  light  administration  (16,24). 
After  the  patient  was  anesthetized,  the  interstitial  CDF  were  placed  in  the 
gland  using  a  template  with  evenly  spaced  holes,  which  was  attached  to  the 
TRUS  unit.  A  17-gauge  plastic  catheter  (Flexi-needle,  Best  Industries,  Inc., 
Springfield,  VA)  containing  a  metal  trocar  was  placed  through  the  template 
and  into  the  prostate.  The  trocar  was  removed  and  replaced  with  the  light 
diffuser.  The  light  energy  delivered  was  prescribed  based  on  in  situ 
measured  light  fluence.  Each  patient  received  a  light  fluence  between  25 
and  100  J/cm2,  determined  by  in  vivo  measurement  using  isotropic 
detectors.  The  maximum  unit  length  source  strength  in  any  one  fiber  was 
limited  to  150  mW/cm.  Table  1  summarizes  the  treatment  delivered  to 
all  patients. 

Measurements  were  taken  at  various  distances  from  the  light  source  before 
and  after  light  treatment  in  various  quadrants  of  each  prostate,  measured  at  2 
cm  from  the  apex  of  the  prostate.  For  the  last  two  patients,  the  light  sources 
were  moved  along  the  catheter  in  several  locations  to  also  quantify  the 
variation  of  optical  properties  in  the  prostate  gland  along  the  catheter. 

Ex  vivo  MLu  measurement.  In  each  patient  a  needle  biopsy  of  prostate 
tissue  from  each  quadrant  was  collected  before  and  after  light  delivery. 
Biopsies  were  immediately  frozen  on  dry  ice,  transported  to  the  laboratory 
and  stored  at  — 80°C  until  the  time  of  assay.  Caution  was  taken  to  protect  the 
samples  from  light.  MLu  was  extracted  from  the  needle  biopsies  using  a 
procedure  based  on  a  previous  report  (25).  Tissue  samples  were  thawed  to 
room  temperature,  weighed  and,  depending  on  the  amount  of  tissue 
available,  two  to  four  of  the  needle  biopsies  were  combined  for  analysis. 
Only  biopsies  collected  at  the  same  time  point,  i.e.  either  pre-  or  post-PDT, 
were  combined.  A  combined  tissue  weight  of  —5-25  mg  was  sought. 
Specimens  were  cut  into  shorter  lengths,  placed  in  a  2  mL  capped 
polypropylene  tube  and  homogenized  (Polytron  1200)  in  400  pL  of  phos¬ 
phate  buffer  (24  m M,  pH  7.5).  Homogenates  were  mixed  with  400  pL  of 
chloroform,  then  400  pL  of  methanol  was  added.  After  centrifugation  (3500 
rpm,  15  min)  the  organic  layer  was  collected  and  200  pL  was  transferred  to 
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Figure  2.  Measurement  of  light  fluence  rate  (mW/cm2)  at  different 
distances,  x,  from  a  point  source  in  liquid  optical  phantoms:  (a)  p'  =  3.6 
cm-1  and  pa  =  0.02, 0. 1  and  0.5  cm-1  and  (b)  p'  =  7.2  cm-1 ,  pa  =  0.02,  0.08 
and  0.40  cm-1.  Symbols  represent  measurements  with  an  isotropic  detector. 
The  solid  lines  are  the  best  fit  with  the  resulting  fit  optical  properties  (from 
top  to  bottom):  (a)  p'  =  3.7,  3.8  and  3.4  cm-1  and  pa  =  0.03, 0.097  and  0.53 
cm-1  and  (b)  pa  =  7.5,  7.1  and  7.3  cm-1,  pa  =  0.022,  0.074  and  0.38  cm-1, 
respectively. 


a  cuvette.  The  fluorescence  of  the  homogenized  sample  was  measured  by 
a  spectrofluorometer  (FluoroMax-3,  Jobin  Yvon,  Inc.,  Edison,  NJ)  with  Xox 
of  474  nm  and  X,cm  of  740  nm  (emission  scan  range  from  650  to  850  nm). 
MLu  concentration  in  the  tissue  was  calculated  on  the  basis  of  the  increase  in 
fluorescence  resulting  from  the  addition  of  a  known  amount  of  MLu  to  each 
sample  after  its  initial  reading.  Data  are  presented  as  nanograms  of  MLu  per 
milligram  of  tissue. 

Correlation  between  pa  and  MLu  concentration.  To  correlate  the  rela¬ 
tionship  between  pa  measured  in  vivo  with  the  MLu  measured  from  ex  vivo 
biopsy,  we  plot  out  the  subset  of  pa  data  (24  points)  that  had  corresponding  ex 
vivo  biopsies.  Notice  that  although  the  ex  vivo  biopsies  were  taken  from  the 
same  prostate,  they  are  not  necessarily  at  the  same  location  as  that  used  for 
pa  measurement  because  the  prostate  template  was  not  used  for  the  biopsy. 
Once  the  parameters  for  a  linear  relationship  between  pa  and  MLu  concentra¬ 


tion  are  established,  all  measured  pa  data  (121  points)  were  used  to  examine 
the  variation  of  MLu  concentration  intra-  and  interprostate. 


RESULTS 

Figure  2  compares  the  fit  and  measurement  of  light  fluence  rate  at 
different  distances  from  a  point  source  in  two  optical  phantoms 
made  of  different  concentrations  of  liposyn  and  green  ink  mixtures. 
The  solid  lines  are  the  result  of  nonlinear  fitting  using  the  differential 
evolution  algorithm.  The  resulting  optical  properties  are  shown.  The 
extrapolated  p'  agreed  with  the  true  values  to  within  ±6%;  the 
extrapolated  pa  agreed  with  the  true  values  to  within  ±11%. 

Figure  3  shows  the  variation  of  measured  light  fluence  rate 
distribution  (solid  lines)  and  associated  fit  (symbols)  for  (1)  different 
locations  in  the  same  prostate  gland  and  (2)  the  same  location  in 
a  prostate  before  and  after  light  delivery.  In  Fig.  3a,  an  example  of 
the  light  fluence  rate  radial  distribution  for  the  right  lower  (RLQ), 
right  upper  (RUQ)  and  left  upper  (LUQ)  quadrants  of  one  prostate  is 
shown.  In  this  particular  patient,  the  extrapolated  optical  properties 
were  very  different  between  RLQ  (pa  =  0.23  cm-1,  p'  =  7.3  cm-1) 
and  RUQ  (pa  =  0.44  cm-1,  p'  =  12  cm-1).  A  large  variation  was  also 
observed  among  different  patients  (see  Fig.  5).  When  evaluated 
across  all  patients,  the  extrapolated  attenuation  coefficients  varied 
between  0.9  and  6.7  cm-1,  whereas  the  extrapolated  p'  varied 
between  1 . 1  and  44  cm-1 .  As  a  result  of  the  heterogeneity  of  optical 
properties,  the  light  fluence  rates  per  unit  source  strength  at  0.5  cm 
from  the  point  source  varied  between  0.2  and  4.1  cm-2  among  the 
patients.  However,  no  significant  change  of  optical  properties  was 
observed  before  and  after  PDT  treatment  in  the  same  site  for  most 
but  not  all  patients  (Fig.  3b). 

Figure  4  shows  the  temporal  dependence  of  the  light  fluence  rate 
at  the  center  of  the  four  treatment  quadrants.  The  large  variation  of 
light  fluence  rate  was  due  to  the  movement  of  the  isotropic  detector. 
Notice  that  we  did  not  observe  any  significant  light  fluence  rate  in 
urethra  in  all  patients  treated  so  far. 

The  overall  variation  of  optical  properties  (pa,  p'  and  peff)  and 
the  optical  penetration  depth  (5  =  l/pefF)  before  (first  bar)  and  after 
(second  bar)  PDT  for  each  patient  studied  is  shown  in  Fig.  5.  The 
error  bars  represent  the  standard  deviation  among  the  four  different 
quadrants  in  the  same  patient.  No  error  bar  is  present  in  cases 
where  only  one  measurement  is  available.  These  results  are  also 
summarized  in  Table  2. 

The  interprostatic  heterogeneity  of  light  fluence  rate  at  0.5  cm 
from  the  point  source  per  unit  source  strength  is  shown  in  Fig.  6. 
The  x  axis  identifies  each  patient  measured.  The  first  bar  of  the  pair 
was  measured  before  PDT,  and  the  second  bar  of  the  pair  was 
measured  after  PDT.  The  error  bars  represent  the  standard 
deviation  among  the  four  different  quadrants  in  the  same  patient. 
The  solid  line  is  an  average  with  the  range  (dashed  lines) 
corresponding  to  standard  deviation.  The  mean  fluence  rate  per  unit 
source  strength  at  0.5  cm  for  all  sites  (excluding  bleeding)  was 
1.5  ±  1.1  cm-2  (sample  size  121). 

The  relationship  between  the  measured  pa,  and  the  tissue 
concentration  of  MLu,  measured  ex  vivo  after  drug  extraction  from 
biopsy  specimens  is  shown  in  Fig.  7a.  pa  were  measured  in  situ  and 
biopsies  were  collected  before  (O)  and  after  (*)  PDT.  For 
comparison,  we  have  also  plotted  the  measured  pa  of  MLu  vs 
MLu  concentration  (+)  in  a  pure  liposyn  phantom.  Figure  7b  shows 
the  relationship  between  MLu  concentration  measured  from  tissue 
biopsies  (bars)  and  extrapolated  from  pa  (symbols)  vs  the  injected 
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Figure  4.  Light  fluence  rate  vs  time  measured  by  the  in  vivo  dosimetry 
system  for  the  four  quadrants  for  Patient  7  at  732  nm  wavelength.  Large 
disturbances  of  light  fluence  rate  were  caused  by  movement  of  detector 
positions,  per  observation  by  the  operators. 
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Figure  3.  Measured  light  fluence  rate  per  unit  source  strength  (<J>/S)  at 
distances  along  the  catheter,  x,  from  the  point  source  measured  in  vivo  in 
human  prostate  gland  for  Patient  13.  Lines  are  measured  data  and  symbols 
are  fits.  (There  are  too  many  measured  points  to  express  the  measured  data 
clearly  as  symbols.)  (a)  Light  fluence  rates  in  the  right  lower  quadrant  (O), 
right  upper  quadrant  (a)  and  left  upper  quadrant  (*)  of  the  same  prostate 
before  PDT.  (b)  Light  fluence  rates  before  (O)  and  after  (+)  light  treatment 
in  the  right  lower  quadrant  of  the  prostate  gland.  The  optical  properties  are 
O — pa  =  0.23  cnT1,  p'  =  7.3  cm-1,  4>(0.5)/S  =  1.1  cm-2  and  h  =  0.5;  x — 
pa  =  0.44  cm-1,  p'  =  12.0  cm-1,  <j)(0.5 )/S  =  0.78  cm-2  and  h  =  0.7;  * — pa  = 
0.25  cm-1,  p'  =  11.6  cm"1,  (j>(0.5)/S  =  1.3  cm-2  and  h  =  0.7;  4 — pa  =  0.25 
cm-1,  p'  =  6.6  cm-1,  4>(0.5)/S  =  1.0  cm-2  and  h  =  0.5  cm. 


MLu  dosage  to  patients.  Drug  concentration  was  calculated  from 
pa  using  the  linear  relationship  determined  in  Fig.  7a:  c  (ng/mg)  = 
(pa  (cm-1)  —  0.227)/0.0658.  Error  bars  represent  the  standard 
deviation  among  measurements  in  different  locations  within 
a  prostate  or  different  prostates.  The  correlation  coefficient  (R2  = 
0.63)  of  the  linear  relationship  between  pa  and  tissue  concentration 
of  MLu  suggests  good  association,  although  the  exact  numerical 
values  for  the  linear  equation  may  contain  error  due  to  the  spread  of 
the  data. 


DISCUSSION 

The  main  objective  of  this  study  was  to  evaluate  the  optical 
properties  of  prostate  tissue  in  the  human  prostate  using  732  nm 
light.  Comparisons  were  made  before,  during  and  after  MLu- 
mediated  PDT.  Differences  within  and  among  patients  were  also 
studied. 

Intrapatient  heterogeneity  of  optical  properties  was  observed. 
Figure  3  is  a  representative  example  of  this  heterogeneity  at  732 
nm.  There  was  a  difference  in  the  pa  (0.23-0.4  cm-1)  andVs  (6.6— 
12  cm-1)  values  at  different  locations  within  the  same  prostate 
gland  (Fig.  3a).  As  a  result,  the  fluence  rate  per  source  strength  at 
0.5  cm  from  the  point  source  varied  from  0.8  to  1.3  cm-2.  In 
contrast,  the  optical  properties  (pa  and  p')  did  not  change 
substantially  before  or  after  PDT  light  delivery.  Thus,  the  light 
fluence  rate  at  the  same  site  of  the  same  prostate  did  not  vary  much 
before  and  after  PDT  light  delivery  (Fig.  3b). 

In  general,  we  found  a  lack  of  change  in  the  optical  properties  in 
a  particular  site  within  the  prostate  before  and  after  PDT.  This  was 
confirmed  by  the  absence  of  any  time  dependence  of  the  fluence 
rate  in  the  middle  of  the  four  quadrants  (Fig.  4).  The  results  shown 
in  Fig.  4  are  typical  of  temporal  fluence  rate  measurements  among 
patients  in  the  study.  It  shows  very  little  variation  of  optical 
properties  for  most  sites,  except  for  LUQ  where  the  fluence  rate 
seems  to  decrease.  These  measurements  were  performed  in  such 
a  way  that  the  same  site  was  used  before  and  after  PDT.  Our 
general  conclusion  is  consistent  with  published  results  at  630  nm 
for  Photofrin-mediated  PDT  (1 1)  and  at  732  nm  for  MLu-mediated 
PDT  (12)  in  the  canine  prostate  model,  where  little  time 
dependence  of  light  fluence  rate  was  observed.  However,  because 
of  large  heterogeneities,  variation  in  optical  property  measurements 
before  and  after  PDT  was  found  (see  Fig.  5). 

There  are  several  potential  reasons  why  no  significant  light 
fluence  rate  was  measured  in  the  urethra  (Fig.  4).  First  of  all,  it 
could  be  due  to  measurement  uncertainty.  The  light  treatment  was 
performed  one  quadrant  at  a  time  rather  than  all  four  quadrants 
together.  Because  the  range  of  the  light  penetration  in  human 
prostate  is  relatively  short  (0.4  cm),  only  those  CDF  (usually  in  the 
RUQ  and  the  LUQ)  near  the  urethra  will  contribute  to  the  urethra. 
Because  the  position  of  the  isotropic  detector  in  urethra  was 


Figure  5.  Variation  of  (a)  pa,  (b)  p',  (c)  g<,n-  and  (d)  8  in  patients  before 
(open  bars)  and  after  (solid  bars)  PDT.  The  height  of  each  bar  is  the  mean 
value  measured  in  each  patient.  The  error  bars  specify  the  standard 
deviation  among  the  four  quadrants  of  a  prostate.  The  solid  horizontal  line 
is  the  average  of  all  measurements,  and  the  dashed  lines  are  the  standard 
deviation  of  the  average. 


optimized  only  once  during  the  treatment  of  one  quadrant,  it  is 
possible  that  we  had  not  obtained  the  peak  light  fluence  rate  for  the 
entire  treatment.  (The  position  of  the  urethral  isotropic  detector 
could  not  be  optimized  multiple  times  because  one  can  only  move 
the  position  of  the  urethra  isotropic  detector  in  one  direction.) 
Second,  the  light  fluence  could  actually  be  lower  in  the  urethra. 
Besides  the  fact  that  we  have  avoided  placing  light  sources  through 
the  urethra,  we  have  noticed  increased  absorption  in  the  urethra  due 
to  high  concentrations  of  MLu  in  urine,  which  will  reduce  the  light 
fluence  rate  reaching  the  isotropic  detector  and  thus,  the  urethra 
itself.  No  serious  urethral  toxicity  has  been  observed  in  all  the 
patients  treated  so  far. 
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Figure  5.  Continued. 


The  accuracy  of  the  method  used  to  determine  optical  properties 
from  measurements  around  a  point  source  was  dominated  by  the 
uncertainty  of  the  distance  between  the  light  source  and  detector 
catheters  ( h ,  see  Fig.  lb).  This  is  because  the  inserted  catheters  are 
not  exactly  parallel  to  each  other  and  an  error  of  1-2  mm  is 
possible,  despite  the  use  of  ultrasound  guidance  to  position  the 
catheters  correctly.  A  detailed  error  analysis  is  included  in  the 
Appendix.  To  reduce  the  positioning  error,  the  later  measurements 
(starting  with  Patient  9)  were  made  using  a  motorized  probe,  with 
a  positioning  accuracy  of  at  worst  0.2  mm.  The  use  of  the 
motorized  probe  also  increased  the  number  of  measurement  points 
from  typically  four  points  per  measurement  to  approximately  800 
points  per  measurement.  We  then  allowed  the  distance  h  to  be  one 
of  the  free  parameters  of  the  fit  (ga>  fi's  and  h),  and  h  was  allowed  to 
be  different  from  the  known  separation  based  on  the  template  by  1— 
2  mm.  Using  this  method  our  phantom  studies  confirmed  that  one 
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Table  2.  Summary  of  optical  properties  before  and  after  PDT  in  human  prostates.  The  values  in  the  parentheses  are  the  standard  deviation  of  the  mean 
values  measured  from  different  locations  of  the  same  prostate.  No  standard  deviation  is  listed  if  only  one  data  point  is  available 


Before  PDT  After  PDT 


Patient  number 

Pa  (cm  ‘) 

fls  (cm  ') 

5  (cm) 

<j )/S  (cm  2) 

Pa  (cm  ‘) 

r4(cm  ') 

8  (cm) 

<| >/S  (cm  2) 

2 

0.09 

29.8 

0.34 

3.34 

0.09 

43.7 

0.29 

3.78 

3 

0.15 

22.0 

0.31 

2.15 

0.07 

33.4 

0.37 

4.07 

4 

0.43  (0.28) 

7.69  (4.76) 

0.41  (0.14) 

0.97  (0.81) 

0.51 

1.67 

0.63 

0.36 

5 

0.21 

11.8 

0.37 

1.44 

0.13 

7.18 

0.60 

1.48 

6 

0.27  (0.27) 

10.5  (11.2) 

0.50  (0.05) 

1.74  (1.78) 

0.19  (0.20) 

18.9  (18.4) 

0.45  (0.06) 

3.18  (3.32) 

7 

— 

— 

— 

— 

0.30  (0.08) 

23.7  (13.9) 

0.24  (0.11) 

1.09  (0.40) 

9 

0.53  (0.36) 

6.61  (4.51) 

0.41  (0.09) 

0.77  (0.42) 

0.64  (0.25) 

7.00  (5.59) 

0.33  (0.10) 

0.54  (0.23) 

10 

0.63  (0.32) 

4.62  (2.87) 

0.42  (0.10) 

0.56  (0.29) 

0.19  (0.05) 

9.27  (4.47) 

0.54  (0.31) 

1.33  (0.64) 

11 

0.67  (0.17) 

6.39  (3.18) 

0.32(0.10) 

0.51  (0.18) 

0.83  (0.45) 

5.45  (3.89) 

0.38  (0.16) 

0.47  (0.30) 

12 

0.71  (0.43) 

8.99  (6.51) 

0.32  (0.12) 

0.61  (0.35) 

0.30  (0.06) 

20.2  (4.8) 

0.28  (0.08) 

0.98  (0.05) 

13 

0.27  (0.14) 

18.5  (11.6) 

0.30  (0.07) 

1.46(0.72) 

0.26  (0.09) 

17.0  (8.8) 

0.31  (0.07) 

1.42  (0.56) 

14 

0.72  (0.11) 

3.37  (1.37) 

0.39  (0.11) 

0.40  (0.01) 

— 

— 

— 

— 

can  determine  4>(0.5)  and  pa  to  within  ±20%  for  most  cases.  The 
error  for  Peff,  is  much  smaller  (±5%)  and  is  insensitive  to  the 
positioning  error  in  h  (see  Appendix). 

Some  of  the  measurements  cannot  be  used  to  extrapolate  optical 
properties  because  of  bleeding.  This  happens  when  blood  surround¬ 
ing  one  of  the  catheters  (either  source  catheter  or  the  detector 
catheter)  blocks  the  light  from  reaching  the  detector.  We  used  a 
standard  of  <J>/S  >0.1  cm-2  to  determine  if  there  was  sufficient 
signal  to  analyze  the  data.  The  average  success  rate  for  the  first 
seven  patients  was  very  poor  (approximately  35%).  This  was 
increased  to  64%  for  later  patients  using  a  motorized  probe,  pri¬ 
marily  because  more  measurements  can  be  made  in  a  shorter  period 
of  time. 

The  overall  variation  of  optical  properties  (pa,  p',  peff  and  5=1/ 
pCff)  before  and  after  PDT  for  different  patients  is  shown  in  Fig.  5 
to  demonstrate  interpatient  heterogeneity.  The  heterogeneity  in 
optical  properties  was  the  major  reason  that  the  light  fluence  rate 
a  fixed  distance  away  from  the  CDF  varied  from  patient  to  patient. 
However,  the  heterogeneity  of  optical  properties  in  human  prostate 
is  somewhat  smaller  than  that  observed  in  canine  prostate  at  732 
nm  (12).  Overall  pa  varied  between  0.07  and  1.62  cm-1,  and  p' 
varied  between  1.1  and  44  cm-1.  pCff  varied  between  0.91  and  6.7 
cm-1,  corresponding  to  an  optical  penetration  depth  (5  =  1/fieff)  of 
0.2-1. 1  cm.  The  mean  values  of  peff  and  5  were  2.9  ±  0.7  cm-1 
and  0.4  ±  0.1  cm,  respectively.  This  penetration  depth  is 
substantially  larger  than  that  of  0.1-0.25  cm  predicted  for  630 
nm  (11)  but  is  smaller  than  0.5-3  cm  observed  in  normal  canine 
prostate  at  732  nm  (12).  The  most  probably  explanation  is  that 
canine  prostate  has  different  grandular/structure  content  than 
that  of  human  prostate.  Whereas  the  mean  reduced  scattering 
coefficient  in  canine  prostate  was  3.6  ±  4.8  cm1  (12),  it  was  14  ± 
1 1  cm-1  in  human  prostrate  at  the  same  wavelength  (732  nm).  The 
increased  reduced  scattering  coefficient  resulted  in  increased  peff, 
or  a  reduction  of  optical  penetration  depth,  assuming  that  the  pa 
remains  the  same.  In  addition,  the  increase  in  pas  in  human  pros¬ 
tates  compared  with  that  in  canine  prostates  contributed  to  the 
reduction  of  the  penetration  depth  in  humans  compared  with  that  in 
canines.  Human  prostates  measured  in  this  study  had  cancer 
present  and  had  been  subjected  to  previous  radiation  therapy.  The 
canine  prostate  is  a  normal  tissue  model.  Hypoxia  has  been  des¬ 
cribed  in  patients  with  prostate  cancer  (26),  whereas  one  would 
expect  normal  oxygen  levels  in  normal  canine  prostate  tissue.  As 


a  result,  there  may  be  more  deoxyhemoglobin  in  these  prostates 
than  in  normal  prostates,  resulting  in  higher  pas  and  shorter  pene¬ 
tration  depths  at  the  wavelength  range  630-800  nm. 

The  interprostatic  difference  in  fluence  rate  per  unit  source 
strength  at  0.5  cm  from  the  point  source  was  shown  in  Fig.  6.  We 
chose  0.5  cm  because  this  is  the  middle  point  for  1  cm  spaced 
interstitial  loading.  The  mean  value  was  1,5  ±  1.1  cm-2  for  all 
prostates  studied  (sample  size  121).  We  observed  interpatient 
fluence  rate  heterogeneity  as  a  result  of  heterogeneity  of  optical 
properties. 

The  tissue  concentration  of  MLu  and  the  pa  in  the  human 
prostate  at  732  nm  as  shown  in  Fig.  7a  demonstrates  a  linear 
relationship  with  measured  MLu  concentration  (correlation  co¬ 
efficient  R2  =  0.63).  This  relationship  is  described  by  the  equation 
pa  =  0.066c  +  0.23.  The  relatively  large  spread  of  the  data  can  be 
explained  by  the  large  heterogeneity  of  MLu  distribution  in 
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Figure  6.  Measured  mean  fluence  rate  per  power.  <)>/S,  at  0.5  cm  from  the 
point  source  at  732  nm  in  different  patients  before  (open  bars)  and  after 
(solid  bars)  PDT.  The  error  bars  reflect  intraprostate  variation  of  c|>/S.  The 
solid  horizontal  line  is  the  average  of  all  measurements,  and  the  dashed 
lines  are  the  standard  deviation  of  average:  1.5  ±  1.1  cm-2. 
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Inject  MLu  dose  and  drug-light  interval 


Figure  7.  (a)  In  vivo  absorption  coefficient  (pa)  at  732  nm  as  a  function  of 
measured  ex  vivo  MLu  concentration  (in  ng/mg)  from  the  same  human 
prostate.  The  symbols  are  measured  points  in  vivo :  O — before  PDT,  * — 
after  PDT.  The  solid  line  is  a  linear  fit  to  the  measured  in  vivo  data.  The 
symbol  “+”  represents  the  measured  pa  of  MLu  vs  MLu  concentration  in 
a  pure  liposyn  phantom.  The  dashed  line  is  a  linear  fit  to  the  phantom  data 
with  ps'  =  4  cm-1,  (b)  Extrapolated  MLu  concentration  from  ex  vivo 
fluorescence  (bars)  and  from  pa  (symbols)  as  a  function  of  injected  MLu 
concentration:  open  bars  (O) — before  PDT,  solid  bar  (□) — after  PDT.  Drug 
concentration  calculated  from  pa  uses  the  linear  relationship  determined  in 
(a):  c  (ng/mg)  =  (pa  (cm-1)  -  0.227)/0.0658.  The  thin  and  thick  error  bars 
correspond  to  the  standard  deviation  of  the  ex  vivo  biopsy  and  in  vivo 
absorption  measurements,  respectively.  Notice  that  (b)  used  all  measured 
pa  data  (121  points),  whereas  (a)  used  a  subset  of  pa  data  (24  points)  that  had 
corresponding  ex  vivo  biopsies.  (Insufficient  data  existed  to  extrapolate 
tissue  concentration  from  in  vivo  absoiption  measurement  for  1  mg/kg  MLu 
injection  at  6  h  drug-treatment  time  interval  before  PDT  due  to  bleeding.) 


different  locations  of  a  human  prostate  (27).  Because  the  location 
of  the  in  situ  optical  measurements  and  the  tissue  biopsies  were  not 
exactly  the  same,  some  discrepancy  between  the  drug  levels 
measured  by  the  two  techniques  is  expected.  Absorption  by  many 
naturally  occurring  chromophores,  such  as  melanins,  cytochromes, 


fiavoproteins,  and  in  particular,  hemoglobin  and  deoxyhemoglobin, 
is  small  (around  732  nm)  (10-13);  as  such,  it  is  not  surprising  that 
the  pa  at  732  nm  varies  approximately  linearly  with  the  MLu 
concentration.  Given  the  heterogeneity  in  the  patient  data,  we  used 
a  liposyn  phantom  to  examine  the  relationship  between  MLu 
concentration  and  the  Pa  and  to  confirm  the  slope  of  our  clinical 
observations.  Data  based  on  this  phantom  are  shown  in  Fig.  7a  by 
the  “+”  symbols.  The  absence  of  other  absorbers  in  the  phantom 
decreases  pa.  The  slope  of  the  best  fit  line  (0.062  cm-1/(ng/mg))  is 
consistent  with  that  observed  in  vivo.  The  linear  relationship 
between  MLu  drug  concentration  and  the  pa  at  732  nm  is  consistent 
with  that  observed  in  canine  prostate  at  732  nm,  but  a  larger 
intersect  (pa  =  0.23  cm-1)  is  found  in  humans  than  in  dogs  (pa  = 
0.08  cm-1)  (12).  This  may  be  caused  by  differences  in  the 
oxygenation  of  cancerous  prostate  tissue  (humans)  and  normal 
prostate  tissue  (dogs)  because  the  attenuation  coefficient  of 
oxyhemoglobin  at  732  nm  is  substantially  less  than  that  of 
deoxyhemoglobin.  If  the  tissue  is  less  oxygenated,  the  attenuation 
coefficient  will  increase.  The  conclusion  that  more  deoxyhemo¬ 
globin  existed  in  human  prostate  is  confirmed  by  in  vivo  diffuse 
absorption  spectroscopy  in  the  same  prostate  patient  population 
(28).  For  the  two  patients  studied,  the  contribution  to  the  pa  from 
oxyhemoglobin  and  deoxyhemoglobin  varied  between  0.1  and 
0.2  cm-1  at  732  nm. 

The  effect  of  injected  drug  dose  and  drug-light  interval  on  tissue 
concentration  of  MLu  is  shown  in  Fig.  7b.  It  shows  significant 
intra-  and  interpatient  variation  of  drug  concentration.  The 
extrapolated  mean  tissue  concentration  from  in  vivo  absorption 
(symbols)  generally  agrees  with  that  from  ex  vivo  biopsy.  Figure  7b 
was  drawn  using  all  121  absorption  measurements  available  as 
opposed  to  Fig.  7a,  which  used  a  subset  of  data  (24  points)  that  had 
corresponding  ex  vivo  biopsies.  The  large  discrepancy  between 
drug  concentration  measured  in  vivo  and  ex  vivo  in  the  two  patients 
receiving  0.5  mg/kg  with  a  24  h  drug-treatment  interval  is  likely 
a  consequence  of  measurement  uncertainties  for  in  vivo  absorption 
in  this  early  group  of  patients.  The  detector  positions  for  these 
patients  were  determined  manually  using  a  ruler,  which  introduced 
larger  errors  in  the  final  determination  of  pa.  In  addition,  because 
p„  for  the  earlier  patients  were  small  (less  MLu),  the  relative 
uncertainty  in  pa  increased.  The  large  discrepancy  in  drug 
concentration  measured  in  vivo  and  ex  vivo  in  the  group  receiving 
1  mg/kg  with  a  6  h  drug-treatment  interval  can  be  explained  by  the 
small  number  of  samples,  i.e.  patients,  in  this  group  (»  =  1). 
However,  some  preliminary  conclusions  can  be  made  from  these 
data.  It  appears  that  a  reduction  in  drug  concentration  is  observed 
after  PDT,  which  may  be  consistent  with  photobleaching.  The 
patients  who  received  higher  injected  MLu  doses  generally  had 
higher  local  drug  concentrations  in  the  prostate.  If  one  estimates  the 
MLu  concentration  for  2  mg/kg  at  24  h  in  vivo  as  twice  the  MLu 
concentration  for  1  mg/kg  at  24  h,  then  the  drug  concentration 
seems  to  decrease  systematically  with  incubation  time,  i.e. 
drug-treatment  interval,  by  examining  the  data  for  2  mg/kg  at  24 
h,  2  mg/kg  at  6  h  and  2  mg/kg  at  3  h  in  Fig.  7b.  The  heterogeneity 
in  MLu  concentration  is  sufficiently  large  to  potentially  require  in 
situ  monitoring  of  MLu  distribution  for  future  PDT  treatments  of 
the  human  prostate. 

When  one  examines  the  variation  of  drug  concentration  in  a 
single  patient,  it  is  clear  that  the  drug  concentration  can  easily 
change  by  three  times  when  measured  from  different  locations  of 
the  same  prostate  (data  presented  in  reference  [27]).  This  strong 
variation  suggests  that  the  mean  intraprostatic  variation  is  usually 
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larger  than  the  mean  interprostatic  variation,  as  shown  in  Fig.  7b, 
where  the  difference  between  bars  (e.g.  before  PDT)  for  different 
groups  is  smaller  than  the  range  of  the  error  bar,  representing  the 
standard  deviation  of  intraprostatic  variation  in  the  same  group; 
most  of  this  variation  is  caused  by  intraprostatic  variation. 

CONCLUSIONS 

We  have  measured  the  optical  properties  and  MLu  concentration  of 
human  prostates  for  732  nm  light.  Substantial  inter-  and  intra¬ 
patient  heterogeneity  was  observed.  The  mean  light  penetration 
depth  in  the  human  prostates  is  about  0.4  cm,  at  least  two  times 
smaller  than  those  found  in  the  normal  canine  prostates.  The  pa  in 
vivo  at  732  nm  has  been  found  to  be  proportional  to  the  tissue 
concentration  of  MLu,  determined  from  ex  vivo  biopsies.  A  real¬ 
time  dosimetry  measurement  and  feedback  system  for  monitoring 
light  fiuences  during  treatment  should  be  considered  for  such  future 
PDT  studies  to  account  for  this  heterogeneity. 

Acknowledgements — This  work  is  supported  by  grants  from  National  In¬ 
stitute  of  Health  (NIH),  R21  CA88064-01  and  P01  CA87971-01,  and  De¬ 
partment  of  Defense  (DOD),  DAMD1703-1-0132.  We  thank  the  referees 
for  their  helpful  critiques  and  constructive  suggestions  for  the  manuscript. 

APPENDIX 

Error  analysis  of  the  optimization  technique  to  extrapolate 
optical  properties 

There  are  several  uncertainties  in  determining  the  optical  properties 
besides  uncertainty  in  measuring  the  light  fluence  rate,  which 
accounts  for  approximately  ±5%.  One  is  the  uncertainty  in  the 
position  0)  of  the  detector,  which  can  be  1-2  mm.  Another 
uncertainty  is  the  uncertainty  in  the  distance  ( h )  between  the  source 
catheter  and  the  isotropic  detector  catheter,  which  could  be  different 
from  the  distances  specified  on  the  template  by  1-2  mm  as  well.  Both 
errors  introduce  an  uncertainty  in  distance  (Fig.  lb),  r  =  yx2  +  h2, 
which  can  be  1-2  mm.  From  the  phantom  measurements,  one  can 
estimate  the  relative  uncertainty  introduced  by  the  positioning  errors 
by  the  ratio  of  the  slope  of  the  fluence  rate  to  fluence  rate  d$/dr/§dr. 
Because  the  fluence  rate  at  distance  r  from  a  point  source  is 
proportional  to  e~l''"r/r,  this  value  decreases  with  increasing  distance 
r  and  follows  (p eff+l/r)dr.  Using  this  relationship  one  gets  a  relative 
uncertainty  in  fluence  rate  of  50%  for  r = 5  mm  and  39%  for  r = 7  mm 
assuming  that  dr  —  2  mm  and  peff  <  0.5  cm-1.  To  separate  the 
two  sources  of  errors  that  contribute  to  an  error  in  r,  one  writes  dr  = 
x/rdx+h/rdh.  The  first  error,  x,  is  a  random  error  and  approaches  dr = 
0  for  x  =  0  for  any  dx  because  dr  =  xdx/r.  Thus,  the  overall  error 
caused  by  uncertainty  of  x  for  determining  pa  and  pj  should  be  small. 
This  is  also  because  it  is  unlikely  that  all  points  are  off  in  the  same 
direction  thereby  introducing  a  systematic  error.  We  have  minimized 
this  error  by  using  a  computerized  step  motor  so  that  the  error  in  x  was 
less  than  0.2  mm.  In  this  case,  the  error  in  <j)  should  be  less  than  2% 
even  for  x  =  r.  The  second  source  of  error  can  be  large  but  it  is 
a  systematic  error. 

Table  A  shows  a  numerical  simulation  of  the  errors  in 
pa  and  p'  caused  by  the  errors  in  h  for  two  conditions  of 
pa  and  p'.  The  standard  deviation  of  the  fit  is  calculated  by 

J2, Li  ((measurement,-  —  fit;)/fit,)2/(n  —  1).  One  can  use  the 
minimum  error  to  determine  the  correct  h  value  to  within  0.05 
cm.  In  Table  A,  this  corresponds  to  the  fit  for  h  =  0.55  cm.  The 


Table  A.  Error  analysis  of  deviations  in  pa,  p(,  pefr  and  ij>/S(0.5  cm) 
caused  by  the  uncertainty  of  h  (nominal  h  =  0.6  cm) 


h  (cm) 

pa  (cm  ’) 

P's  (cm  ' 

)  Peff  (cm  *) 

<j>/S  (0.5  cm) 

Standard 
deviation 
of  fit  (%) 

o 

II 

=2 

cm-1,  p(  = 

3.6  cm-1 

,  Peff  =  2.3  cm" 

1 

0.40 

1.07 

1.17 

1.93 

0.28 

27.6 

0.45 

0.81 

1.76 

2.07 

0.35 

17.2 

0.50 

0.65 

2.48 

2.20 

0.43 

11.9 

0.55 

0.53 

3.37 

2.32 

0.53 

10.9 

0.60 

0.44 

4.57 

2.45 

0.66 

12.7 

0.65 

0.37 

5.99 

2.57 

0.81 

15.2 

0.70 

0.30 

7.98 

2.69 

1.00 

17.8 

0.75 

0.25 

10.65 

2.81 

1.26 

20.3 

0.80 

0.20 

14.30 

2.93 

1.58 

22.4 

pa  =  0.4  cm  \  p(  = 

7.2  cm-1 

l.  Peff  =  2.9  cm" 

1 

0.40 

0.70 

2.75 

2.41 

0.43 

15.2 

0.45 

0.56 

3.91 

2.57 

0.54 

8.7 

0.50 

0.46 

5.41 

2.73 

0.68 

5.1 

0.55 

0.38 

7.34 

2.90 

0.83 

4.8 

0.60 

0.30 

10.33 

3.04 

1.09 

7.0 

0.65 

0.24 

14.31 

3.19 

1.39 

9.6 

0.70 

0.19 

19.96 

3.34 

1.80 

12.1 

0.75 

0.14 

28.09 

3.49 

2.35 

14.3 

0.80 

0.11 

39.95 

3.64 

3.10 

16.3 

physical  separation  is  h  =  0.6  cm.  We  do  not  understand  why 
the  optimal  h,  which  gives  the  best  pa  and  p',  is  about  0.05  cm 
shorter  than  the  physical  separation  of  h  =  0.6  cm.  But  we  have 
demonstrated  that  one  can  extract  h  to  within  0.05  cm.  Using 
this  additional  constraint  and  allowing  h  to  deviate  from  the 
nominal  value  will  reduce  the  final  uncertainty  to  less  than  ±6% 
for  pa  and  less  than  ±11%  for  p'.  Assuming  that  the  h  can  still 
deviate  from  the  true  h  (0.55  cm)  by  ±0.05  cm,  we  estimated 
that  one  can  determine  pa  to  an  uncertainty  of  less  than  ±20% 
and  p'  to  an  uncertainty  of  less  than  ±40%,  provided  that  the 
optimal  h  is  within  0.5  mm  of  the  true  value  of  h.  Clearly,  if  one 
allows  up  to  2  mm  uncertainty  in  h,  then  one  can  introduce  an 
unacceptable  error  for  pa  and  p'  (see  Table  A)  with  up  to 
±215%  and  ±450%  errors,  respectively.  The  error  in  pcff  is 
much  smaller  (±30%)  over  the  same  range  of  h.  To  reduce  this 
error  further,  we  have  used  the  value  of  the  standard  deviation  of 
the  fit  as  a  guide  to  determine  the  optimal  separation  h. 
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Abstract 

It  is  desirable  to  quantify  the  distribution  of  the  light  fluence  rate,  the  optical  properties,  the  drug  concentration,  and  the  tissue 
oxygenation  for  photodynamic  therapy  (PDT)  of  prostate  cancer.  We  have  developed  an  integrated  system  to  determine  these  quan¬ 
tities  before  and  after  PDT  treatment  using  motorized  probes.  The  optical  properties  (absorption  (pa),  transport  scattering  (/;'),  and 
effective  attenuation  (pefr)  coefficients)  of  cancerous  human  prostate  were  measured  in-vivo  using  interstitial  isotropic  detectors. 
Measurements  were  made  at  732  nm  before  and  after  motexafin  lutetium  (MLu)  mediated  PDT  at  different  locations  along  each 
catheter.  The  light  fluence  rate  distribution  was  also  measured  along  the  catheters  during  PDT.  Diffuse  absorption  spectroscopy 
measurement  using  a  white  light  source  allows  extrapolation  of  the  distribution  of  oxygen  saturation  (St02),  total  blood  volume 
([Hb]t),  and  MLu  concentration.  The  distribution  of  drug  concentration  was  also  studied  using  fluorescence  from  a  single  optical 
fiber,  and  was  found  to  be  in  good  agreement  with  the  values  determined  by  absorption  spectroscopy.  This  study  shows  significant 
inter-  and  intra-prostatic  variations  in  the  tissue  optical  properties  and  MLu  drug  distribution,  suggesting  that  a  real-time  dosimetry 
measurement  and  feedback  system  for  monitoring  these  values  during  treatment  should  be  considered  in  future  PDT  studies. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Prostate  adenocarcinoma  is  the  most  common  malig¬ 
nancy  in  men.  In  2003,  it  was  estimated  that  220,900  cases 
of  prostate  adenocarcinoma  were  diagnosed  in  the  United 
States  [1],  Although  the  availability  of  serum  prostate- 
specific  antigen  (PSA)  measurement  as  a  screening  tool 
has  resulted  in  earlier  detection  of  the  disease  [2],  prostate 
cancer  still  accounted  for  28,900  deaths  in  2003  [1]. 

Photodynamic  therapy  (PDT)  is  an  emerging  treat¬ 
ment  modality  based  on  the  interaction  of  light,  a  pho- 
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tosensitizing  drug,  and  oxygen  [3].  PDT  has  been  a 
proposed  treatment  for  a  variety  of  malignancies  and 
premalignant  conditions.  PDT  has  been  approved  by 
the  US  Food  and  Drug  Administration  for  the  treat¬ 
ment  of  microinvasive  lung  cancer,  obstructing  lung 
cancer,  and  obstructing  esophageal  cancer.  The  prostate 
gland  is  an  organ  that  appears  to  be  a  good  target  for 
interstitial  PDT.  Tumors  of  the  prostate  are  often  con¬ 
fined  to  the  prostate  itself  and  brachytherapy  techniques 
used  for  the  placement  of  radioactive  seed  implants  can 
be  adapted  for  the  placement  of  interstitial  optical  fibers 
[4]-  Several  preclinical  studies  have  evaluated  the  feasi¬ 
bility  of  delivering  PDT  to  the  prostate  via  this  intersti¬ 
tial  approach  [5-10],  A  trial  of  interstitial  prostate  PDT 
in  humans  has  been  reported  [11],  Nathan  et  al.  treated 
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14  men  with  locally  recurrent  prostate  cancer  using 
meso-tetrahydroxyphenyl  chlorin  (mTHPC)-mediated 
interstitial  PDT.  The  light  treatment  was  directed 
against  regions  from  which  biopsies  showed  cancer  or 
which  were  suspicious  on  imaging  studies. 

Motexafin  lutetium  (MLu)  is  a  water-soluble  second 
generation  synthetic  photoactive  drug  that  has  a  Q-band 
absorption  peak  at  732  nm  [12,13],  Based  on  a  feasibility 
and  toxicity  study  in  a  canine  prostate  model  [9],  we 
have  started  a  phase  I  study  of  MLu-mediated  PDT 
for  prostate  cancer  [14].  In  the  canine  study,  comprehen¬ 
sive  treatment  of  the  prostate  gland  was  achieved  with 
MLu-mediated  PDT  using  an  interstitial  approach. 
The  development  of  this  light  delivery  technique  has 
necessitated  an  improved  understanding  of  PDT  dosi¬ 
metry,  critical  for  determining  the  efficacy  of  the  PDT 
treatment. 

Explicit  PDT  dosimetry  includes  quantifying  the  light 
and  tissue  optical  properties,  the  drug  concentration, 
and  the  tissue  oxygenation.  The  light  fluence  (expressed 
in  J/cm2)  is  a  measure  of  light  energy  imparted  to  tissue. 
The  total  fluence  in  tissues  is  a  function  not  only  of  the 
incident  light  delivered  by  the  laser  but  also  of  scattered 
light.  Often  clinical  PDT  treatments  are  prescribed  in 
terms  of  the  incident  light  delivered  from  the  laser  rather 
than  the  total  fluence  of  light  the  tissues  receive  which  is 
a  combination  of  scattered  and  incident  light.  Dosimetry 
systems  using  isotropic  light  detectors  have  been  devel¬ 
oped  to  measure  both  incident  and  scattered  light 
[15,16].  These  systems  allow  us  to  measure  and  therefore 
prescribe  a  consistent  total  fluence  to  the  tissues. 

Several  investigators  have  attempted  to  characterize 
the  optical  properties  of  prostate  tissue  in  animals  [17- 
19]  and  in  humans  [20-22]  to  more  reliably  predict  the 
in  vivo  light  distribution.  Using  diffusion  theory  for  a 
point  source,  the  absorption  and  transport  scattering 
coefficients  of  a  particular  tissue  can  be  determined, 
yielding  the  effective  attenuation  coefficient,  which  pro¬ 
vides  a  measure  of  light  penetration  in  that  tissue  [18]. 
This  measurement  is  a  critical  factor  in  planning  intersti¬ 
tial  light  source  placement.  To  include  the  drug  concen¬ 
tration  in  the  evaluation  of  PDT  efficacy,  in  situ 
measurements  of  photosensitizer  fluorescence  emission 
are  made  in  the  prostate  using  a  single  optical  fiber,  orig¬ 
inally  developed  for  surface  application  by  Diamond 
et  al.  [23],  We  have  modified  their  technique  by  replacing 
the  flat  cut  fiber  with  a  side  fire  fiber  to  introduce  light 
interstitially  to  the  target  tissue  more  efficiently.  It  has 
also  been  shown  that  one  can  determine  MLu  concen¬ 
tration  using  diffuse  absorption  spectra  [24].  The  MLu 
concentration  in  prostate  tissue  can  be  determined  using 
a  slightly  modified  technique  via  an  interstitial  approach 
[25].  The  MLu  tissue  concentrations  determined  from 
absorption  spectra  can  be  compared  to  those  obtained 
using  fluorescence  spectra.  Using  the  same  diffuse 
absorption  spectra,  it  is  also  possible  to  determine  the 


concentrations  of  deoxyhemoglobin  ([Hb])  and  oxyhe¬ 
moglobin  ([Hb02])  [24-27],  To  determine  the  tissue  oxy¬ 
genation,  the  oxygen  saturation  ratio  St02  =  [Hb02]/ 
([Hb]  +  [Hb02])  and  total  hemoglobin  concentration 
[Hb]t  =  [Hb]  +  [Hb02]  can  be  calculated. 

In  this  study,  we  review  our  preliminary  observations 
measuring  the  interstitial  distribution  of  light  fluence 
rate,  optical  properties,  tissue  MLu  concentration,  and 
tissue  oxygenation  in  two  patients  with  prostate  cancer. 
Measurements  were  made  before  and  after  PDT  treat¬ 
ment  in  which  the  entire  prostate  was  treated.  Intrapro¬ 
static  as  well  as  interprostatic  differences  were  evaluated. 

2.  Materials  and  methods 

2.1.  Patient  selection,  surgical  and  PDT  procedure 

A  Phase  I  clinical  trial  of  PDT  with  MLu  in  patients 
with  locally  recurrent  prostate  carcinoma  was  initiated 
at  the  University  of  Pennsylvania  in  1999.  The  protocol 
was  approved  by  the  Institutional  Review  board  of  the 
University  of  Pennsylvania,  the  Clinical  Trials  and  Sci¬ 
entific  Monitoring  Committee  (CTSRMC)  of  the  Uni¬ 
versity  of  Pennsylvania  Cancer  Center,  and  the  Cancer 
Therapy  Evaluation  Program  (CTEP)  of  the  National 
Cancer  Institute.  Each  patient  who  signed  the  informed 
consent  document  underwent  an  evaluation,  which  in¬ 
cluded  an  MRI  of  the  prostate,  bone  scan,  laboratory 
studies  including  PSA,  and  a  urological  evaluation. 
Approximately  two  weeks  prior  to  the  scheduled  treat¬ 
ment  a  transrectal  ultrasound  (TRUS)  was  performed 
for  treatment  planning.  A  urologist  drew  the  target  vol¬ 
ume  (the  prostate)  on  each  slice  of  the  ultrasound 
images.  These  images  were  spaced  0.5  cm  apart  and  were 
scanned  with  the  same  ultrasound  unit  used  for  treat¬ 
ment.  A  built-in  template  with  a  0.5-cm  grid  projected 
the  locations  of  possible  light  sources  relative  to  the 
prostate.  A  treatment  plan  was  then  prepared  to  deter¬ 
mine  the  location  and  length  of  light  sources.  Cylindrical 
diffusing  fibers  (CDF’s)  with  active  lengths  1,  2,  3, 4,  and 
5  cm  were  used.  The  sources  were  spaced  one  centimeter 
apart  and  the  light  power  per  unit  length  was  the  same 
for  all  CDF’s.  The  length  of  a  light  source  at  a  particular 
position  was  selected  to  cover  the  full  length  of  the  pros¬ 
tate.  In  the  final  plan,  the  prostate  was  divided  into  four 
quadrants.  Four  isotropic  detectors  were  used,  each 
placed  in  the  center  of  one  quadrant.  A  fifth  isotropic 
detector  was  placed  in  a  urethral  catheter  to  monitor 
the  light  fluence  in  the  urethra. 

The  patients  were  anesthetized  in  the  operating  room 
with  general  anesthesia  to  minimize  patient  movement 
during  the  procedure.  Transrectal  ultrasound-guided 
biopsies  for  MLu  measurements  were  obtained  prior 
to  light  delivery.  The  same  ultrasound  unit  used  to  per¬ 
form  the  pretreatment  TRUS  was  used  to  guide  needle 
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placement  in  the  operating  room.  A  template  was  at¬ 
tached  to  the  ultrasound  unit  and  was  matched  to  the 
0.5-cm  grid  used  for  treatment  planning.  Four  detector 
catheters  (one  for  each  quadrant)  were  inserted  into 
the  prostate.  These  detectors  were  kept  in  place  during 
the  entire  light  delivery  period.  The  pre-planned  treat¬ 
ment  catheters  for  light  sources  were  then  inserted  0.5 
or  0.7  cm  away  from  the  detector  catheters.  These 
source  catheters  were  used  for  light  delivery  as  well  as 
optical  properties  and  absorption  spectra  measure¬ 
ments.  The  optical  properties  and  absorption  spectra 
of  the  prostate  in  each  quadrant  were  then  measured  be¬ 
fore  light  delivery  using  these  existing  detectors.  A  point 
light  source  was  inserted  into  one  source  catheter.  Fluo¬ 
rescence  measurements  were  made  through  the  detector 
catheter  in  each  quadrant  using  a  single  optical  fiber  act¬ 
ing  as  both  a  light  source  and  a  detector.  The  fiber  intro¬ 
duced  460  nm  light  and  collected  fluorescence  light 
above  700  nm  at  right  angles  from  the  optical  axis  of 
the  beveled  fiber  tip.  The  distribution  measurements  of 
optical  properties,  absorption  spectra  and  fluorescence 
were  completed  in  approximately  10  min.  Light  delivery 
was  then  performed  one  quadrant  at  a  time  by  inserting 
the  CDFs  into  the  source  catheters.  The  treatment  time 
for  each  quadrant  was  dependent  on  the  detector  read¬ 
ing  in  that  quadrant.  The  light  fluence  rate  distributions 
were  measured  during  PDT  along  the  catheters.  Cumu¬ 
lative  fluences  of  50  and  100  J/cm2  were  delivered  to  pa¬ 
tients  12  and  13,  respectively.  After  light  delivery,  the 
optical  properties  and  absorption  and  fluorescence  of 
the  prostate  in  all  four  quadrants  were  measured  again. 
The  light  sources  and  detectors  were  then  removed  and 
post-treatment  biopsies  were  performed. 

2.2.  Computer  controlled  step  motors 

Two  step  motors  were  used  to  control  the  movement 
of  the  light  source  and  the  isotropic  detector  during 
optical  property  measurements  (see  Fig.  1  (a)).  Each  step 
motor  and  translational  stage  had  a  maximum  range  of 
20  cm  and  a  maximum  speed  of  12.5  mm/s.  The  step  mo¬ 
tor  produced  400  pulses  per  rotation,  which  translated 
to  a  resolution  of  0.0025  mm.  Control  software  was 
developed  to  integrate  the  movement  of  the  step  motor 
with  data  acquisition  of  isotropic  detectors  as  well  as 
the  spectrometer  for  fluorescence  and  absorption  spectra 
measurements.  For  the  integrated  system,  the  position¬ 
ing  accuracy  was  0.1  mm.  The  data  acquisition  system 
was  programmed  to  acquire  data  every  0.05  mm  of 
detector  movement. 

To  determine  the  relative  position  between  the  light 
source  and  the  isotropic  detector  for  optical  properties 
or  absorption  spectra  measurements,  software  was 
developed  to  automatically  reset  the  position  of  the  peak 
fluence  rate  as  the  zero  position.  This  point  corre¬ 
sponded  to  the  point  with  the  shortest  distance  between 


the  line  defined  by  the  detector  catheter  and  the  point 
source.  A  similar  scan  was  used  for  a  white  light  source 
to  determine  the  relative  distances  between  the  light 
source  and  the  locations  where  the  diffuse  absorption 
spectra  were  taken,  typically  at  x  =  -0.4,  -0.2,  0,  0.2, 
0.4,  0.6,  0.8,  1.0,  1.2,  1.4,  and  1.6  cm.  The  position  of 
the  light  source  was  defined  as  the  distance  of  the  source 
from  the  tip  of  catheter,  intended  to  be  placed  at  the 
apex  of  the  prostate  under  ultrasound  guidance. 


2.3.  Measurement  of  optical  properties  at  the  treatment 
wavelength 


The  transport  scattering  (p[)  and  absorption  (pa) 
coefficients  characterize  the  scattering  and  absorption 
properties  of  tissue.  Within  the  diffusion  approximation, 
the  light  fluence  rate  </>  at  a  distance  r  from  a  point 
source  of  source  strength,  S ,  can  be  expressed  as  [28] 


4>  = 
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where  S  is  the  strength  of  the  point  source;  is  the 
fluence  rate  at  position  r;  /refr  =  ^/3  •  pa  ■  p's  is  the  effec¬ 
tive  attenuation  coefficient  in  tissues  and  is  applicable 
for  a  wider  range  of  pa  and  p's  combinations  than  the  tra¬ 
ditional  definition  of  /refr  =  y/3  ■  pa  ■  (p's  +  jua)  [29], 
r  —  \]x2  +  h2,  where  x  and  h  are  the  parallel  and  perpen¬ 
dicular  distance  from  the  center  of  the  point  source  (see 
Fig.  1(b)).  The  two  free  parameters  (pa  and  p!s)  are  inher¬ 
ently  separable  in  that  for  a  point  source  with  a  given 
source  strength,  the  magnitude  of  the  fluence  rate  near 
the  light  source  ( x  =  0)  is  determined  by  p[  only  and 
the  slope  of  the  spatial  decay  of  the  light  fluence  rate 
is  determined  by  /4ff  only. 

In  theory,  measurements  of  <j>  at  two  different  dis¬ 
tances  r  from  the  point  source  with  source  strength  S 
are  sufficient  to  determine  both  pa  and  Measure¬ 
ments  at  multiple  sites  allow  evaluating  the  variation 
of  these  optical  characteristics  within  the  prostate  vol¬ 
ume.  Since  Eq.  (1)  is  a  non-linear  equation  of  two  free 
parameters  pa  and  p!s,  we  used  the  differential  evolution 
algorithm  developed  by  Storn  et  al.  [30],  This  algorithm 
is  simple,  and  converges  faster  and  is  more  robust  than 
adaptive  simulated  annealing  or  the  annealed  Nelder 
and  Mead  approach  [30],  We  modified  the  algorithm 
to  require  that  all  free  parameters  (pa  and  /«'  in  this  case) 
are  positive  [31].  The  deviations  between  measurement 
and  fit  are  represented  by  standard  deviation, 

fib  —  meas,)2/(n  —  1),  where  n  is  the  number 
of  data  points  involved  in  the  fit. 


2.4.  Calibration  of  the  isotropic  detector 

Isotropic  detectors  described  by  Marijnissen  and  Star 
[32]  were  used  as  detectors  in  this  study.  Each  isotropic 
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Fig.  1.  (a)  Photograph  and  (b)  Schematic  of  the  measurement  geometry,  illustrating  the  setup  for  measurement  of  distribution  of  light  fluence  rate, 
optical  properties,  and  diffuse  absorption  spectra.  Two  step  motors  are  used,  one  for  the  light  source  and  the  other  for  the  detector.  The  distribution 
of  light  fluence  rate  was  determined  by  moving  the  detector  in  a  catheter  during  initial  PDT  treatment.  The  experimental  setup  for  the  optical 
properties  and  absorption  spectra  measurement  is  identical  except  a  different  light  source  (732  nm  vs.  white  light)  was  used.  The  distribution  of 
optical  properties  and  absorption  spectra  was  determined  by  moving  both  a  point  source  (by  z  from  the  end  of  catheter)  and  a  detector  (by  x  from  the 
point  source  location)  before  and  after  PDT  treatment.  The  fluorescence  distribution  is  achieved  using  a  single  side  firing  fiber  used  as  both  a  light 
source  and  light  detector,  (c)  3D  graph  of  the  arrangement  of  catheters  in  a  typical  prostate  treatment.  The  left  upper  quadrant’s  treatment  (blue)  and 
detector  (red)  catheters  are  shown  in  color  for  emphasis.  The  remaining  catheters  are  for  treatment  of  the  other  quadrants.  The  detector  catheters  for 
these  quadrants  are  not  shown. 


detector  was  made  of  an  optical  fiber  with  a  spherical  tip 
made  of  Ti02  (a  scattering  material).  The  isotropic 
detectors  were  made  by  Rare  Earth  Medical  (now  Car- 
dioFocus,  Norton,  MA)  and  had  an  isotropy  of  better 
than  ±30%  from  any  direction  except  for  angles  within 
30°  of  the  optical  fiber  attachment  point.  Each  detector 
fiber  was  connected  to  a  photodiode  via  a  SMA  connec¬ 
tor.  The  measured  photovoltage  (V)  from  the  isotropic 
detector  was  converted  to  light  fluence  rate  using: 

4>  =  a(V-b),  (2) 


where  a  is  the  conversion  factor  and  b  characterizes  the 
leakage  of  the  photodiode.  This  calibration  was  per¬ 
formed  under  collimated  732-nm  laser  light  in  air.  When 
the  isotropic  detector  was  used  to  measure  fluence  rate 
in  tissue  (or  in  a  liquid  optical  phantom),  a  tissue  (or 
water)  correction  factor  of  2.0  was  used  [32].  In  our  case, 
the  1  mm  inner  diameter  catheter  is  filled  with  air,  but 
the  detector  tissue  correction  factor  is  still  2.0  as  if  the 
air  is  replaced  by  the  outer  most  medium  (tissue).  The 
calibration  of  the  isotropic  detectors  was  checked  to 
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be  accurate  to  within  5%  using  an  integrating  sphere  be¬ 
fore  each  individual  measurement. 

A  12-channel  light  dosimetry  system  developed  in  the 
Department  of  Radiation  Oncology  at  the  University  of 
Pennsylvania  was  used  for  all  in  situ  dosimetry  measure¬ 
ments.  Five  different  isotropic  detectors  were  used  and 
the  conversion  factors  were  a  =  129-161  and  a  =  63-72 
mW/cm2/V  for  the  four  isotropic  detectors  with  0.5 
mm  scattering  tips  and  the  isotropic  detector  with  1 
mm  scattering  tip  used  in  the  urethra,  respectively.  For 
all  isotropic  detectors,  b  was  found  to  be  0.020  V.  The 
response  was  found  to  be  linear  over  the  light  fluence 
rate  range  (0-1200  mW/cm2)  to  within  5%. 

2.5.  Absorption  and  fluorescence  spectroscopy 
measurements 

To  obtain  absorption  spectra,  a  method  similar  to  the 
optical  property  measurement  described  above  was 
used.  In  this  case,  the  laser  and  photodiode  detector 
were  replaced  with  a  white  light  source  and  a  spectro¬ 
graph,  respectively.  To  calibrate  the  detector,  the  detec¬ 
tion  fiber  and  source  fiber  were  both  placed  in  the 
integrating  sphere,  and  a  reference  spectrum  was  ob¬ 
tained  by  the  CCD.  The  fluence  rate  measured  at  each 
wavelength  was  related  to  the  source  power  at  that 
wavelength  by  a  constant  factor  determined  by  the 
geometry  of  the  integrating  sphere.  This  factor  was  cal¬ 
culated  based  on  independent  measurements  at  three 
wavelengths  using  a  calibrated  detector.  With  this  cali¬ 
bration,  the  detector  measures  the  ratio  of  fluence  rate 
to  source  power.  The  fluence  rate  spectra  between  650 
and  800  nm  were  fit  using  a  nonlinear  fitting  algorithm 
under  the  assumption  that  the  reduced  scattering  spec¬ 
trum  had  the  form  /('  =  A/Tb,  where  A  is  expressed  in 
nm.  The  free  parameters  in  the  fit  were  A,  b,  and  the 
absorption  coefficients  at  8  selected  wavelengths.  Once 
A  and  b  were  determined  by  this  fit,  the  absorption  coef¬ 
ficients  at  the  remaining  wavelengths  were  determined 
by  a  second  fit  in  which  the  value  of  p's  was  fixed  at  its 
predetermined  value  for  each  wavelength  [25].  The 
measurement  uncertainty  of  the  absorption  spectra  has 
been  examined  in  optical  phantoms  made  of  intralipid, 
MLu,  and  black  ink,  each  with  known  absorption  spec¬ 
tra.  A  comparison  between  the  absorption  spectra  ob¬ 
tained  by  our  measurement  and  analysis  method  and 
the  absorption  spectrum  reconstructed  from  individual 
spectra  of  the  phantom  components  shows  that  our 
algorithm  can  reconstruct  pa  with  an  uncertainty  of 
approximately  5%  for  the  wavelength  range  where 
absorption  is  appreciable,  in  the  phantom  case  between 
650-800  nm  (data  not  shown). 

Fluorescence  spectra  were  acquired  using  an  optical 
fiber  with  a  beveled  tip,  which  emits  and  collects  light 
at  right  angles  to  its  optical  axis.  Excitation  light  was 
provided  by  a  460-nm  light-emitting  diode  (LED),  and 


passed  through  a  dichroic  beamsplitter,  which  directed 
the  fluorescence  collected  by  the  same  fiber  back  into 
the  spectrograph. 

Both  the  absorption  and  fluorescence  spectra  were 
analyzed  using  the  singular  value  decomposition  fitting 
algorithm  developed  by  Finlay  et  al.  [33]  to  determine 
the  contributions  to  the  spectra  of  known  absorbers  or 
fluorophores.  In  addition  to  the  spectra  of  known  com¬ 
ponents,  the  basis  set  included  a  61 -term,  exponentially 
weighted  Fourier  series  designed  to  account  for  the  pres¬ 
ence  of  unknown  absorbers  or  fluorophores  [33].  In  the 
case  of  absorption  spectra,  the  known  spectra  were 
those  of  oxy-  and  deoxyhemoglobin,  MLu  and  water. 
Because  the  absorption  coefficient  at  each  wavelength 
is  determined  absolutely,  the  absolute  concentrations 
of  the  various  absorbers  can  be  determined  quantita¬ 
tively  from  their  contributions  to  the  complete  absorp¬ 
tion  spectrum. 

In  the  fluorescence  case  the  basis  spectra  were  the  flu¬ 
orescence  of  MLu  and  the  background  fluorescence 
originating  in  the  catheter.  The  fitting  algorithm  allows 
the  separation  of  these  two  components,  allowing  us 
to  determine  an  MLu  contribution  free  from  back¬ 
ground  contamination.  The  MLu  fluorescence  was  nor¬ 
malized  by  dividing  by  the  catheter  background  to 
account  for  variations  in  lamp  output.  A  single  conver¬ 
sion  factor  of  12.9  mg/kg  between  normalized  fluores¬ 
cence  signal  at  the  peak  wavelength  (unitless)  and 
MLu  concentration  (mg/kg)  was  established  by  compar¬ 
ing  the  signal  obtained  in  one  prostate  in  vivo  to  the 
MLu  concentration  at  the  same  position  determined 
by  absorption  spectroscopy.  Once  this  factor  was  deter¬ 
mined,  the  fluorescence  could  be  analyzed  independent 
of  absorption  measurements  in  each  quadrant. 

2.6.  Photodynamic  therapy  and  in-vivo  light 
measurements 

MLu  (Pharmacyclics,  Inc.,  Sunnyvale,  CA)  was 
administered  at  a  dose  of  2  mg/kg  intravenously  3  h 
prior  to  light  administration  [12,13].  This  drug-light 
interval  was  chosen  because  preclinical  studies  in  other 
model  systems  demonstrated  the  greatest  antitumor  effi¬ 
cacy  with  this  timing  [12,34].  A  15-W  diode  laser,  model 
730  (Diomed,  Ltd.,  Cambridge,  United  Kingdom)  was 
used  as  the  732  nm  light  source. 

Interstitial  CDFs  were  placed  in  the  gland  using  a 
template  with  evenly  spaced  holes,  which  was  attached 
to  the  TRUS  unit.  For  each  light  source,  a  17-gauge 
plastic  catheter  (Flexi-needle  from  Best  Industries, 
Inc.,  Springfield,  VA)  containing  a  metal  trocar  was 
placed  through  the  template  and  into  the  prostate.  The 
trocar  was  removed  and  replaced  with  the  light  diffuser. 

The  light  energy  delivered  was  prescribed  based  on  in 
situ  measured  light  fluence.  Each  patient  received  a  light 
fluence  between'  50-100  J/cm2.  The  maximum  unit 
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length  source  strength  in  any  one  fiber  was  limited  to 
150  mW/cm.  Measurements  were  taken  at  the  above 
mentioned  distances  from  the  light  source  before  and 
after  light  treatment  in  the  four  quadrants  of  each  pros¬ 
tate.  The  light  sources  were  moved  along  the  catheter  in 
several  locations  to  quantify  the  variation  of  optical 
properties  in  the  prostate  gland  along  the  catheter. 


3.  Results 

Table  1  summarizes  the  treatment  parameters  and 
the  measurements  done  for  the  two  patients.  We  have 
made  four  types  of  measurements:  (1)  Distribution  of 
light  fluence  rate;  (2)  Distribution  of  optical  properties 
at  treatment  wavelength;  (3)  Distribution  of  absorp¬ 
tion  spectra;  (4)  Distribution  of  fluorescence  spectra. 
Tissue  oxygenation  and  MLu  concentration  can  be  ex¬ 
tracted  from  absorption  spectra  in  vivo.  MLu  concen¬ 
tration  in  vivo  is  also  extracted  from  the  fluorescence 
spectra. 

3.1.  Distribution  of  fluence  rate 

Fig.  2  shows  the  light  fluence  rate  measured  in  three 
quadrants  during  PDT  treatment  of  patient  13.  Three 
curves  are  shown.  The  first,  labeled  ‘RUQ’  was  taken 
by  a  detector  in  the  right  upper  quadrant  of  the  prostate 
during  interstitial  illumination  of  that  quadrant.  The 
light  fluence  rate  for  RUQ  (solid  line)  between  1.5  and 
2.5  cm  was  about  5  times  as  high  as  the  rest  of  region 
(e.g.,  between  3  and  4.5  cm).  This  was  caused  by  the 
loading  pattern  of  the  light  sources,  in  which  two  1-cm 
CDF’s  were  placed  at  the  apex  of  the  prostate.  The  sec¬ 
ond  curve,  labeled  ‘LLQ’  shows  the  light  fluence  rate 
that  was  measured  in  the  left  lower  quadrant  during  illu¬ 
mination  of  the  left  upper  quadrant,  and  demonstrates 
the  possible  extent  of  light  spread  between  quadrants. 
The  third  curve  labeled  ‘RLQ’  shows  light  fluence  rate 
that  was  taken  in  the  right  lower  quadrant  during  treat¬ 
ment  of  that  quadrant.  The  symbols  correspond  to  the 
predicted  light  fluence  rate  from  a  point  source  with 
source  strength  of  35  mW  based  on  the  optical  proper¬ 
ties  measured  along  the  detector  catheter.  (The  use  of 
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Fig.  2.  Distribution  of  light  fluence  rate  for  732  nm  light  in-vivo 
during  PDT.  The  profile  measured  in  the  LLQ  was  acquired  during 
illumination  of  the  LUQ  and  indicates  the  spread  of  light  between 
adjacent  quadrants.  The  open  circles  indicate  the  fluence  rate  at  0.5  cm 
from  the  source  predicted  by  the  diffusion  theory  using  the  optical 
properties  measured  at  various  points  in  the  RLQ.  The  geometry  of 
each  measurement  is  shown  in  the  diagrams  below  the  figure.  In  each 
case,  an  end-view  of  the  prostate  is  shown.  The  measurement  channel 
in  each  case  is  marked  by  an  ‘x’,  and  the  filled  circles  indicate  the 
channels  delivering  illumination.  The  length  in  cm  of  the  cylindrical 
diffuser  in  each  channel  is  indicated  in  the  diagram. 


a  cylindrical  theoretical  model  is  not  feasible  due  to 
the  varying  optical  properties  along  the  catheters,  thus 
point  sources  are  used.)  The  corresponding  loading  pat¬ 
terns  of  CDF’s  (‘o’)  and  the  locations  of  the  detector  (‘x’) 
are  shown  in  the  insert.  The  number  next  to  each  source 
channel  indicates  the  length  in  cm  of  the  corresponding 
CDF.  All  linear  light  sources  started  from  the  apex  of 
the  prostate  (z  =  0  cm). 

We  did  not  measure  significant  light  fluence  rate  in 
urethra  for  the  12  patients  treated  so  far. 


Table  1 

Variation  of  tissue  optical  properties,  tissue  concentration  of  MLu,  and  tissue  oxygenation  among  the  two  patients  before  and  after  PDT 


Patient  number 

Light  fluence  (J/cm2) 

Measurement  conditions 

#<u  (cm  ') 

#4  (cm  *) 

f‘etr  (cm  ') 

MLua  (ng/mg) 

[Hb],  (pM) 

sto2  (%) 

12 

50 

Before  PDT 

0.3-1 .6 

1.2-18 

1. 9-5.4 

2.5-5.0 

_ 

_ 

After  PDT 

0.25-0.5 

5.9-39 

2.5-5.0 

0.4-1. 7 

- 

- 

13 

100 

Before  PDT 

0.11-0.9 

7.5—40 

2.2-6.3 

1. 4-9.2 

51-310 

79-88 

After  PDT 

0.13-0.6 

5.5-34 

2.2~5.2 

1. 1-8.1 

55-370 

67-89 

The  range  specifies  the  minimum  and  maximum  values  of  the  measured  quantities  at  different  locations.  Both  patients  received  2  mg/kg  MLu 
intravenously  at  a  drug-treatment  interval  of  3  h. 

a  For  patient  #12,  the  MLu  concentration  is  determined  from  the  fluorescence  measurement. 
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3.2.  Optical  properties  at  treatment  wavelength 

A  typical  profile  scan  from  a  point  source  is  shown  in 
Fig.  3.  The  data  were  fitted  using  Eq.  (1)  to  obtain  the 
optical  properties.  The  optical  properties  derived  from 
both  patients  are  shown  in  Fig.  4.  These  data  show 
the  variation  of  measured  absorption  coefficients  and 
effective  absorption  coefficient  vs.  positions  along  the 
catheters  in  the  prostate  gland  of  patient  #12  (a  and 
b)  and  patient  #13  (c  and  d).  Due  to  time  constraints 
we  did  not  measure  LUQ  after  PDT  for  patient  #12. 
We  did  not  get  any  results  (before  and  after  PDT)  for 
LLQ  of  patient  #12  and  #13  due  to  bleeding.  The  effec¬ 
tive  attenuation  coefficients  varied  between  1.9  and  6.3 
cm-1  while  the  absorption  coefficients  varied  between 
0.1  and  1.6  cm-1.  As  a  result  of  the  heterogeneity  of 
optical  properties,  the  light  fluence  rates  per  unit  source 
strength  at  0.5  cm  from  the  point  source  varied  between 
0.2  and  0.6  cm-2  between  the  two  patients. 

3.3.  Absorption  spectra 

A  typical  absorption  spectrum  acquired  in  vivo  is 
shown  in  Fig.  5(a),  along  with  the  components  of  the 
spectrum  as  determined  by  linear  fitting.  The  dominant 
absorbers  in  the  wavelength  region  of  our  measurement 
are  hemoglobin  in  its  oxygenated  and  deoxygenated 
forms,  MLu,  and,  to  a  lesser  extent,  water.  In  the  cases 
presented  here,  the  contribution  of  the  Fourier  synthesis 
was  smaller  than  the  contributions  of  known  absorbers, 
indicating  that  the  majority  of  absorption  in  tissue  was 
accurately  accounted  for  by  our  basis  set.  From  the  con¬ 
tributions  of  oxy-  and  deoxyhemoglobin  to  the  meas- 


Fig.  3.  Measured  light  fluence  rate  per  unit  source  strength  (0/S)  at 
732  nm  versus  distances  along  the  catheter,  x,  from  the  point  source 
measured  in-vivo  in  human  prostate  gland  for  patient  #13.  Line  is 
measured  data  and  symbols  are  fits.  The  optical  properties  are:  o  - 
/<„  =  0.23  cm-1,  /i'  =  7.3  cm-1,  <j>(0.5)IS  =1.1  cm-2,  h  =  0.5.  Similar 
measurements  were  made  for  different  wavelengths  for  the  absorption 
spectra,  although  the  measurements  are  now  made  at  only  1 1  points 
between  -0.4  and  1.6  cm  with  0.2  cm  steps. 


ured  absorption  spectrum,  we  determined  the  total 
hemoglobin  concentration  [Hb]t,  the  sensitizer  concen¬ 
tration  [MLu],  and  the  hemoglobin  oxygen  saturation 
St02  given  by 

St02  =  [Hb02]/([Hb]  +  [HbOJ).  (3) 

The  values  of  these  three  parameters  as  functions  of 
position  within  a  typical  prostate  quadrant  (RUQ,  pa¬ 
tient  13),  are  plotted  in  Fig.  6(a).  The  concentration  of 
MLu  is  given  in  ng  mg-1  and  has  been  scaled  by  a  factor 
of  10  for  comparison  with  the  hemoglobin  concentration 
and  saturation. 

3.4.  Fluorescence  spectra 

The  concentration  of  MLu  via  its  intrinsic  fluores¬ 
cence  emission  was  also  measured  around  750  nm.  A 
typical  fluorescence  spectrum  is  shown  in  Fig.  5(b). 
The  singular  value  decomposition  fitting  algorithm  is 
able  to  separate  the  component  arising  from  the  back¬ 
ground  fluorescence  of  the  catheter  and  fiber  assembly 
(labeled  ‘Bkgnd’)  from  the  MLu  component.  The  small 
amplitude  of  the  Fourier  component  indicates  that  the 
fluorescence  is  dominated  by  these  two  contributions. 
The  background  fluorescence  provides  a  measure  of 
the  lamp  intensity,  and  is  used  to  normalize  the  MLu 
fluorescence.  The  MLu  distribution  determined  by  fluo¬ 
rescence  spectroscopy  is  overlaid  on  the  corresponding 
distribution  determined  by  absorption  spectroscopy  in 
Fig.  6(b).  The  spatial  distribution  of  MLu  determined 
by  the  two  methods  is  similar. 


4.  Discussion 

The  main  objective  of  this  work  was  to  demonstrate 
the  feasibility  of  measuring  the  distribution  of  important 
dosimetric  parameters  for  PDT  in  vivo,  namely  the  tis¬ 
sue  optical  properties,  tissue  concentration  of  drug, 
and  tissue  oxygenation.  Measurements  were  made 
before,  during  and  after  PDT.  Differences  within  and 
between  patients  were  also  evaluated. 

4.1.  Distribution  of  light  fluence  rate  in  vivo 

Our  data  show  that  there  is  heterogeneity  of  light  flu¬ 
ence  in  different  regions  of  the  prostate  (Fig.  2).  The 
data  also  show  that  light  delivery  to  one  quadrant  of 
the  prostate  may  lead  to  the  delivery  of  measurable  light 
fluence  to  other  regions  of  the  prostate.  For  the  region 
of  high  fluence  rate,  there  were  three  CDFs  contributing 
to  light  fluence  rate,  due  to  the  contributions  from  the 
two  1-cm  long  CDFs  and  the  one  4-cm  long  CDF. 
For  the  region  of  low  fluence  rate,  only  the  4-cm  long 
linear  fiber  was  contributing  to  the  light  fluence  rate. 
A  similar  explanation  can  also  be  used  for  the  fluence 
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Fig.  4.  In-vivo  distribution  of  (a)  absorption  and  (b)  effective  attenuation  coefficients  at  732  nm  in  the  human  prostate  for  patient  #12.  (c)  absorption 
and  (d)  effective  attenuation  coefficients  at  732  nm  in  the  human  prostate  for  patient  #13. 


rate  distribution  in  the  LLQ.  The  light  fluence  rate  for 
the  RLQ  was  much  larger  than  that  for  the  RUQ,  which 
can  be  explained  by  dilferences  in  optical  properties.  In 
the  RLQ,  two  similar  length  CDF’s  were  used  (3  and  4 
cm).  If  the  tissue  optical  properties  were  uniform  along 
the  catheter,  the  profile  would  have  been  uniform.  How¬ 
ever,  the  profile  of  light  fluence  rate  is  not  uniform  be¬ 
cause  the  optical  properties  change  along  the  catheter. 
The  predicted  light  fluence  rate  for  a  35  mW  source 
based  on  measured  optical  properties  (#a  and  is 
indicated  by  symbols.  The  similarity  in  shape  between 
the  measured  and  the  predicted  light  fluence  rate  indi¬ 
cates  the  variation  of  light  fluence  rate  was  due  to  the 
variation  of  optical  properties. 

4.2.  Distribution  of  optical  properties  at  treatment 
wavelength 

Fig.  3  shows  a  typical  example  of  measured  fluence 
rate  per  source  strength  from  a  point  source.  The  meas¬ 
ured  profile  usually  contains  800  data  points  at  0.05  mm 


steps.  The  fitting  assumes  that  the  distance  h  between  the 
source  and  detector  catheter  is  also  a  free  variable, 
which  is  allowed  to  vary  up  to  0.2  cm  from  the  known 
separation  determined  from  the  template  positions.  We 
used  the  same  method  to  obtain  the  absorption  spectra 
from  a  white  light  source  using  a  somewhat  longer  step 
size  (0.2  cm)  in  the  range  of  —0.4  to  1.6  cm  relative  to 
the  point  source. 

The  intra-prostatic  tissue  heterogeneity  of  the  optical 
properties  at  732  nm  is  demonstrated  in  Fig.  4.  For  the 
same  patient,  the  effective  attenuation  coefficient,  pe ff, 
varied  by  up  to  3  times  between  different  quadrants  of 
the  same  prostate  (Fig.  4(d)).  This  large  variation  of 
optical  properties  resulted  in  a  large  difference  in  light 
fluence  rate  between  the  right  upper  quadrant  (RUQ) 
and  the  right  lower  quadrant  (RLQ)  (see  Fig.  2).  Within 
the  same  quadrant  of  a  prostate,  /rerr  can  change  sub¬ 
stantially  as  well.  Since  each  CDF  has  uniform  light 
strength  along  the  catheter,  the  variation  of  optical 
properties  can  result  in  large  variation  of  light  fluence 
rate  along  the  catheter,  as  demonstrated  by  the  case  of 
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Fig.  5.  Components  of  (a)  absorption  spectra  and  (b)  fluorescence 
spectra  acquired  prior  to  PDT  from  the  LUQ  of  patient  #13.  In  each 
case,  the  measured  data  are  labeled  as  symbols  (‘o’)  and  the 
components  are  labeled  as  lines.  The  components  include:  MLu, 
water  (H20),  deoxyhemoglobin  (Hb),  and  hemoglobin  (HbCL).  The 
‘Fourier’  component  is  a  Fourier  series  designed  to  account  for 
unknown  absorbers  or  fluorophores. 


Fig.  6.  In  vivo  distribution  of  (a)  St02,  blood  volume  (pM),  and  MLu 
concentration  determined  using  the  absorption  spectra  and  (b)  MLu 
concentration  as  determined  by  absorption  spectra  (triangles)  and 
fluorescence  (circles)  measurements  for  RUQ  in  patient  13. 


the  RLQ.  Since  our  model  assumes  the  optical  proper¬ 
ties  of  the  medium  to  be  homogeneous,  the  measured 
optical  properties  has  limitation  in  that  it  should  be  an 
average  of  tissue  optical  properties  within  the  maximum 
distance  between  the  source  and  the  detector,  typically  2 
cm  (see  Fig.  3). 

The  measurement  standard  deviation  of  pa,  p's  and  pea 
is  7%,  20%,  and  5%,  respectively,  when  h  is  optimized  in 
the  fitting  [35],  Details  of  the  error  analysis  are  included 
in  [35], 

Since  our  model  assumes  infinite  medium,  we  make 
sure  that  the  measurement  points  are  at  least  1  cm  away 
from  the  boundaries  of  prostate  to  minimize  the  bound¬ 
ary  effect.  The  starting  position  of  the  point  source  is 
moved  1  cm  from  the  end  of  the  catheter  before  optical 
properties  measurement.  The  reduction  of  the  light  flu- 
ence  rate  at  the  tissue  boundary  has  been  characterized 
in  a  solid  prostate  phantom  to  be  less  than  2%  at  1  cm 
from  the  boundary  at  732  nm  for  the  range  of  optical 
properties  used  in  the  study  [36], 


We  have  demonstrated  in  canine  prostate  that  the 
absorption  coefficient  at  732  nm  is  approximately  pro¬ 
portional  to  the  tissue  concentration  of  MLu  [19].  This 
linear  relationship  is  assumed  to  hold  for  human  pros¬ 
tate  as  well.  Since  pa  varied  by  up  to  2.5  times  along 
some  catheters  (Fig.  4(c),  RUQ,  after  PDT),  the  tissue 
concentration  of  MLu  should  vary  by  the  same  magni¬ 
tude  along  the  catheter.  This  is  consistent  with  the  distri¬ 
bution  of  MLu  concentration  measured  by  absorption 
spectra  and  fluorescence  in  the  same  location  (see  Fig. 
6(b)). 

4.3.  Distribution  of  absorption  and  fluorescence  spectrum 

Fig.  5(a)  illustrates  how  one  can  extract  the  concen¬ 
trations  of  oxyhemoglobin  (Hb02),  deoxyhemoglobin 
(Hb),  and  MLu  from  the  measured  absorption  spec¬ 
trum.  Because  the  absorption  spectrum  is  determined 
by  fitting  data  using  Eq.  (1)  at  each  wavelength 
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independently  without  the  knowledge  of  the  known 
spectral  components  [25],  and  because  the  sum  of 
known  components  accurately  fits  the  extracted 
spectrum  with  little  absorption  accounted  for  by  the 
Fourier  components,  we  are  very  confident  of  extracted 
MLu,  Hb02  and  Hb  concentrations.  Since  the  average 
optical  penetration  depth  (<5  =  l/,ueff)  at  732  nm  is  about 
0.4  cm  in  human  prostate  and  the  separation  between 
the  detector  and  the  light  source  (0.5-1. 5  cm)  is  gener¬ 
ally  larger  than  <5,  the  extracted  quantities  reflect  a  mac¬ 
roscopic  average  in  a  tissue  dimension  of  0.5-1. 5  cm, 
depending  on  the  separation  between  the  detector  and 
the  light  source.  Fig.  5(b)  shows  the  corresponding  anal¬ 
ysis  of  a  typical  fluorescence  spectrum.  The  MLu  con¬ 
centration  obtained  from  each  spectrum  is  normalized 
to  the  background  signal  and  therefore  the  measured 
MLu  concentration  is  insensitive  to  variations  in  lamp 
intensity  and  fiber  coupling  efficiency.  The  inherent 
absorption  by  tissue  is  much  greater  at  the  excitation 
wavelength  (460  nm)  than  the  wavelengths  over  which 
we  measure  absorption  spectra.  The  generation  of  fluo¬ 
rescence  is  therefore  confined  to  a  smaller  region  of  tis¬ 
sue  than  that  sampled  by  the  absorption  spectra 
measurement. 

Fig.  6  illustrates  the  variation  in  total  hemoglobin 
concentration,  hemoglobin  saturation,  and  MLu  con¬ 
centration  within  one  quadrant  of  a  typical  prostate. 
The  variation  in  measured  [Hb],  and  [MLu]  likely  re¬ 
sult  from  variations  in  vascular  density  or  perfusion, 
limiting  the  supply  of  drug  and  blood  to  some  regions 
of  the  tissue.  The  agreement  in  the  shape  of  the  distri¬ 
butions  of  MLu  determined  by  absorption  and  fluores¬ 
cence  spectroscopy  demonstrates  the  ability  of  both 
methods  to  give  quantitative  information  concerning 
the  distribution  of  drug  in  vivo.  While  the  two  methods 
give  similar  information,  each  has  its  advantages.  The 
absorption  spectra  measurement  gives  not  only  the 
drug  distribution  but  hemodynamic  information  as 
well.  The  fluorescence  measurement,  on  the  other  hand, 
samples  a  much  smaller  volume  of  tissue  than  the 
absorption  spectra  measurement,  and  samples  many 
more  points  than  is  practical  for  absorption  spectra 
measurements,  allowing  a  higher-resolution  map  of 
the  spatial  drug  distribution. 

Table  1  shows  the  variation  of  optical  properties  be¬ 
fore  and  after  PDT.  For  patient  12,  we  have  observed 
substantial  reduction  of  MLu  tissue  concentration  after 
PDT  treatment,  presumably  due  to  photobleaching.  A 
smaller  reduction  of  MLu  tissue  concentration  was  ob¬ 
served  for  patient  13.  The  absorption  coefficient,  which 
has  been  shown  to  be  approximately  linear  to  MLu  tis¬ 
sue  concentration,  also  decreased  after  PDT  treatment, 
more  for  patient  12  than  for  patient  13.  The  effective 
attenuation  coefficient  also  varied  after  PDT,  although 
the  magnitude  of  the  change  was  generally  much 
smaller. 


5.  Conclusions 

We  have  demonstrated  inter-  and  intra-prostate  var¬ 
iation  of  optical  properties  and  MLu  tissue  concentra¬ 
tion.  The  variation  of  optical  properties  can  be  used  to 
explain  the  observed  variation  of  light  fluence  rate.  We 
observed  that  the  optical  penetration  depth  varied  be¬ 
tween  0.15-0.5  cm  for  732  nm  light  within  one  prostate. 
The  tissue  concentration  of  MLu  varied  between  1.1  and 
8  ng/mg  within  one  prostate.  These  studies  confirm  sub¬ 
stantial  inter-organ  and  intra-organ  variations  in  optical 
properties  and  drug  concentration  in  the  prostates.  Gi¬ 
ven  this  heterogeneity,  a  real-time  dosimetry  measure¬ 
ment  and  feedback  system  for  monitoring  light 
fluences  during  treatment  should  be  considered  for  inter¬ 
stitial  prostate  PDT  studies. 
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Abstract 

We  have  developed  a  method  to  quickly  determine  tissue  optical  properties 
(absorption  coefficient  /xa  and  transport  scattering  coefficient  fi's)  by  measuring 
the  ratio  of  light  fluence  rate  to  source  power  along  a  linear  channel  at  a  fixed 
distance  (5  mm)  from  an  isotropic  point  source.  Diffuse  light  is  collected  by 
an  isotropic  detector  whose  position  is  determined  by  a  computer-controlled 
step  motor,  with  a  positioning  accuracy  of  better  than  0.1  mm.  The  system 
automatically  records  and  plots  the  light  fluence  rate  per  unit  source  power  as  a 
function  of  position.  The  result  is  fitted  with  a  diffusion  equation  to  determine 
Ma  and  jtt'.  We  use  an  integrating  sphere  to  calibrate  each  source-detector 
pair,  thus  reducing  uncertainty  of  individual  calibrations.  To  test  the  ability 
of  this  algorithm  to  accurately  recover  the  optical  properties  of  the  tissue,  we 
made  measurements  in  tissue  simulating  phantoms  consisting  of  Liposyn  at 
concentrations  of  0.23,  0.53  and  1.14%  (/x'  =  1.7-9. 1  cm-1)  in  the  presence 
of  Higgins  black  India  ink  at  concentrations  of  0.002,  0.012  and  0.023% 
(Ma  =  0.1-1  cm-1).  For  comparison,  the  optical  properties  of  each  phantom 
are  determined  independently  using  broad-beam  illumination.  We  find  that  /xa 
and  /x'  can  be  determined  by  this  method  with  a  standard  (maximum)  deviation 
of  8%  (15%)  and  18%  (32%)  for  /xa  and  /x',  respectively.  The  current  method 
is  effective  for  samples  whose  optical  properties  satisfy  the  requirement  of  the 
diffusion  approximation.  The  error  caused  by  the  air  cavity  introduced  by  the 
catheter  is  small,  except  when  /xa  is  large  (/xa  >  1  cm-1).  We  presented  in  vivo 
data  measured  in  human  prostate  using  this  method. 

(Some  figures  in  this  article  are  in  colour  only  in  the  electronic  version) 


1.  Introduction 

The  in  vivo  determination  of  tissue  optical  properties  has  been  an  area  of  extensive  research. 
The  optical  properties  include  determination  of  tissue  absorption  coefficient  (/xa),  scattering 
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coefficient  ( /xs )  and  scattering  anisotropy  (g).  The  optical  absorption  and  scattering  properties 
of  tissue  can  be  used  to  calculate  the  fluence  distribution  of  light  during  light-based  treatments 
such  as  photodynamic  therapy.  In  addition,  the  optical  absorption  of  tissue  can  be  used  to 
monitor  changes  in  the  volume  and  oxygenation  of  blood  and  the  perfusion  of  tissues  (Liu  et  al 
1995,  Kienle  et  al  1996,  Hull  et  al  1999,  Rolfe  2000)  and  the  distribution  of  exogenous 
absorbers  (Doombos  et  al  1999,  Solonenko  et  al  2002).  The  scattering  coefficient  can  be 
used  to  extract  information  about  cellular  structure  of  tissue  (Mourant  et  al  1998,  Gurjar  et  al 
2001).  In  the  near-infrared  region  (X  =  600-800  nm)  where  tissue  scattering  far  exceeds  tissue 
absorption,  the  diffusion  approximation  can  be  used  (Star  1997).  Under  this  approximation, 
only  the  linear  term  of  the  anisotropy  in  radiance  is  considered.  A  reduced  scattering  coefficient 
(/x'  =  /xs(l  —  g))  is,  therefore,  sufficient  to  describe  all  the  tissue  scattering  properties. 

Many  techniques  (e.g.  relative  fluence  rate  versus  depth  under  uniform  light  illumination) 
can  determine  the  effective  attenuation  coefficient,  but  cannot  separate  the  effects  of  tissue 
scattering  and  absorption.  jxa  and  /x'  can  be  determined  separately  using  ex  vivo  measurements 
of  the  absorbing  and  scattering  components  of  the  sample;  however,  these  methods  are  not 
practical  for  determining  the  optical  properties  in  vivo.  Existing  in  vivo  methods  of  determining 
the  optical  properties  of  tissue  rely  on  measurement  of  the  diffuse  reflectance  on  tissue 
surface.  For  CW  reflectance  measurements  (Farrell  and  Patterson  1992,  Nichols  etal  1997, 
Hull  et  al  1998,  Solonenko  et  al  2002,  Doombos  et  al  1999,  Swartling  et  al  2003),  the  tissue 
is  illuminated  by  a  pencil-beam  CW  light  source  and  the  diffuse  reflectance  is  recorded  at 
different  radial  distances  from  the  source.  Time-resolved  measurements  use  subnanosecond 
pulses  from  a  laser  (Patterson  et  al  1989,  Pogue  and  Patterson  1994,  Kienle  and  Patterson 
1997a,  1997b,  Coquoz  et  al  2001,  Torricelli  et  al  2001).  After  a  pulse  passes  through  tissue, 
its  time  dispersion  can  be  measured.  In  frequency  domain  methods,  the  source  is  sinusoidally 
modulated  and  modulation  amplitude  and  phase  shift  of  the  detected  signal  are  measured  to 
obtain  information  about  optical  properties. 

The  diffuse  reflectance  techniques  outlined  above,  which  can  determine  /xa  and  /x' 
simultaneously  in  vivo,  cannot  be  used  interstitially.  Our  goal  is  to  develop  an  interstitial 
method  that  can  be  used  in  vivo  to  quickly  determine  both  the  absorption  and  the  reduced 
scattering  coefficients  of  tissue  using  a  spatial  CW  method.  For  this,  we  have  developed  a 
device  to  quickly  determine  tissue  optical  properties  by  measuring  the  ratio  of  light  fluence 
rate  to  the  source  power  along  a  linear  channel  at  a  fixed  distance  (5  mm)  from  an  isotropic 
point  source.  Diffuse  light  is  collected  by  an  isotropic  detector  whose  position  is  determined 
by  a  computer-controlled  step  motor,  with  a  positioning  accuracy  of  better  than  0.1  mm.  The 
result  is  fitted  with  a  diffusion  equation,  using  a  nonlinear  optimization  algorithm,  to  determine 
/xa  and  ytx' .  This  method  has  been  applied  to  in  vivo  optical  property  measurements  in  human 
prostate. 


2.  Broad  beam  set-up 

2.1.  Description  of  the  broad  beam  set-up 

We  determined  the  optical  properties  (/xa  and  /x')  of  each  optical  phantom  by  measuring  the 
fluence  rate  as  a  function  of  depth  under  broad  beam  illumination.  We  used  an  isotropic 
detector  manufactured  by  CardioFocus,  Inc  (West  Yarmouth,  MA)  that  consists  of  an  optical 
fibre  with  a  1-mm-diameter  spherical  scattering  bulb  at  the  tip  to  measure  the  fluence  rate  at 
different  depths.  The  light  collected  by  this  detection  fibre  was  measured  and  digitized  by 
an  in  vivo  dosimetry  system  (to  be  described  later).  The  detector  position  was  controlled  by 
a  computerized  positioning  system  (Velmex,  Inc,  East  Bloomfield,  NY).  The  measurement 
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Figure  1.  Experimental  set-up  for  broad  beam  measurements.  The  optical  properties  of  the 
Liposyn  and  ink  used  for  the  optical  phantom  were  characterized  independently  using  a  broad 
parallel  beam  incident  on  the  phantom.  We  used  a  1  mm  scattering  tip  isotropic  detector  to 
measure  the  fluence  rate  at  different  depths  in  the  phantom  using  a  step  motor  that  has  an  accuracy 
of  0.1  mm.  The  measurements  were  made  under  five  different  conditions:  one  for  Liposyn 
solution,  three  for  Liposyn  solution  plus  three  ink  concentrations  with  known  optical  properties, 
and  one  for  pure  water.  The  ratio  between  Liposyn  (or  Liposyn  with  ink)  and  water  was  fit  to 
an  exponential  function:  <p  =  kexp(— /xeffd),  where  /xeff  is  the  effective  attenuation  coefficient. 
This  ratio  eliminates  the  effect  of  inverse-square  law  and  the  detector  sensitivity  variation  between 
water  and  air  for  the  isotropic  detector. 


set-up  is  shown  in  figure  1 .  For  each  optical  phantom  (Liposyn  and  ink)  two  scans  were  made, 
under  the  same  broad  beam  illumination,  one  in  the  optical  phantom,  the  other  in  pure  water. 
The  ratio  of  the  two  scans  was  used  to  correct  for  the  detector  response  difference  between 
air  and  water.  The  simultaneous  measurements  were  made  for  nine  different  conditions:  three 
Liposyn  concentrations  for  each  of  the  three  ink  concentrations.  In  addition,  optical  properties 
for  a  Liposyn  phantom  without  ink  were  measured. 

Analysis  of  the  depth  dependence  of  fluence  rate  in  phantoms  illuminated  from  above  by 
a  broad  beam  requires  a  one-dimensional  solution  to  the  diffusion  equation.  In  this  case,  the 
ratio  between  fluence  rates  measured  in  Liposyn  with  ink  and  that  measured  in  water  was  fit 
to  an  exponential  function: 

<1 b  =  kc~^d  (1) 

where  /Xcij  is  the  effective  attenuation  coefficient,  /x'  for  the  Liposyn  and  /xa  for  water  were 
determined  from  /xeff  and  the  /xa  of  ink  deduced  from  the  known  ink  concentration.  Since  /x' 
should  be  independent  of  ink  concentration  for  a  particular  Liposyn  concentration,  we  took 
the  average  of  the  three  extrapolated  values  of  /t'  (one  for  each  ink  concentration)  as  the 
/x'  of  each  Liposyn  phantom.  The  /xa  for  different  ink  concentrations  was  then  determined 
from  the  best  fit  value  of  /xeff  and  the  mean  value  of  /x'  using  the  relation  /xcff  =  N/'3/xa/x's. 
Although  /xcff  =  ,/3/x  a/x(  is  an  approximation  of  the  expression  /xcff  =  ^3/x  a(/x(  +  /xa),  Nakai 
et  al  (1997)  have  shown  that  this  expression  is  valid  over  a  wider  range  of  optical  properties, 
especially  in  cases  of  high  absorption.  The  values  of  /xa  determined  by  this  method  and  those 
expected  based  on  the  known  ink  concentrations  agreed  to  within  8%.  The  uncertainty  in 
/xcff,  as  determined  by  propagation  of  uncertainty  based  on  the  1 .5%  uncertainty  in  our  fluence 
measurements,  is  less  than  0.1%.  The  true  uncertainty  of  this  measurement  is,  therefore,  likely 
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Table  1.  A  summary  of  ju.a  and  fx's  as  a  function  of  percentage  concentration  for  scattering  and 
absorbing  materials  used  in  the  liquid  optical  phantom.  c%  is  the  percentage  concentration  of  the 
scattering  (or  absorption)  material  in  volume  dissolved  in  water  (in  volume). 

Brand 

Wavelength  (nm) 

u's(.cm  ') 

Ma  (cm  *) 

Source 

10%  Intralipid 

633 

M's  =  (11.2±0.1I)xc% 

- 

Moes  etal  (1989) 

Nutralipid 

630 

fi's  =  17.9  xc%  +  0.56 

0.0026 

Driver  et  al  (1989) 

10%  Intralipid 

400-1100 

Us  =  ll(X/730)-2-4  x  c% 

- 

van  Staveren  et  al  (1997) 

10%  Intralipid 

460-890 

U's  =  4.4(A/730)-2-33  x  c% 

0.001-0.015 

Flock  etal  (1992) 

10%  Intralipid 

700 

Us  =  7.33  x  c% 

- 

Mourant  et  al  (1997) 

30%  Liposyn 

730 

u's  =  8.1  x  c% 

0.020 

This  study 

Higgins  Ink 

594 

- 

=  123  x  c% 

Madsen  etal  (1992) 

Higgins  Ink 

730 

- 

Uz  =  42.99  x  c% 

This  study 

limited  by  the  uncertainty  in  measuring  the  volumes  of  phantom  components.  To  minimize 
this  error,  we  used  a  pipettor  which  has  an  accuracy  of  0.01  ml  to  accurately  determine  ink 
and  Liposyn  volumes  of  2000  ml. 

2.2.  Liquid  tissue-simulating  phantom 

The  tissue  simulating  phantoms  are  made  of  separate  scattering  and  absorbing  components. 
The  phantoms  were  placed  in  a  plastic  container  that  was  painted  in  black  and  was  large  enough 
(18.2  x  14.6  x  7.7  cm3)  to  avoid  scattering  from  the  boundary.  This  type  of  phantom  has  been 
described  in  the  literature  (Madsen  et  al  1992).  The  scattering  media  used  are  phospholipid 
emulsions  (Liposyn  III,  30%  Abbott  Lab,  North  Chicago,  IL).  The  scattering  coefficient  and 
concentration  of  Liposyn  are  related  by  the  expression, 

/tx'  =8.1  x  (c%IL)  cm-1  (2) 

where  (c%IL)  is  the  lipid  concentration  used.  This  formula  was  obtained  for  the  specific 
batch  of  Liposyn  used  in  the  experiments  reported  here  using  the  broad  beam  measurement 
described  above. 

We  used  Higgins  black  India  ink  #441 8  (Higgins,  Bellwood,  IL)  as  the  absorbing  medium. 
The  absorption  coefficient  dependence  on  ink  concentration  is  given  by: 

p,a  =  42.99  x  (c%ink)  cm-1  (3) 

where  (c%ink)  is  the  ink  concentration  in  per  cent  volume.  This  formula  was  obtained  from 
transmission  measurement  of  pure  ink  diluted  in  water. 

Figure  2  shows  the  results  of  all  experiments  using  the  broad  beam  method  for  nine  tissue 
simulating  phantoms  with  Liposyn  concentrations  of:  (A)  0.23%,  (B)  0.53%  and  (C)  1.14% 
and  ink  concentrations  of:  0.002%,  0.012%  and  0.023%.  The  best  fit  scattering  coefficients 
were:  1.73  cm-1  (with  absorption  coefficients  of  0.10,  0.48  and  1.00  cm-1);  4.19  cm-1 
(with  absorption  coefficients  of  0.10,  0.49  and  0.99  cm-1)  and  9.14  cm-1  (with  absorption 
coefficients  of  0.10, 0.50  and  0.99  cm-1). 

The  optical  properties  of  liquid  tissue-simulating  phantoms  composed  of  Intralipid  and 
ink  have  been  extensively  studied  for  different  wavelengths  by  several  investigators.  A 
summary  of  the  measured  scattering  and  absorption  coefficients  measured  by  other  people 
(Moes  et  al  1989,  Driver  et  al  1989,  Flock  et  al  1992,  van  Staveren  et  al  1997,  Mourant  et  al 
1997,  Madsen  et  al  1992)  and  the  results  of  this  study  are  given  in  table  1. 
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Figure  2.  Optical  properties  characterization  using  broad  beam  method.  Results  are  shown  for 
three  tissue  simulating  phantoms  with  Liposyn  concentrations  of:  (A)  0.23%,  (B)  0.53%  and 
(C)  1.14%.  For  each  Liposyn  concentration  phantom,  we  added  three  different  concentrations 
of  black  ink:  0.002%,  0.012%  and  0.023%.  Symbols  represent  measurements  with  an  isotropic 
detector.  Solid  lines  are  the  best  fit.  The  resulting  optical  properties  were:  (A)  /a'  =  1 .78  cm-1  and 
=  0. 1 0, 0.49  and  1 .01  cm- 1 ,  (B)  /n'  =  4.36  cm" 1  and  n*  =  0. 1 0, 0.50  and  0.99  cm" 1 ,  (C)  n's  = 
9.06  cm-1  and  /ra  =  0.10, 0.50  and  1.00  cm-1.  See  text  for  details. 


3.  Interstitial  set-up 

3.1 .  Description  of  the  interstitial  set-up  for  phantom  measurement 

The  experimental  set-up  for  the  parallel-catheter  measurement  system  is  shown  in  figure  3. 
We  constructed  a  device  consisting  of  3  parallel  catheters  positioned  at  3  distances  (3,  5  and 
7  mm)  from  a  central  catheter.  Only  one  separation  is  required  to  determine  and  /r',  but 
we  examined  the  accuracy  of  the  extrapolation  for  the  three  separations  independently  to  find 
the  most  suitable  separation.  An  isotropic  point  source  was  placed  in  the  middle  catheter  and 
connected  to  a  730  nm  diode  laser  (Diomed  730,  Cambridge,  UK).  An  isotropic  detector  was 
placed  in  one  of  the  parallel  detector  channels.  This  detection  fibre  was  connected  to  a  light 
dosimetry  system  (described  below).  The  detector’s  position  was  controlled  by  computer- 
controlled  step  motors  (Velmex,  Inc  East  Bloomfield,  NY),  allowing  the  detector  to  be  moved 
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Figure  3.  (A)  Picture  of  the  optical  property  device  consisting  of  4  parallel  catheters  positioned  at 
3  different  distances  (3,  5  and  7  mm)  from  the  central  catheter.  The  light  source  is  placed  in  the 
centre  catheter,  while  the  detector  is  moved  along  each  catheter,  positioned  at  different  distances 
from  the  light  source.  (B)  Top  view  of  the  optical  property  device  pictured  in  (A).  (C)  Schematics 
of  the  light  source  and  detector  placement.  The  distance  between  the  light  source  and  the  detector 
is  h.  The  light  source  is  placed  at  a  distance  x  from  the  surface  of  the  phantom,  while  the  detector  is 
moved  along  the  catheter.  The  distance  from  the  centre  of  the  detector  to  the  point  source  is  given 
by  r  =  -Jx2  +h 2.  (D)  Diagram  of  cathteter  positioning  during  prostate  PDT.  The  cathteters  are 
placed  at  fixed  distatnces  ( h )  through  a  template  and  into  the  prostate.  The  light  source  is  placed 
in  one  of  the  cathteters  and  the  isotropic  detector  is  placed  in  the  other  catheter. 


to  different  distances  from  the  light  source.  Each  data  set  was  obtained  by  scanning  the  detector 
along  its  catheter  while  the  source  remained  fixed  and  acquiring  fluence  rate  measurements  at 
0.05  mm  intervals  along  its  movement.  The  data  acquisition  time  for  a  scan  of  10  cm  distance 
with  2000  data  points  is  about  8  s  since  the  speed  of  the  step  motor  is  12.5  mm  s_1.  The 
data  acquisition  board  has  a  maximum  data  transfer  speed  of  300  kilosample/s,  which  can  be 
adjusted  to  match  the  data  acquisition  rate  of  250  samples/s  in  the  application. 

Using  the  diffusion  approximation,  the  light  fluence  rate  4>  per  source  power  5  at  a  distance 
r  from  a  point  source  can  be  expressed  (Jacques  1998): 

*t.  —  jfeff  p-Pcfrr  _  -r*eirr  rA\ 

S  -  4 nrn,  ~  4nr 6  W 

where  S  is  the  source  power  of  the  point  source  (in  mW),  <p(r)  is  the  fluence  rate  in  mW  cm-2 
at  r.  r  =  %/r  +  h2,  where  x  and  h  are  parallel  and  perpendicular  distances  from  the  centre  of 
the  point  source  (figure  3(C)). 
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Figure  4.  Measured  light  fluence  rate  per  unit  source  strength  at  distances  along  the  catheter,  x, 
from  the  point  source  measured  in  vivo  in  human  prostate  gland.  Lines  are  measured  data  and 
symbols  are  the  fits. 


3.2.  Description  of  the  interstitial  set-up  for  in  vivo  patient  measurement 

Optical  properties  (/xa,  p's  and  5=1  /p.cff)  were  also  measured  in  1 1  patients  with  locally 
recurrent  prostate  carcinoma  using  the  interstitial  set-up.  A  template  with  evenly  spaced  holes 
5  mm  apart  was  used  for  positioning  of  the  catheters  inside  the  prostate  gland  under  ultrasound 
guidance.  One  point  source  and  one  detector  were  introduced  in  two  parallel  catheters 
(figure  3(D)).  Detectors  were  placed  at  5  or  7  mm  (h)  away  from  the  light  sources.  The 
uncertainty  of  distance  between  the  light  source  catheter  and  the  detector  catheter  (h),  which 
can  be  off  1-2  mm  from  the  position  determined  by  the  template,  can  introduce  errors  in 
determining  the  optical  properties  of  the  phantom.  Because  of  this,  the  optimization  algorithm 
is  designed  to  include  the  separation  h,  as  a  separate  fitting  parameter.  Optical  properties  of 
the  prostate  were  determined  applying  the  diffusion  theory  (equation  (4))  to  the  fluence  rates 
measured  at  several  distances  (5-50  mm)  from  the  light  source.  The  isotropic  detectors  were 
calibrated  under  collimated  732  nm  laser  light  in  air  as  described  in  section  3.3. 

Figure  4  shows  the  variation  of  the  measured  light  fluence  rate  distribution  (solid  lines) 
in  different  quadrants  of  the  prostate  and  symbols  are  the  fits.  There  are  800  measured  points 
for  each  scan  with  a  resolution  of  0.05  mm  in  the  range  of  (—2,  2)  cm.  The  measured  optical 
properties  in  this  particular  case  varied  from  0.23  to  0.4  cm-1  for  the  absorption  coefficient 
and  from  6.6  to  12  cm-1  for  the  scattering  coefficient  for  different  locations  in  prostate. 

Table  2  summarizes  the  measured  optical  properties  in  human  prostate  for  1 1  patients  (Zhu 
et  al  2005a).  Five  patients  were  measured  using  a  few  points  with  manual  positioning  and  six 
were  measured  using  a  motorized  probe  (*  indicates  motorized  probe  measurements).  The 
first  column  lists  the  absorption  coefficients,  the  second  column  lists  the  scattering  coefficients 
and  the  last  column  lists  the  optical  penetration  depth,  for  each  patient.  The  values  in  the 
parenthesis  are  the  standard  deviations  of  the  mean  values  measured  from  different  locations 
in  the  same  prostate  gland.  Whenever  no  standard  deviation  is  listed,  only  one  data  point  was 
available. 

3.3.  Calibration  of  detectors 

We  used  an  isotropic  detector  made  of  an  optical  fibre  with  a  0.5-mm-diameter  scattering  tip 
(CardioFocus  West  Yarmouth,  MA)  to  measure  the  fluence  rate.  The  fibre  is  connected  to  one 
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Table  2.  Summary  of  optical  properties  measured  in  human  prostate.  The  values  in  parentheses  are 
the  standard  deviation  of  the  mean  values  measured  from  different  locations  in  the  same  prostate. 
No  standard  deviation  is  listed  if  only  one  data  point  is  available.  (*  represents  the  measurements 
done  with  the  motorized  probe.) 


Patient  number 

Ma( cm  ’) 

cm  ') 

8  (cm) 

1 

0.09 

29.8 

0.34 

2 

0.15 

22.0 

0.31 

3 

0.43  (0.28) 

7.69  (4.76) 

0.41  (0.14) 

4 

0.21 

11.8 

0.37 

5 

0.27  (0.27) 

10.5(11.2) 

0.50  (0.05) 

6  * 

0.53  (0.36) 

6.61  (4.51) 

0.41  (0.09) 

7* 

0.63  (0.32) 

4.62  (2.87) 

0.42(0.10) 

8  * 

0.67  (0.17) 

6.39(3.18) 

0.32(0.10) 

9  * 

0.71  (0.43) 

8.99(6.51) 

0.32(0.12) 

to* 

0.27  (0.14) 

18.5(11.6) 

0.30(0.07) 

11  * 

0.72  (0.11) 

3.37(1.37) 

0.39(0.11) 

port  of  a  light  dosimetry  system,  which  consists  of  an  array  of  independent  photodiodes,  each 
connected  to  an  SMA-style  fibre  optic  connector.  The  photovoltage  generated  by  these  diodes 
is  amplified  and  recorded  by  an  analog-to-digital  (AD)  data  acquisition  board  (DataTranslation, 
Marlboro,  MA).  The  isotropic  detector  was  calibrated  in  a  15.2  cm  diameter  integrating  sphere 
using  730  nm  light.  The  measured  photovoltage  (V)  from  the  isotropic  detector  was  converted 
to  light  fluence  rate  using  the  expression, 


<h  =  A(V-B), 


(5) 


where  A  (mW  cm-2  V-1)  is  the  conversion  factor  for  fluence  rate  and  B  (V)  characterizes 
the  leakage  of  the  photodiode.  Since  the  calibration  is  performed  in  air,  when  the  isotropic 
detector  is  used  in  tissue  a  correction  factor  of  a  —  1.9  was  used.  This  value  was  measured 
using  the  response  of  the  isotropic  detector  in  and  out  of  water  medium  fpr  the  same  incident 
fluence  rate  (Vulcan  et  al  2000,  Marijnissen  and  Star  2002,  Zhu  et  al  2005a,  2005b). 

It  is  possible  to  determine  the  light  fluence  rate  per  unit  power,  <p/S,  without  independent 
calibrations  of  the  isotropic  detector  (for  <p)  and  the  light  source  power  (S).  An  integrating 
sphere  with  a  built-in  detector  port  and  two  input  ports  (one  for  the  light  source  and  the  other 
for  the  isotropic  detector)  was  used.  A  custom-made  baffle  blocked  the  direct  light  from  the 
light  source  to  either  the  built-in  detector  or  the  isotropic  detector.  The  integrating  sphere  was 
calibrated  for  the  light  fluence  rate  and  the  light  source  power  (figure  5)  as  follows.  The  fluence 
rate  in  air  4>0  is  plotted  against  the  power  reading  I  recorded  by  the  built-in  power  meter  and 
then  fitted  to  obtain  the  constant  a  (figure  5(A)),  such  that  cj>o  =  al.  The  calibration  constant  b 
is  then  obtained  from  the  fit  of  the  power  reading  from  the  integrating  sphere  as  a  function  of 
the  actual  power  (figure  5(B))  such  that  S  =  bl.  The  fluence  rate  and  the  uncalibrated  detector 
signal  ( F )  are  related  by  (po  =  gFo  in  air  and  (p  =  agF  in  tissue,  where  a  is  the  tissue  correction 
factor  and  g  is  a  constant  that  converts  the  light  fluence  rate  from  the  detector  signal  (with 
leakage  correction).  One  can  determine  the  ratio  <p/S  from  a  simultaneous  measurement  of 
the  point  source  with  the  power  S  and  the  same  isotropic  detector  in  the  integrating  sphere 
using: 


(6) 


where  Fo  is  the  isotropic  detector  reading  in  the  integrating  sphere  in  air  for  the  same  point 
source  with  power  S,  P  is  the  ratio  of  the  integrating  sphere  calibration  coefficients  a/b,  a  is 
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Figure  5.  (A)  Fluence  rate  calibration.  The  constant  a  is  determined  from  the  fit  of  the  fluence 
rate  in  air  to  that  of  the  power  meter  reading.  (B)  The  power  calibration  determines  the  constant  f>, 
from  the  fit  of  the  power  reading  from  the  integrating  sphere  to  that  of  the  power  meter. 


the  tissue  calibration  factor  and  F  is  the  detector  signal  in  the  phantom  measurement,  a  =  1 .9 
for  our  isotropic  detectors,  calibrated  in  air  and  used  in  water,  ft  =  0.172  cm-2  for  our 
integrating  sphere  (see  figure  5). 

3.4.  Fitting  algorithm 

A  Matlab-based  program1  using  a  graphical  user  interface  (GUI)  was  developed  to  analyse  the 
measured  data.  First,  the  data  are  read  into  the  program  from  the  selected  files.  Each  profile  is 
adjusted  to  account  for  the  difference  in  sensitivity  between  in-air  and  in-water  measurements 
and  divided  by  the  source  power  according  to  equation  (6).  The  program  displays  the  measured 
profile  and  its  reflection  in  the  y-axis  on  a  common  plot  along  with  the  difference  between  the 
two.  The  user  has  the  option  of  applying  an  offset  to  the  positions  recorded  by  the  motorized 

1  The  program  is  available  for  download  at  www.xrt.upenn.edu/radiation_physics/research/index.html. 
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positioner  to  place  the  peak  of  the  profile  at  x  =  0.  When  this  is  achieved,  the  profile  and  its 
reflection  match,  and  the  difference  between  them  is  minimized. 

Each  profile  is  fit  using  two  independent  optimization  algorithms.  Both  algorithms  attempt 
to  minimize  the  reduced  x2,  defined  as 


x2  = 


(N  -  Np)  (7) 

where  the  subscripts  m  and  t  denote  the  measured  and  theoretical  values,  respectively.  The 
sum  is  over  the  N  measured  data  points.  The  difference  between  N  and  Np,  the  number  of 
fitting  parameters,  gives  the  number  of  degrees  of  freedom  of  the  fit.  The  uncertainty  a  at 
each  data  point  is  estimated  based  on  the  measured  relationship  between  noise  and  signal  in 
our  photodiode  detectors,  which  is  approximately  linear  with  a  coefficient  of  0.01 5  (figure  7). 
To  determine  the  relationship  between  the  uncertainty  in  the  signal  recorded  by  our  dosimetry 
system  and  the  signal  itself,  we  have  performed  eight  in-air  measurements  of  the  same  point 
source  and  same  parallel  catheters  with  h  =  5  mm.  Figure  7  is  a  plot  of  the  standard  deviation 
in  fluence  rate  as  a  function  of  fluence  rate.  The  symbols  indicate  the  measured  standard 
deviation,  and  the  solid  line  is  a  linear  fit.  The  measured  data  for  <j>  >  600  mW  cm-2  was 
excluded  from  the  fit.  To  account  for  the  round-off  error  associated  with  conversion  from 
analog  to  digital  signal,  we  add  additional  uncertainty  to  each  point  equal  to  the  value  of 
the  least  significant  bit  of  the  12  bit  digital  signal.  In  practice,  this  additional  uncertainty  is 
significant  only  for  weak  signals. 

The  GUI  of  the  fitting  program  for  the  determination  of  optical  properties  is  shown  in 
figure  6.  The  add/delete  data  button  enables  one  to  choose/delete  the  data  one  needs  to 
plot/remove.  The  user  can  enter  the  power,  the  distance  between  the  two  catheters  h, 
calibration  ratio  fi,  the  detector  tissue  correction  factor  a,  as  well  as  the  fitting  range  for  x. 
One  has  the  choice  of  (1)  using  a  fixed  h  or  (2)  using  the  ‘optimize  h’  to  allow  the  program 
to  find  the  optimal  value  of  h.  In  the  text,  we  denoted  h  as  being  the  physical  separation 
between  the  catheters  and  italic  h  as  being  the  optimized  separation.  The  overlay  solution 
box  lets  one  overlay  a  plot  of  the  diffusion  theory  solution  for  a  given  set  of  absorption  and 
scattering  coefficients  and  separation  h.  Results  are  summarized  in  the  middle  lower  portion 
showing  the  current  file  name,  source  power  (in  mW),  distance  h  between  the  catheters,  values 
of  absorption,  scattering,  effective  attenuation  coefficients,  error  and  time  for  fitting  (in  ms). 

In  the  first  fitting,  equation  (4)  is  linearized  by  multiplying  by  r  and  taking  the  natural  log. 
A  corresponding  transformation  is  performed  on  the  measured  data,  which  is  then  subjected 
to  a  linear  fit  to  determine  /ze ff  and  /r',  from  which  /U,a  can  be  determined.  To  explicitly  take  into 
account  uncertainty  in  the  measured  data,  we  have  adopted  the  singular  value  decomposition 
algorithm  of  Press  et  al  (1992).  Because  the  linear  fitting  uses  ln(<t>r)  as  its  independent 
variable,  we  must  calculate  the  uncertainty  in  ln(<J>r),  given  by: 

3(ln(4>r))  1 

O'lnftl >r)  =  - XX - =  X0*-  W 


3$ 


The  second  fit  uses  the  differential  evolution  routine  adapted  by  Zhu  et  al  (2001)  from  the 
work  of  Stom  and  Price  (1997)  with  q,a  and  /z'  as  free  parameters.  This  nonlinear  algorithm 
minimizes  y2  directly. 

The  results  of  both  fits  are  displayed  as  they  are  calculated,  allowing  the  user  to  identify 
regions  of  the  profile  that  contribute  to  poor  fitting  and  exclude  them  from  the  fit.  After  the 
fitting  is  complete,  the  diffusion  theory  expression  is  evaluated  with  r  equal  to  0.5  cm  to 
determine  the  fluence  rate  per  unit  power,  and  the  value  of  x  2  of  each  fit  is  calculated. 

The  measurement  system  described  here  is  intended  for  use  in  the  measurement  of  intact 
human  tissues,  in  which  case  it  is  rarely  possible  to  verify  the  distance  between  individual 
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Figure  6.  The  GUI  of  the  fitting  program  for  the  determination  of  optical  properties.  The  user 
can  enter  the  power,  the  distance  h  calibration  ratio,  and  the  detector  water  correction  factor,  as 
well  as  the  fitting  range  forte.  One  has  the  choice  of  (1)  using  a  fixed  h  or  (2)  using  the  ‘optimize 
h'  for  the  program  to  find  the  optimal  value  of  h.  The  overlay  solution  box  lets  one  overlay  a 
fit  for  a  set  absorption  and  scattering  coefficient,  using  a  constant  h.  Results  summarized  in  the 
middle  lower  portion  showing  the  current  file  name,  source  power  (in  mW),  distance  h  between 
the  catheters,  values  of  absorption,  scattering,  effective  coefficients,  error  and  time  for  fitting 
(in  ms).  The  program  is  available  for  download,  (see  text  for  details). 


Figure  7.  Standard  deviation  of  fluence  rate  as  a  function  of  fluence  rate.  The  symbols  were 
measured  data  obtained  from  eight  measurements  of  the  same  point  source.  The  solid  line  is  a 
linear  fit.  The  measured  data  for  < p  >  600  mW  cm-2  were  excluded  from  the  fit. 


catheters.  To  account  to  variations  in  catheter  distance,  we  have  implemented  a  variable- 
distance  version  of  each  of  the  algorithms  described  above.  In  the  case  of  linear  fitting,  the 
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Figure  8.  Fluence  rate  per  unit  power  measured  in  one  phantom  (0.53%  Liposyn  and  0.012% 
black  ink)  at  source-detector  distances  of  3,  5  and  7  mm.  Fitting  parameters  are  listed  in 
table  3  (B). 


fitting  described  above  is  repeated  inside  a  nonlinear  fitting  algorithm  native  to  the  Matlab 
environment,  which  uses  h,  the  distance  between  catheters,  as  part  of  free  parameters  (^a, 
At',  h).  The  determination  of  optical  properties  is  still  performed  by  the  linear  algorithm. 
The  differential  evolution  algorithm  is  modified  only  slightly,  by  allowing  one  additional  free 
parameter  ( h );  In  this  case,  the  algorithm  optimizes  all  three  parameters  simultaneously. 

3.5.  Results 

Table  3  summarizes  the  optical  properties  (/xa,  /t',  /u.efr)  obtained  using  the  parallel  catheter 
measurements  and  those  obtained  from  broad-beam  measurements  in  the  same  phantoms. 
Results  are  presented  for  three  optical  phantoms  with  Liposyn  concentration  of  (A)  0.23%, 
(B)  0.53%,  (C)  1.14%  and  ink  concentrations  of  0.002%,  0.012%,  and  0.023%  for  three 
different  physical  separations  (h)  between  the  two  catheters.  The  first  column  specifies  the 
distance  between  the  light  source  catheter  and  the  detector  catheter.  The  second  column  shows 
the  optimized  distance  (/t)  between  the  catheters  that  gives  the  best  agreement  between  the 
optical  properties  determined  by  parallel-catheter  and  broad-beam  measurements.  The  values 
of  the  absorption,  scattering  and  effective  attenuation  coefficients  determined  by  parallel 
catheter  measurement  are  shown  in  columns  3-5.  These  values  are  then  compared  to  the 
values  measured  independently  in  a  broad-beam  geometry.  The  per  cent  differences  between 
the  two  are  listed  in  the  last  three  rows. 

Figure  8  shows  the  results  for  the  measurements  done  in  an  optical  phantom  (/x'  =4.19 
and  /xa  =  0.5  cm-1)  but  at  different  distances  from  the  light  source.  From  top  to  bottom,  h  = 
3,  5,  and  7  mm,  respectively. 

Figure  9  plots  the  error  in  the  optical  properties  determined  by  parallel -catheter 
measurements  for  all  optical  phantoms  used  in  the  study  for  physical  separation  of  5  mm. 
The  Liposyn  concentrations  are  0.23%,  0.53%  and  1.14%,  and  ink  concentrations  are  0.002%, 
0.012%  and  0.023%.  We  listed  fittings  using  h  of:  4,  4.5,  5  and  5.5  mm. 
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Table  3.  A  comparison  between  optical  properties  (fiz,  /r.',  Meff)  determined  by  broad-beam  and 
parallel-catheter  methods  for  three  optical  phantoms  with  Liposyn  concentrations:  (A)  0.23%, 
(B)  0.53%,  (C)  1.14%  and  ink  concentrations  (0.002%,  0.012%  and  0.023%)  for  three  different 
separations  ( h )  between  the  catheters. 

Parallel  catheters  Broad  beam  Difference  (%) 


Separation 

h 

Ma 

*4 

Mef r 

M  a 

K 

Meff 

fia 

4 

Meff 

h  =  3  mm 

2.5 

0.10 

2.04 

0.77 

(A) 

0.10 

1.73 

0.72 

2.0 

14.0 

6.9 

0.47 

2.11 

1.72 

0.48 

1.73 

1.58 

5.1 

17.9 

7.5 

0.93 

2.49 

2.64 

1.00 

1.73 

2.28 

8.6 

39.1 

15.8 

h  =  5  mm 

4.5 

0.11 

2.04 

0.82 

0.10 

1.73 

0.72 

8.0 

17.9 

13.9 

0.53 

2.21 

1.87 

0.48 

1.73 

1.58 

9.4 

27.7 

18.4 

1.10 

2.29 

2.75 

1.00 

1.73 

2.28 

9.9 

32.4 

20.6 

h  =  7  mm 

6 

0.11 

1.93 

0.79 

0.10 

1.73 

0.72 

9.0 

7.8 

9.7 

0.55 

1.91 

1.78 

0.48 

1.73 

1.58 

12.7 

6.7 

11.3 

1.18 

1.49 

2.29 

1.00 

1.73 

2.28 

15.2 

16.8 

0.4 

h  =  3  mm 

2.5 

0.10 

4.35 

1.15 

(B) 

0.10 

4.19 

1.11 

1.0 

3.8 

3.6 

0.43 

4.74 

2.47 

0.49 

4.19 

2.49 

12.9 

13.1 

0.8 

0.89 

4.30 

3.38 

0.99 

4.19 

3.53 

10.3 

2.6 

4.2 

h  =  5  mm 

4.5 

0.09 

3.92 

1.05 

0.10 

4.19 

1.11 

7.0 

6.4 

5.4 

0.56 

3.63 

2.46 

0.49 

4.19 

2.49 

13.5 

13.4 

1.2 

1.14 

3.53 

3.47 

0.99 

4.19 

3.53 

15.1 

15.8 

1.7 

h  =  7  mm 

6 

0.08 

3.35 

0.92 

0.10 

4.19 

1.11 

16.0 

20.0 

17.1 

0.59 

2.80 

2.23 

0.49 

4.19 

2.49 

20.6 

33.2 

10.4 

1.35 

2.57 

3.23 

0.99 

4.19 

3.53 

36.5 

38.7 

8.5 

h  =  3  mm 

2.5 

0.09 

6.97 

1.35 

(C) 

0.1 

9.14 

1.64 

13.0 

23.7 

17.7 

0.46 

5.05 

2.63 

0.5 

9.14 

3.7 

8.8 

44.7 

28.9 

0.81 

5.03 

3.49 

0.99 

9.14 

5.2 

18.4 

45.0 

32.9 

h  =  5  mm 

4.5 

0.1 

8.76 

1.62 

0.1 

9.14 

1.64 

0 

4.2 

1.2 

0.51 

7.40 

3.38 

0.5 

9.14 

3.7 

2.0 

19.0 

8.6 

0.92 

8.02 

4.70 

0.99 

9.14 

5.2 

7.3 

12.3 

9.6 

h  =  7  mm 

6 

0.10 

5.92 

1.32 

0.1 

9.14 

1.64 

2.0 

35.2 

19.5 

0.78 

3.97 

3.04 

0.5 

9.14 

3.7 

55.6 

56.6 

17.8 

1.85 

3.27 

4.26 

0.99 

9.14 

5.2 

86.9 

64.2 

18.1 

Figure  10  shows  the  results  of  the  optical  properties  for  three  different  scattering 
coefficients  (A)  1.73,  (B)  4.19  and  (C)  9.14  cm-1  and  ink  concentration  of  0.002,  0.012  and 
0.023%,  giving  an  absorption  coefficient  /xa  of  0.1,  0.5  and  1.0  cm-1.  The  distance  between 
the  light  source  and  the  detector  was  kept  constant,  at  5  mm.  The  scan  for  the  phantom  with 
optical  penetration  depth  of  8  =  1  / ficff  =  1.39  cm  (figure  10(A))  has  a  significant  boundary 
effect  since  the  point  source  was  placed  at  1.5  cm  from  the  boundary.  For  that  reason,  we  only 
fit  the  right  side  of  the  data  further  away  from  the  boundary.  When  the  optical  penetration 
depth  is  shorter  than  the  distance  to  the  boundary,  this  effect  becomes  negligible.  We  have  also 
plotted  in  figure  10(A),  the  value  of  4>/S  calculated  using  the  optical  properties  determined 
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Figure  9.  Error  in  the  optical  properties  determined  by  parallel-catheter  measurements  for  all 
optical  phantoms  used  in  the  study  for  physical  separation  of  5  mm.  The  Liposyn  concentrations 
are  0.23%,  0.53%  and  1.14%.  Catheter  separations  h  of  4  mm,  4.5  mm,  5  mm  and  5.5  mm  were 
plotted. 


by  the  diffusion  theory  fit,  but  calculated  using  the  higher-order  P3  approximation  (Hull  and 
Foster  2001)  (dotted  lines)  since  the  transport  albedo  ( a '  =  /x'/(/xa  +  /4.'))  can  be  significantly 
smaller  than  0.9. 

4.  Interstitial  set-up:  advanced 

4.1 .  Validation  of  the  diffusion  approximation 

For  cases  where  the  transport  albedo  of  the  turbid  medium  is  small  (a'  <  0.8),  the  diffusion 
approximation  is  known  to  fail  (Star  1997).  An  example  corresponding  to  the  worst  case  of 
albedo  ( a!  =  0.64)  in  our  study  is  shown  in  figure  1 1 ,  which  plots  the  product  of  fluence  rate 
and  radial  distance  as  a  function  of  radial  distance  for  a  Monte  Carlo  simulation  (solid  line) 
of  a  point  source  in  an  optical  phantom  with  =  1.0  cm-1  and  —  1.79  cm-1  (/xcff  = 
2.32  cm”1).  The  corresponding  solution  of  the  diffusion  equation  is  shown  by  the  dashed  line. 
Our  Monte  Carlo  algorithm  was  implemented  in  Matlab  using  the  implicit  capture  variance 
reduction  technique  described  by  Prahl  et  al  (1989)  and  implemented  in  the  commonly  used 
MCML  code  by  Wang  etal  (1995).  The  MC  solution  tends  to  deviate  more  from  the  diffusion 
solution  at  shorter  distances  ( r  Sj  0.3  cm).  The  slopes  of  the  two  solutions  are  slightly 
different  for  larger  distances.  As  a  result,  the  best  fit  to  this  data  using  diffusion  theory,  which 
is  indistinguishable  from  the  data  itself  on  the  scale  shown,  gives  an  artificially  higher  value 
for  jtia  (1.18  cm-1),  p/s  (2.14  cm”1)  and  /reff  (2.7  cm”1),  respectively.  For  comparison,  we  have 
also  plotted  the  results  of  a  Monte  Carlo  simulation  that  includes  the  effects  of  a  cylindrical 
catheter  surrounding  the  isotropic  source.  A  catheter  diameter  of  1.1  mm  and  an  index  of 
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Figure  10.  Measured  fluence  rate  per  unit  source  power  for  different  scattering  media  using  0.23, 
0.53  and  1.14%  Liposyn  concentrations  and  0.002,  0.012  and  0.023%  ink  concentrations  at  light- 
detector  distance  of  0.5  cm.  (A)  ft',  =  1.73  cm"1,  (B)  jr'  =  4.19  cm"1,  (C)  n't  =  9.14cm_1.  The 
solid  lines  are  the  fit  using  diffusion  theory.  The  optical  properties  determined  using  the  diffusion 
theory  are  shown  next  to  each  fit.  The  dotted  lines  represent  the  fluence  rate  predicted  by  the  P3 
theory  for  these  optical  properties. 


refraction  mismatch  between  the  catheter  (air)  and  tissue  of  1 .4  were  assumed.  In  this  case, 
the  presence  of  the  air  cavity  makes  less  difference  to  the  final  result  than  the  breakdown  of 
the  diffusion  approximation.  This  is  expected  because  /rcff  in  this  case  is  relatively  small. 
However,  even  in  cases  such  as  this,  where  the  diffusion  theory  is  clearly  beginning  to  fail,  the 
optical  properties  (/xa  and  yu,')  determined  by  the  two-catheter  method  are  still  within  20%  of 
the  true  values.  The  limitation  of  the  diffusion  approximation  can  be  addressed  by  the  use  of 
higher-order  approximations  (e.g.  P3  theory)  than  the  diffusion  theory.  In  figure  1 0(A),  we  have 
shown  that  the  use  of  P3  approximation  (dotted  line)  does  improve  the  agreement  between  the 
measurement  and  theory  in  most  cases.  The  apparent  disagreement  between  measurement  and 
P3  calculation  for  the  lowest  curve  results  from  the  fact  that  the  P3  calculation  used  the  optical 
properties  determined  using  the  diffusion  theory  (ixa  =  1.1  cm"1  and  fi's  =  2.29  cm"1),  which 
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Figure  11.  Effect  of  small  transport  albedo  on  diffusion  theory-based  fitting.  The  Monte  Carlo 
simulated  data  set  for  an  isotropic  emitter  in  an  infinite  medium  with  /ia  of  1.0  cm-1  and  ix's  of 
1.79  cm-1,  measured  by  a  detector  at  h  =  5  mm,  is  indicated  by  the  solid  line.  The  diffusion 
theory  solution  for  these  optical  properties  is  shown  by  the  dashed  line.  The  best  fit  of  the  diffusion 
theory  to  the  simulated  data  is  indistinguishable  from  the  data,  but  gives  /ra  of  1.18  and  fi's  of 
2.14  cm-1.  For  comparison,  the  dotted  line  indicates  the  Monte  Carlo  simulation  for  the  same 
optical  properties  with  the  source  embedded  in  a  1 . 1  mm  diameter  catheter. 

deviate  significantly  from  the  true  optical  properties  {p,s  =  1.0  cm-1  and  fi's  —  1.73  cm-1). 
As  the  albedo  of  this  phantom  was  only  0.63,  it  is  not  surprising  that  the  diffusion  theory  failed 
to  recover  its  optical  properties  accurately. 

4.2.  Monte  Carlo  simulation  of  the  air  gap  effect 

As  shown  in  figure  3(C),  the  experimental  set-up  introduces  two  air  cavity  columns,  one 
surrounding  the  detector  and  one  surrounding  the  light  source.  These  catheters  have  an  outer 
diameter  of  1.1  mm  and  an  inner  diameter  of  1.0  mm.  The  light  source  and  the  isotropic 
detector  each  have  an  outer  diameter  of  0.5  mm.  To  evaluate  the  effect  of  these  air  cavities, 
we  have  performed  Monte  Carlo  simulations.  The  optical  properties  used  were  n  =  1.4,  g  = 
0.9,  (p.a,  fis )  =  (0.10,  91.3),  (1.00, 17.9)  and  (1.01, 91.3)  cm-1,  respectively.  The  Monte  Carlo 
simulations  were  performed  in  cylindrical  coordinates.  The  simulated  volume  was  divided 
into  annular  bins  of  thickness  0.025  cm  and  height  0.025  cm.  We  have  simplified  the  catheter 
as  an  air  cavity  of  1.1  mm  diameter  in  the  tissue  phantom,  with  a  light  source  at  its  centre. 
Photons  launched  from  the  isotropic  source  were  propagated  without  absorption  or  scattering 
to  the  edge  of  the  source  catheter.  Refraction  at  the  boundary  was  accounted  for  both  for  the 
escaping  photons  and  for  any  photons  that  re-entered  the  cavity.  To  maintain  the  cylindrical 
symmetry  of  the  system,  we  have  ignored  the  air  cavity  surrounding  the  detector. 

The  air  cavity’s  main  effect  is  reducing  the  distance  light  must  travel  in  the  scattering 
medium  between  the  source  and  detector.  In  cases  with  small  jua,  the  fluence  rate  changes 


A  device  for  optical  properties 


2307 


Figure  12.  Effect  of  the  air  gap  introduced  by  the  source  catheter  on  fitting  results.  For  each  of 
three  sets  of  optical  properties  (indicated  in  the  legend),  we  plot  the  Monte  Carlo  simulated  fluence 
rate  measured  by  a  detector  at  h  =  5  mm  from  a  source  embedded  in  a  1.1  mm  diameter  catheter 
(solid  line)  and  that  in  an  infinite  medium  along  a  line  separated  by  h  =  5  mm  from  the  source 
(dotted  line). 


slowly  with  radial  distance,  so  the  air  cavity  has  little  effect.  In  contrast,  when  /ra  is  large,  the 
air  cavity  effect  is  much  more  pronounced.  In  figure  12,  we  plot  the  fluence  rates  predicted  by 
Monte  Carlo  simulations  with  (solid  line)  and  without  (dashed  line)  an  air  cavity  for  various 
sets  of  optical  properties.  In  the  cases  where  /xa  is  small  (/ua  =  0.1  cm-1),  the  air  cavity  effect 
is  negligible.  When  /xa  is  large,  however,  the  effect  becomes  more  significant,  and  it  changes 
the  shape  of  the  curve.  The  air  cavity  can  be  partially  accounted  for  by  reducing  the  value 
of  h.  For  [i cff  <  4  cm-1  and  h  =  5  mm,  the  best  average  value  of  h  is  4.5  mm  for  the  range  of 
optical  properties  studied. 

5.  Discussion  and  conclusion 

The  main  objective  of  this  study  was  to  create  a  device  that  can  assess  the  optical  properties 
(scattering  and  absorption  coefficients)  in  vivo  by  interstitial  measurements.  We  tested 
this  device  in  tissue-simulating  phantoms  with  different  optical  properties.  During  in  vivo 
measurement,  the  scanning  distance  is  typically  5  cm,  so  each  measurement  of  optical 
properties  takes  4  s.  This  is  the  time  required  to  obtain  a  useful  data  set.  Extensive 
commissioning  of  the  device  has  been  performed  to  ensure  the  accuracy  of  the  measurement  at 
this  speed.  We  compared  the  results  of  our  measurements  with  optical  properties  determined 
by  an  ex  vivo  method. 

Our  characterization  of  scattering  properties  forLiposyn  30%  yields  8.1  x  (c%IL),  which 
is  consistent  with  the  literature,  considering  the  variations  among  brands  and  batches  of  lipid 
solutions.  An  extensive  study  was  made  by  Madsen  et  al  (1992)  for  several  brands  of  India 
ink.  The  ink  concentrations  used  were  0.01-1%  and  the  total  attenuation  coefficient  was 
determined  as  a  function  of  ink  concentration  with  a  value  of  123  cm-I/%  for  Higgins  ink 
at  594  nm.  This  is  larger  than  our  value  (42.99  cm-l/%)  for  Higgins  ink  at  730  nm,  but 
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can  easily  be  attributed  to  difference  among  batches.  We  can  determine  /xa  and  /x'  of  the 
phantom  medium  to  an  uncertainty  of  better  than  5%.  This  error  is  determined  in  part  by 
comparing  extrapolated  results  from  broad  beam  measurements  against  known  ink  and  Liposyn 
concentrations.  Our  independent  method  of  determining  the  optical  properties  is  shown  in 
figure  2.  The  optical  properties  (/xa  and  /x')  were  determined  for  each  optical  phantom  under 
broad  beam  illumination.  The  black  curve  is  the  fit  to  the  data  that  excludes  the  air-phantom 
interface  and  the  background.  The  uncertainty  in  /xetf  can  be  estimated  from  the  uncertainty  in 
our  measured  data  (see  figure  7),  giving  a  relative  uncertainty  of  less  than  0.1%,  much  smaller 
than  the  uncertainty  in  the  measurement  of  the  volumes  of  the  ink  and  Liposyn  components  of 
the  phantom.  The  accuracy  of  determining  /x'  using  the  broad  beam  technique  is  within  8% 
for  the  lowest  Intralipid  concentration  of  0.23%  and  5%  for  the  0.53  and  1.14%  Liposyn.  The 
relative  error  of  determining  /x'  is  less  then  5%  with  a  maximum  of  8%. 

Comparing  the  uncertainties  of  optical  properties  obtained  with  different  source-detector 
separations  (in  table  3),  it  is  clear  that  the  best  result  is  obtained  using  h  =  5  mm.  For  the 
small  separation  (h  =  3  mm)  the  uncertainty  is  larger  because  when  the  source-to-detector 
distance  is  small,  the  effect  of  air  cavity  introduced  by  the  catheters  is  increased  due  to  the 
reduction  of  the  amount  of  scattering  material  between  the  catheters.  In  addition,  diffusion 
theory  breaks  down  when  the  detector  is  near  the  source  (r  <  lu,  where  /tr  is  the  transport 
mean  free  path  given  by  l/(/x')).  For  the  large  separation  (h  —  7  mm),  the  uncertainty  for 
phantoms  with  large  /xe tr  starts  to  increase  greatly  because  of  decrease  of  the  light  fluence  rate, 
resulting  in  larger  uncertainty  in  extrapolating  optical  properties  (see  figure  8).  Figure  8  shows 
the  fitting  results  for  an  optical  phantom  with  /xa  =  0.49  cm-1  and  /x'  =  4.19  cm-1.  Detailed 
fitting  results  are  shown  in  the  figure  as  well  as  in  table  3,  (B).  While  one  gets  reasonable 
results  of  optical  properties  at  h  =  3  mm,  the  fit  deviates  from  measured  data  near  the  source. 
The  fits  are  good  for  h  =  5  and  7  mm.  However,  some  of  the  resulting  optical  properties  at 
h  =  7  mm  deviate  from  the  true  value  by  more  than  30%  (table  3). 

There  are  several  causes  that  give  rise  to  uncertainties  in  determining  the  optical  properties 
of  a  phantom:  (1)  measurement  uncertainties  of  light  fluence  rate,  (2)  uncertainty  in  detector 
positioning  ( x )  and  distance  between  source  catheter  and  detector  catheter  (h),  (3)  air  cavity 
introduced  by  the  catheter,  and  (4)  limitations  of  the  diffusion  theory. 

The  first  source  of  error  is  the  uncertainty  of  light  fluence  rate  measurement.  We  compared 
the  data  from  eight  identical  measurements  of  the  same  point  source  and  plotted  the  standard 
deviation  as  a  function  of  fluence  rate.  The  linear  fit  (shown  in  figure  7)  gives  the  standard 
deviation  cr(0)  =  0.0150  as  a  linear  function  of  fluence  rate  0.  This  indicates  that  the  random 
error  of  our  system  corresponds  to  a  relative  uncertainty  of  1.5%,  which  is  not  the  limiting 
factor  in  our  determination  of  optical  properties. 

By  analysing  the  measured  data  using  different  distances  h,  we  determined  the  optimal 
value  of  h,  i.e.  that  which  gives  the  best  agreement  with  the  known  optical  properties  of  the 
sample.  Often,  this  optimal  h  is  smaller  than  the  measured  centre-to-centre  distance  between 
channels.  The  major  cause  of  error  in  determining  the  optical  properties  was  the  uncertainty 
in  the  distance  between  the  detector  and  the  light  source  (see  figure  9),  where  a  0.5  mm  error 
in  catheter  positioning  could  change  dramatically  the  uncertainty  of  the  fit.  In  figure  9,  it  can 
be  seen  that  the  best  fit  for  the  intended  5  mm  separation  is  given  by  h  =  4.5  mm  for  physical 
separation  of  5  mm.  Here,  we  used  the  standard  deviation  of  the  optical  properties  in  all 
the  phantoms  as  an  indication  of  the  separation  to  be  used.  As  shown  in  figure  9,  the  4  mm 
separation  gave  an  average  error  of  21%  with  a  maximum  error  of  49%,  the  5  mm  separation 
gave  an  average  error  of  17%  with  a  maximum  error  of  49%  and  the  5.5  mm  separation 
gave  an  average  error  of  44%  with  a  maximum  error  of  155%.  As  can  be  clearly  seen  from 
figure  8,  the  4.5  mm  separation  gives  the  lowest  uncertainty  with  an  average  error  of  17%  and 
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maximum  error  of  19%.  Reducing  the  physical  separation  between  the  detector  and  the  light 
source  to  4.5  mm  from  5  mm  partially  accounts  for  the  effect  of  air  cavity  introduced  by  the 
catheter.  Another  cause  of  error  was  the  positioning  of  the  detector  along  the  catheter.  This 
error  was  minimized  by  using  the  motorized  probe,  giving  a  positioning  accuracy  of  at  worst 
0.1  mm.  The  motorized  probe  is  able  to  take  approximately  800  points  per  measurement, 
further  reducing  the  error  measurement  by  increasing  the  number  of  measurement  points. 

The  range  of  validity  of  the  parallel-catheter  method  is  limited  to  cases  where  jU.eff  < 
4  cm-1.  Cases  with  /xcff  >  4  cm-1  produce  large  errors  in  /xa  and  /x'  if  h  is  different  from  the 
optimal  h  by  0.5  mm. 

The  optical  properties  (absorption  Qxa),  transport  scattering  (/x'),  and  effective  attenuation 
(/Zeff)  coefficients)  of  eleven  patients  with  locally  recurrent  prostate  cancer  were  measured 
in  situ  using  interstitial  isotropic  detectors  (see  table  3).  Measurements  were  made  at  732  nm 
before  motexafin  lutetium  (MLu)-mediated  PDT  in  four  quadrants.  /xa  and  //,'  varied  between 
0.07  and  1.62  cm-1  (mean  0.37  ±  0.24  cm-1)  and  1.1 — 44  cm-1  (mean  14  ±  11  cm-1), 
respectively.  /xa  was  proportional  to  the  concentration  of  MLu  measured  by  an  ex  vivo 
fluorescence  assay.  /xe ft  varied  between  0.91  and  6.7  cm-1  (mean  2.9  ±  0.7  cm-1), 
corresponding  to  an  optical  penetration  depth  (S  =  1  //xcff)  of  0.1-1 .1  cm  (mean  0.4  ±  0.1  cm). 
These  results  are  in  the  range  of  optical  properties  used  in  phantom  measurement. 

Pathologically,  one  can  attribute  the  difference  in  /x'  to  difference  in  cell  and  tissue 
structure  caused  by,  among  many  possibilities,  differences  in  cell  type,  the  presence  of  scar 
tissue  or  local  inflammation.  Judging  from  the  difference  between  the  mean  value  of  /x'  at 
732  nm  between  human  (Mill  cm-1)  and  dog  prostate  (3.6  i  4.8  cm-1)  (Zhu  et  al  2005a, 
2003),  our  measurement  is  sensitive  enough  to  show  the  glandular  structure  difference  between 
human  and  canine  prostates.  Since  our  measurement  was  made  in  cancerous  prostates  with 
prior  radiation  therapy,  it  is  very  possible  that  there  are  cancerous  cells  and  normal  cells, 
necrotic  cells  due  to  prior  radiation  therapy,  and  local  inflammation  due  to  the  PDT  procedure, 
all  of  which  can  contribute  to  the  heterogeneity  in  /x'.  We  find  that  in  only  14%  of  cases  was 
the  measured  /x'  larger  than  20  cm-1.  These  exceptional  values  of  /x'  may  well  indicate  the 
presence  of  abnormal  cells,  although  more  data  will  be  required  to  determine  a  correlation 
between  large  /x'  and  specific  changes  in  cellular  structure.  There  is  also  the  possibility 
of  errors  in  estimating  /x'  that  contribute  to  the  spread  of  /x'.  However,  our  current  paper 
gives  the  upper  limit  of  error  in  estimating  /x',  which  in  the  worst  case  presented  in  the 
paper,  is  about  155%  (see  figure  9  for  the  case  for  jxa  =  0.99  cm-1,  /x'  =  9.14  cm-1  and 
h  =  5.5  cm). 

Quick  and  accurate  determination  of  the  optical  properties  of  tissue  is  very  important 
in  a  variety  of  diagnostic  and  therapeutic  procedures.  We  developed  a  method  to  quickly 
determine  the  optical  properties  in  tissue  simulating  liquid  phantom  for  txa  between  0.1  and 
1.0  cm-1  and  /x'  between  1.8  and  9.0  cm-1.  Our  device  determines  the  optical  properties  /xa 
and  /x'  with  a  standard  and  maximum  deviation  of  8%  (15%)  and  18%  (32%),  respectively. 
The  high  uncertainty  in  determining  the  scattering  coefficient  comes  from  the  fact  that  the 
diffusion  theory  only  works  for  high  transport  albedo  ( a '  >  0.9).  These  errors  are  due  mainly 
to  uncertainty  of  the  distance  between  the  detector  catheter  and  light  source  catheter,  and  the 
effect  of  the  air  cavity  introduced  by  the  catheter. 
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ABSTRACT 

Among  the  challenges  to  the  clinical  implementation  of  photodynamic  therapy  (PDT)  is  the 
delivery  of  a  uniform  photodynamic  dose  to  induce  uniform  damage  to  the  target  tissue.  As  the 
photodynamic  dose  depends  on  both  the  local  sensitizer  concentration  and  the  local  fluence  rate 
of  treatment  light,  knowledge  of  both  of  these  factors  is  essential  to  the  delivery  of  uniform  dose. 

In  this  paper,  we  investigate  the  distribution  and  kinetics  of  the  photosensitizer  motexafin  lutetium 
(MLu,  Lutrin®)  as  revealed  by  its  fluorescence  emission.  Our  current  prostate  treatment  protocol 
involves  interstitial  illumination  of  the  organ  via  cylindrical  diffusing  fibers  (CDF’s)  inserted  into  the 
prostate  though  clear  catheters.  For  planning  and  treatment  purposes,  the  prostate  is  divided  into 
4  quadrants.  We  use  one  catheter  in  each  quadrant  to  place  an  optical  fiber-based  fluorescence 
probe  into  the  prostate.  This  fiber  is  terminated  in  a  beveled  tip,  allowing  it  to  deliver  and  collect 
light  perpendicular  to  the  fiber  axis.  Excitation  light  is  provided  by  a  465  nm  light  emitting  diode 
(LED)  source  coupled  to  a  dichroic  beamsplitter,  which  passes  the  collected  fluorescence 
emission  to  a  CCD  spectrograph.  Spectra  are  obtained  before  and  after  PDT  treatment  in  each 
quadrant  of  the  prostate  and  are  analyzed  via  a  linear  fitting  algorithm  to  separate  the  MLu 
fluorescence  from  the  background  fluorescence  originating  in  the  plastic  catheter.  A  computer- 
controlled  step  motor  allows  the  excitation/detection  fiber  to  be  moved  along  the  catheter,  building 
up  a  linear  profile  of  the  fluorescence  emission  spectrum  of  the  tissue  as  a  function  of  position. 
We  have  analyzed  spectral  fluorescence  profiles  obtained  in  4  patients  before  and  after  MLu- 
mediated  PDT.  We  find  significant  variation  both  within  individual  prostates  and  among  patients. 
Within  a  single  quadrant,  we  have  observed  the  fluorescence  signal  to  change  by  as  much  as  a 
factor  of  3  over  a  distance  of  2  cm.  Comparisons  of  pre-  and  post-PDT  spectra  allow  a 
quantification  treatment-induced  photobleaching.  Like  the  drug  distribution,  the  extent  of 
photobleaching  varies  widely  among  patients.  In  two  cases,  we  observed  bleaching  of 
approximately  50%  of  the  drug,  while  others  exhibited  negligible  photobleaching. 

1.  INTRODUCTION 

A  continuing  challenge  in  the  clinical  administration  of  photodynamic  therapy  (PDT)  is  the 
initiation  of  uniform  cellular  damage.  PDT  causes  cellular  damage  via  reactions  of  singlet  oxygen 
with  tissue  substrates.  The  production  of  singlet  oxygen  requires  a  sensitizing  drug,  light,  and 
available  molecular  oxygen.  Any  (or  all)  of  these  may  be  distributed  nonuniformly  in  tissue, 
resulting  in  heterogeneities  in  PDT-induced  damage.  In  an  effort  to  understand  the  causes  of 
such  heterogeneity,  we  seek  to  map  the  distribution  of  drug  in  the  prostates  of  patients 
undergoing  motexafin  lutetium  (MLu)-mediated  PDT  by  measuring  the  fluorescence  emission  of 
the  drug.  MLu  is  a  second-generation,  water-soluble  photosensitizer  with  an  absorption 
maximum  around  732  nm.1.  The  data  reported  here  was  acquired  from  patients  enrolled  in  a 
Phase  I  trial  of  MLu-mediated  PDT  for  recurrent  prostate  adenocarcinoma.  The  design  of  this 
study  was  based  on  a  prior  feasibility  study  in  a  canine  model.3 

The  measurement  of  fluorescence  emission  in  vivo  is  complicated  by  the  absorption  and 
scattering  of  light  within  the  sample  being  measured.  Variations  in  absorption  coefficient  may  be 


mistaken  for  variations  in  fluorophore  concentration.  Several  researchers  have  developed 
methods  for  reducing  the  effects  of  background  optical  properties  on  the  measured  fluorescence, 
either  through  specially  designed  optical  probes4, 5  or  by  using  independent  measurements  of 
optical  properties  to  apply  a  correction  to  the  measured  fluorescence  signal.6'11  In  the  current 
paper,  we  use  a  single  optical  fiber  as  the  source  and  detector.  This  minimizes  the  effects  of 
background  optical  properties  by  allowing  us  to  primarily  collect  light  that  has  traveled  a  short 
distance  in  the  tissue. 

2.  METHODS 


2.1  Patient  preparation  and  treatment 

The  in  vivo  results  presented  here  were  obtained  as  part  of  an  ongoing  Phase  I  clinical  trial  of 
MLu  PDT  for  the  treatment  of  recurrent  prostate  carcinoma.12  Under  this  protocol,  clear 
brachytherapy  catheters  are  placed  in  the  prostate.  Irradiation  of  the  prostate  is  accomplished  by 
placing  cylindrical  diffusing  fibers  (CDF's)  of  various  lengths  into  the  catheters.  Several  weeks 
before  the  planned  treatment,  patients,  having  given  informed  consent,  were  examined  by 
transrectal  ultrasound  (TRUS)  to  determine  the  size  and  position  of  the  prostate.  For  planning 
purposes,  the  prostate  was  divided  into  four  quadrants.  The  positions  of  the  catheters  were 
chosen  to  provide  uniform  illumination  of  the  prostate  by  maintaining  a  spacing  of  1  cm  between 
them.  The  length  of  the  CDF  in  each  catheter  was  chosen  to  cover  the  entire  length  of  the 
prostate. 

For  the  purposes  of  monitoring  the  progress  of  treatment,  one  additional  catheter  was  inserted  in 
each  quadrant.  This  catheter  held  an  isotropic  fiber  optic  based  detector  used  to  monitor  the 
local  light  fluence  during  treatment.13  This  detection  catheter  was  also  used  for  fluorescence 
measurements  before  and  after  PDT  treatment,  as  described  below. 

2.2  Fluorescence  measurement  setup 

A  schematic  of  the  fluorescence  spectroscopy  measurement  system  is  shown  in  figure  1 . 
Fluorescence  excitation  light  was  provided  by  a  465-nm  light  emitting  diode  (LED)  source  coupled 
into  an  optical  fiber.  The  beam  exiting  the  fiber  was  collimated  and  reflected  by  a  dichroic 
beamsplitter  (Chroma  Technology,  Rockingham,  VT)  with  a  cutoff  wavelength  of  600  nm.  The 
beam  was  then  refocused  into  a  fiber  optic  based  probe.  The  probe  consisted  of  a  single  optical 


Figure  1 :  Fluorescence 
spectroscopy  setup.  The 
computer  that  acquires 
and  stores  fluorescence 
spectra  also  controls  the 
position  of  the  detection 
fiber  via  a  step-motor 
positioner  (not  shown). 
The  dichroic  allows  a 
single  fiber  to  deliver 
excitation  light  and  collect 
emitted  fluorescence. 


fiber  terminated  in  a  beveled  tip.  This  fiber  emits  and  collects  light  at  right  angles  to  its  axis, 
allowing  it  to  interrogate  tissue  adjacent  to  the  catheter  in  which  it  is  placed.  Fluorescence 
collected  by  the  fiber  is  collimated  and  passed  through  the  dichroic  beamsplitter.  To  further 
discriminate  against  excitation  light,  an  OG  530  (Schott  Glass  Technologies,  Duryea,  PA)  glass 
filter  is  placed  in  the  beam  path,  blocking  light  at  wavelengths  shorter  than  530  nm.  The 
transmitted  light  is  refocused  onto  another  optical  fiber  which  terminates  at  the  focal  plane  of  a 
0.125  m  focal  length  spectrograph  (Acton  Instruments,  Acton,  MA).  The  spectrograph  images  the 
fiber  onto  a  liquid  nitrogen  cooled  CCD  (Princeton  Instruments,  Princeton,  NJ).  The  spectrograph 
and  CCD  were  configured  to  collect  spectra  with  a  pixel  width  corresponding  to  0.43  nm  over  the 
range  from  440  nm  to  940  nm.  The  actual  resolution  of  the  measurement  is  limited  by  the 
spectrograph  resolution  (~5  nm).  The  CCD  shutter  opening  and  data  acquisition  were  triggered 
by  a  transistor-transistor  logic  (TTL)  pulse  from  the  dosimetry  computer. 

Fluorescence  spectra  were  acquired  at  intervals  of  2  mm  along  the  detection  catheter  in  each 
quadrant.  A  total  length  of  40  mm  was  sampled  in  each  quadrant.  The  position  of  the 
fluorescence  probe  was  controlled  by  a  step  motor-driven  positioning  stage  (Unislide,  Velmex, 

Inc.,  East  Bloomfield,  NY)  controlled  by  the  same  computer  that  triggered  the  data  acquisition14. 
To  ensure  repeatable  positioning  of  the  probe  within  the  catheter,  the  probe  was  aligned  to  the 
closed  end  of  the  catheter  prior  to  the  initialization  of  each  scan. 

2.3  Data  analysis 

The  spectra  acquired  by  the  CCD  were  digitized  and  stored  in  binary  format.  For  each  set  of 
measurements,  an  additional  spectrum  was  acquired  with  the  shutter  closed,  to  account  for  the 
dark  charge  accumulated  by  the  CCD  and  for  the  systematic  offset  present  in  the  CCD’s  analog- 
to-digital  converter.  This  signal  was  subtracted  from  each  frame  prior  to  analysis. 

The  analysis  of  data  was  accomplished  using  a  custom-designed  graphical  user  interface  (GUI) 
written  in  the  Matlab  programming  environment.  The  user  interface  of  this  program  is  shown  in 
figure  2.  While  the  essential  data  analysis  and  programming  could  be  accomplished  in  plain  text 
programming,  the  GUI  offers  a  significant  advantage  in  that  it  allows  the  user  to  see  the  fitting 
results  as  they  are  generated,  and  to  identify  poor  fits  or  spectra  that  include  artifactual  data,  and 
to  adjust  the  fitting  range  to  compensate. 

The  spectra  were  analyzed  using  the  singular  value  decomposition  (SVD)  fitting  algorithm 
described  by  Finlay  et  al.  15  and  included  in  the  GUI.  This  algorithm  requires  the  selection  of 
basis  spectra  corresponding  to  the  known  components  of  the  fluorescence  emission  spectrum.  In 
this  case,  we  used  the  fluorescence  of  MLu,  measured  in  a  Liposyn  phantom,  as  our  primary 
basis  spectrum.  To  isolate  the  MLu  emission,  we  subtracted  the  signal  measured  in  the  same 
phantom  prior  to  the  addition  of  MLu.  The  in  vivo  chemical  environment  differs  from  that  of  our 
phantom.  As  a  result,  we  consistently  observe  a  shift  in  the  emission  maximum  of  MLu  from 
approximately  738  nm  in  the  phantom  to  745  nm  in  vivo.  We  have  accounted  for  the  wavelength 
shift  by  digitally  shifting  our  basis  spectrum  prior  to  fitting. 

The  second  component  of  fluorescence  we  observe  in  vivo  and  in  phantoms  arises  from 
fluorescence  in  the  optical  fiber  and  in  the  plastic  catheter.  While  this  background  fluorescence 
cannot  be  eliminated  given  the  constraints  of  our  clinical  protocol,  it  can  be  characterized  by 
measurements  in  fluorophore-free  media,  providing  a  second  basis  spectrum 

The  SVD  algorithm  we  employ  also  includes  a  61-term  Fourier  series15  to  account  for 
fluorescence  of  unknown  origin.  The  Fourier  components  are  given  much  lower  weight  in  the 
fitting  routine  than  the  basis  spectra  of  known  fluorophores  to  restrict  their  application  to 
components  of  the  spectrum  that  cannot  be  fit  by  combinations  of  these  species.  In  the  cases 
presented  here,  the  Fourier  components  constitute  only  a  minor  contribution  to  the  total  fit, 
indicating  that  the  known  fluorophores  adequately  account  for  the  fluorescence  we  observe. 
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Figure  2:  Graphical  user  interface  (GUI)  used  to  analyze  fluorescence  emission  spectra.  The  GUI  provides 
a  front  end  for  the  SVD  algorithm  and  gives  real-time  feedback,  allowing  interactive  fitting. 

The  final  result  of  the  SVD  fitting  is  an  amplitude  for  each  of  the  known  basis  spectra,  namely 
MLu  and  background.  Because  the  background  signal  originates  in  the  measurement  apparatus, 
its  intensity  is  unchanged  from  one  treatment  to  the  next.  We  can  therefore  use  the  intensity  of 
the  background  to  correct  for  variations  in  lamp  output,  fiber  coupling  efficiency,  and  beamsplitter 
alignment.  The  MLu  amplitude  we  report  has  been  divided  by  the  background  amplitude  to 
account  for  these  variations. 

The  amplitude  of  the  MLu  fluorescence  obtained  by  this  method  in  vivo  is  approximately  10  times 
lower  than  that  obtained  by  absorption  spectroscopy  (see  section  3.2).  This  is  the  result  of  the 
fact  that  our  basis  spectrum  was  measured  in  a  phantom  where  the  only  absorption  was  that  due 
to  MLu  itself.  In  the  in  vivo  environment,  other  absorbers  are  present,  which  leads  to  a  reduction 
in  the  measured  intensity.  It  is  possible  that  variations  in  absorption  among  patients  or  among 
different  positions  within  the  same  prostate  could  cause  deviations  from  this  ideal  correction 
factor. 

One  can  measure  the  MLu  concentration  by  an  alternate  method  such  as  absorption 
spectroscopy  (see  section  3.3),  and  use  this  value  to  determine  a  correction  factor  for  each 
prostate  or  even  for  each  quadrant.  However,  in  practice,  the  availability  of  reliable  absorption 
data  is  often  limited  by  localized  bleeding  and  by  time  constraints.  Of  the  four  patients  studied 
thus  far,  we  have  reliable  absorption  and  fluorescence  data  in  multiple  quadrants  both  before  and 
after  treatment  for  only  one  patient,  number  13.  We  have  adopted  a  single  correction  factor  equal 
to  10  obtained  by  comparing  the  fluorescence  amplitude  of  MLu  measured  in  vivo  in  prostate 
number  13  with  absorption  measurements  at  the  same  position  in  the  same  prostate  (see  section 
3.3  below).  This  conversion  yields  a  concentration  of  12.9  mg/kg  for  a  spectrum  whose  MLu 
concentration  has  a  peak  intensity  of  1 .0  after  normalization.  This  factor  is  in  approximate 
agreement  with  the  results  of  phantom  experiments  in  which  ink  was  added  to  the  phantom  to 
simulate  absorption  by  biological  chromophores. 


2.4  Monte  Carlo  Simulation  and  Analytic  Modeling 


To  asses  the  effects  of  optical  property  variations  on  measured  fluorescence,  we  have  performed 
Monte  Carlo  simulations  in  which  the  source  was  modeled  as  a  pencil  beam  and  the  detector  as 
a  circular  aperture  with  a  radius  of  500  microns  and  an  acceptance  numerical  aperture  of  0.22. 
The  relative  changes  in  fluorescence  resulting  from  changes  in  the  absorption  coefficients  at  the 
excitation  (pax)  and  emission  (pam)  wavelengths  and  in  a  simultaneous  change  in  the  scattering 
coefficient  (p/)  are  summarized  in  table  1 .  In  each  case,  the  change  shown  is  the  average  over 
several  simulations  in  which  one  variable  was  changed  and  the  others  were  held  constant. 


These  simulations  indicate  that  changes  in 
absorption  at  the  emission  wavelength  will  have 
negligible  effect  on  the  measured  fluorescence, 
since  pa  at  the  emission  wavelength  (745  nm)  is 
generally  much  smaller  than  that  at  the  excitation 
wavelength  (465  nm).  Note  that  a  tenfold  increase 
in  p.am  caused  less  than  10%  change  in 
fluorescence.  Even  at  the  excitation  wavelength,  a 
doubling  of  pa  caused  less  than  a  20%  change  in 
measured  fluorescence.  Similarly,  a  doubling  of 
the  scattering  coefficient  at  both  wavelengths 
caused  only  a  22%  change  in  fluorescence.  Correction  for  individual  variations  in  optical 
properties  is  therefore  unlikely  to  be  necessary  unless  the  variations  in  pa  are  very  great. 


Parameter 

Change 

Change  in 
fluorescence 

Pax  (cm1) 

2  ->  4 

-19% 

Mam  (cm1) 

0.05  0.5 

-8.1% 

Psx’  (cm"1) 

10 20 

+  22  % 

Psx’  (cm1) 

5  ->  10 

Table  1:  Effect  of  changing  optical 
properties  on  the  measured  fluorescence 
signal  in  a  pencil-beam  geometry  in  an 
infinite  medium. 


Based  on  these  results,  we  designed  a  set  of  experiments  in  tissue-simulating  phantoms 
containing  Liposyn  as  a  scatterer  and  MLu.  In  the  first,  the  pa  of  the  phantom  was  varied  by 
varying  the  MLu  concentration  (pam  =  0-0.6  cm'1,  pax  =  0-1.6  cm"1),  while  p/  was  held  constant 
(psm'  =  4  cm"1,  psx'  =  6  cm"1).  In  the  second,  the  pa  was  varied  by  adding  black  ink  (pam  =  0-0.4 
cm"1,  pax  =  0-3  cm"1  psm'  =  6  cm"1,  psx  =  10  cm"1).  In  each  case,  the  fluorescence  was  analyzed  as 
described  above  and  the  resulting  MLu  signal  was  normalized  to  the  known  MLu  concentration. 
The  MLu  fluorescence  measured  in  these  phantoms  is  shown  in  figure  3(a)  as  function  of 
excitation-wavelength  pa.  The  data  points  from  the  two  phantoms  do  not  overlap  because  the 
values  of  ps'  and  pam  for  the  two  phantoms  were  different.  The  change  in  fluorescence  induced 
by  changing  pa  appears  to  be  larger  than  that  predicted  by  Monte  Carlo  simulation. 


As  an  alternative  to  Monte  Carlo  simulation,  the  effects  of  absorption  and  scattering  on  measured 
fluorescence  in  semi-infinite  media  can  be  modeled  using  the  forward-adjoint  fluorescence 
Scheme  proposed  by  Crilly  et  a/.16.  Briefly,  this  method  models  the  forward  propagation  of 
excitation  light  from  the  source  and  the  time-reversed,  or  adjoint,  propagation  of  positional 
importance  from  the  detector.  The  positional  importance  is  defined  as  the  probability  that  a 
photon  emitted  at  a  point  is  eventually  captured  by  a  detector.  The  volume  integral  of  the  product 
of  the  excitation  fluence  rate  and  the  importance  is  proportional  to  the  measured  signal.  Finlay 
and  Foster9  have  derived  an  analytic  solution  to  this  model  for  the  case  of  an  isotropic  point 
source  and  an  isotropic  detector  in  an  infinite  homogeneous  medium,  which  for  the  diffusion 
approximation  takes  the  form: 
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where  rscj  is  the  distance  between  the  source  and  detector,  D  is  the  diffusion  constant,  //eff  is  the 

effective  attenuation  coefficient,  the  subscripts  x  and  m  denote  the  excitation  and  emission 
wavelengths,  respectively,  and  F0  is  the  intrinsic  fluorescence  of  the  sample.  This  formula 


Figure  3:  (a)  MLu  fluorescence  measured  in  two  sets  of  phantoms,  plotted  as  a  function  of  excitation- 
wavelength  absorption  coefficient.  See  text  for  details.  The  solid  lines  indicate  the  values  predicted  by  eq.  1 , 
and  (b)  the  corresponding  fluorescence  corrected  for  optical  effects  (Fo)  found  from  eq.  1. 

assumes  that  both  the  source  and  detector  are  points.  In  the  measurements  reported  here,  a 
single  fiber  serves  as  both  source  and  detector;  however  it  is  finite  in  extent  and  non-isotropic. 

To  approximate  the  true  physical  situation,  we  assume  that  the  light  beam  exiting  or  entering  our 
probe  is  a  pencil  beam,  which  can  be  represented  by  a  point  source  at  a  distance  of  1  transport 
mean  free  path  in  the  z  direction.17  To  account  for  the  finite  extent  of  the  probe,  we  introduce  a 
lateral  shift  in  the  x  direction  between  the  source  and  detector  pencil  beams.  We  have  applied 
equation  1  to  the  case  of  the  phantoms  described  above.  The  results  are  shown  by  the  solid 
lines  in  figure  3(a).  In  this  case  we  used  a  lateral  shift  of  0.45  cm.  It  should  be  emphasized  that 
this  lateral  shift  was  determined  empirically,  and  does  not  necessarily  correspond  to  any  physical 
dimension  of  the  probe. 

To  correct  for  absorption  and  scattering  effects,  we  solve  equation  1  for  F0.  The  data  originally 
shown  in  figure  3(a),  corrected  by  this  method,  are  shown  in  figure  3(b).  This  correction  reduces 
the  nearly  10-fold  variation  seen  in  figure  3(a)  to  slightly  more  than  a  factor  of  2.  This  method  can 
be  used  in  the  phantom  case  only  because  the  optical  properties  of  the  phantom  were  known  a 
priori.  To  implement  such  a  correction  in  vivo  requires  an  accurate  measurement  of  the  local  in 
vivo  optical  properties.  This  is  difficult  in  a  clinical  situation  because  (a)  it  is  often  the  case  that 
localized  bleeding  prevents  accurate  measurement  of  absorption  and  (b)  when  a  measurement  of 
(ia  is  available,  the  uncertainty  in  pa  may  be  as  large  as  the  correction  being  made,  rendering  the 
uncertainty  in  the  corrected  measurement  too  large  to  make  the  correction  meaningful. 
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Figure  4:  SVD  analysis  of  a  typical 
fluorescence  emission  spectrum.  The 
spectrum  is  separated  into  MLu  and 
background  components  and  a  small 
residual  composed  of  a  Fourier  series. 


3.  RESULTS 


3.1  MLu  concentration 

We  have  measured  fluorescence  spectra  for  4  patients  thus  far.  The  results  of  spectral  analysis 
of  a  typical  fluorescence  spectrum  are  shown  in  figure  4.  The  noisy  line  indicates  the  original 
data,  and  the  dotted  and  dashed  lines,  the  components  determined  by  the  SVD  algorithm.  The 
sum  of  these  is  represented  by  the  solid  line  that  closely  matches  the  data.  The  sum  of  the 
contributions  of  the  terms  of  the  Fourier  series,  labeled  ‘residual’  on  the  plot,  is  much  smaller  in 
amplitude  than  the  contributions  of  the  background  fluorescence  and  MLu,  indicating  that  these 
known  components  accurately  model  the  majority  of  the  measured  fluorescence 

3.2  Spatial  distribution  of  MLu 

Figure  5  illustrates  the  distributions  of  MLu  within  three  of  the  prostates  we  have  measured. 
Profiles  measured  immediately  before  (dashed  lines)  and  immediately  after  (solid  lines)  PDT 
treatment  are  similar  in  shape  and  amplitude.  Only  the  profiles  from  patient  12  indicate  significant 
reduction  in  MLu  concentration. 

In  one  patient  (13),  we  have  obtained  fluorescence  profiles  in  all  four  quadrants  of  the  prostate 
both  before  and  after  PDT  treatment.  These  profiles  are  shown  in  figure  6.  As  in  the  cases 
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Figure  5:  MLu  fluorescence  profile  as  a  function  of  position  within  the  prostate  for  patients  12  (a),  15 
(b  and  c)  and  17  (d)  The  profiles  measured  before  (dashed  line)  and  after  (solid  line)  PDT  treatment 
are  similar  in  shape  and  amplitude.  The  patient  number  and  quadrant  from  which  each  data  set  was 
taken  are  indicated  in  the  title  of  each  frame. 
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Figure  6:  MLu  fluorescence  profiles  acquired  in  the  four  quadrants  of  a  single  prostate  before 
(dashed  lines)  and  after  (solid  lines)  PDT  treatment.  The  quadrant  from  which  the  data  was  taken  is 

indicated  in  the  title  of  each  frame. 

shown  in  figure  4,  the  profiles  taken  before  and  after  treatment  show  similar  shape  and 
concentration,  with  the  exception  of  the  RUQ  (panel  b),  which  shows  a  similar  shape  but  reduced 
concentration  compared  to  the  pre-treatment  profile.  The  measurements  in  the  upper  quadrants 
(LUQ  and  RUQ)  were  made  1 .5  cm  anterior  to  those  in  the  lower  quadrants  (LLQ  and  RLQ, 
respectively),  while  those  in  the  right  quadrants  (RUQ  and  RLQ)  were  shifted  2  cm  laterally  from 
those  in  the  left  quadrants  (LUQ  and  LLQ,  respectively). 

Among  the  four  prostates  we  have  measured,  we  consistently  observe  variations  of  a  factor  of  3 
in  fluorescence  signal  over  distances  as  short  as  a  few  millimeters.  In  the  majority  of  cases,  we 
find  that  the  spatial  distribution  of  sensitizer  measured  prior  to  PDT  treatment  is  similar  to  that 
acquired  after  treatment.  The  variation  in  MLu  concentration  within  individual  prostates  is  as  great 
as  the  variation  among  patients,  as  indicated  in  Table  2.  Each  of  the  four  patients  listed  were 
given  MLu  at  a  concentration  of  2  mg/kg  body  weight  3  hours  prior  to  the  first  set  of 
measurements. 


The  maximum  fluorescence  observed  in  each  patient  varied  from  5  to  13  mg/kg,  significantly 
higher  than  the  injected  concentration.  This  is  consistent  with  the  results  of  ex  vivo 
measurements  made  on  biopsies  of  the  same  tissues  (data  not  shown).13  In  one  case  (patient 
#12)  we  observed  significant  photobleaching  of  the  sensitizer  during  treatment.  In  another 
(patient  #  13),  we  observe  photobleaching  by  approximately  a  factor  of  two  in  one  quadrant,  but 
no  significant  photobleaching  in  other  quadrants. 


Patient 

MLu  Concentration 

After  (mg/kg) 

#12 

2.5 -5.0 

0.4-1. 7 

#13 

2-10.5 

1-11 

#15 

1-16 

2-13 

#17 

3.4-9. 4 

0.7-12.4 

Table  2:  Range  in  MLu  concentration  measured  by  fluorescence 
spectroscopy  before  and  after  PDT  treatment.  Patient  #13  exhibited 
photobleaching  in  one  quadrant  only. 


3.3  Comparison  with  absorption  spectra 

For  one  patient,  we  have  measured  profiles  of  MLu  concentration  using  not  only  fluorescence 
spectra  but  also  diffuse  transmission  spectra.  These  were  acquired  by  placing  a  white  light 
source  in  a  catheter  parallel  to  the  detection  catheter  and  measuring  the  fluence  rate  spectrum  in 
the  detection  catheter  with  an  isotropic  detector,  as  reported  previously.14, 18  A  series  of  such 
spectra  can  be  fit  using  photon  diffusion  theory  to  determine  the  absorption  and  scattering 
spectra  of  the  prostate  as  a  function  of  position  along  the  catheter.  The  resulting  absorption 
spectra  can  then  be  fit  using  the  SVD  algorithm  described  above  to  determine  the  concentrations 
of  oxy-  and  deoxy-hemoglobin  and  MLu.  In  figure  7,  we  plot  the  concentrations  of  MLu 
determined  by  absorption  and  fluorescence  spectroscopy.  The  comparison  between  the  MLu 
fluorescence  amplitude  and  the  MLu  concentration  determined  by  absorption  measurement  in 
this  prostate  was  the  basis  for  the  conversion  between  MLu  amplitude  and  absolute  concentration 
described  in  section  2.3.  It  is  therefore  expected  that  the  two  measurements  agree  in  absolute 
concentration. 

The  agreement  in  the  spatial  distributions  measured  by  absorption  and  fluorescence  is  striking, 
and  confirms  that  the  fluorescence  and  absorption  spectroscopy  measurements  are  in  fact 
measuring  the  same  distribution.  In  the  case  shown  here,  the  shape  of  the  MLu  distribution 
remains  similar  before  and  after  PDT,  however  the  MLu  concentration  as  reported  by  both 
methods  decreases  by  nearly  a  factor  of  two,  especially  in  the  region  around  2  to  3  cm,  in  the 
center  of  the  prostate. 


RUQ  RUQ 


Figure  7:  MLu  fluorescence  (open  circles)  and  absorption  (filled  squares)  profiles  as  functions  of  position 
within  the  prostate  of  patient  #13.  The  profiles  measured  before  (panel  a)  and  after  (panel  b)  PDT  treatment 
are  similar  in  shape,  but  indicate  some  photobleaching,  especially  in  the  center  of  the  prostate. 


4.  DISCUSSION 


We  have  demonstrated  the  ability  of  fluorescence  spectroscopy,  coupled  with  a  linear  fitting 
algorithm,  to  quantitatively  measure  the  concentration  and  distribution  of  MLu  in  the  human 
prostate.  There  may  be  significant  variations  in  optical  properties  among  prostates13  and  within  a 
single  prostate14.  The  theoretical  work  presented  in  section  2.4  indicates  that  if  the  reliability  and 
accuracy  of  absorption  measurements  can  be  improved,  these  effects  can  be  corrected  for  each 
individual  prostate.  Given  the  constraints  of  our  current  ability  to  make  absorption  measurements, 
however,  we  have  adopted  a  single,  uniform  correction  factor.  The  fact  that  the  MLu 
concentration  profiles  measured  by  two  independent  spectroscopic  techniques  agree  (see  figure 
7)  indicates  that  this  variation  will  not  significantly  hinder  fluorescence  spectroscopy. 

While  the  absorption  and  fluorescence  measurements  reported  here  agree  in  the  MLu 
concentration,  the  two  measurements  are  not  interchangeable.  First,  absorption  spectroscopy 
measurements  require  multiple  spectra  at  each  point  of  interest,  and  necessarily  average  over  a 
large  volume  of  tissue  to  make  a  measurement.  4  This  limits  the  resolution  of  absorption 
spectroscopy,  and  makes  it  more  time-consuming  than  fluorescence  spectroscopy.  While 
fluorescence  spectroscopy  lacks  the  sensitivity  of  absorption  measurements  to  scattering  and 
hemoglobin  dynamics,  it  offers  significant  advantages.  First,  because  each  measurement 
requires  only  the  acquisition  of  a  single  spectrum,  fluorescence  measurements  can  be  made  very 
quickly,  making  it  practical  to  increase  the  number  of  measurement  points  in  each  scan.  This  is 
reflected  in  figure  6,  where  the  fluorescence  measurement  points  are  more  densely  spaced  and 
cover  a  larger  range  than  the  absorption  measurement  points.  Second,  the  fluorescence 
measurement  averages  a  small  volume  of  tissue  due  to  the  limited  penetration  of  465  nm  light. 
This  allows  a  finer  resolution  than  can  be  achieved  using  absorption  measurement. 

Finally  and  perhaps  most  importantly  from  a  clinical  point  of  view,  the  fluorescence  measurement 
requires  only  a  single  catheter  to  acquire  a  profile,  while  the  absorption  measurements  require 
two.  In  our  clinical  experience,  we  have  often  encountered  cases  where  the  pooling  of  blood 
around  one  or  more  catheters  renders  the  absorption  spectroscopy  data  uninterpretable.  In  many 
of  these  cases,  fluorescence  spectra  could  still  be  measured  and  analyzed. 

The  goal  of  the  measurements  described  here  is  to  map  the  MLu  distribution  in  the  human 
prostate.  The  primary  motivation  for  this  effort  lies  in  the  fact  that  the  effective  application  of 
photodynamic  therapy  depends  on  the  creation  of  reactive  singlet  oxygen.  This  requires 
sufficient  tissue  oxygenation,  sufficient  intensity  of  treatment  light,  and  sufficient  concentration  of 
the  sensitizing  drug.  Our  efforts  at  optimizing  PDT  delivery  have  focused  on  the  latter  two  factors. 
The  goal  of  light  and  drug  optimization  is  to  deliver  an  optimal  photodynamic  dose  (defined  as  the 
product  of  light  fluence  rate  and  drug  concentration)  uniformly  over  the  target  volume.  As  the 
clinician  has  only  limited  control  over  the  drug  concentration  in  vivo  this  optimization  must  be 
accomplished  by  adjusting  the  delivered  light  intensity.  The  optimal  light  intensity  distribution  can 
be  determined  only  if  the  drug  distribution  can  be  measured. 

In  theory,  a  generic  drug  distribution  could  be  developed  and  used  for  all  patients.  However,  the 
variation  among  prostates  makes  this  impractical.  In  addition,  we  observe  variations  in  MLu 
fluorescence  within  individual  prostates  that  is  equal  to  or  greater  than  the  variation  among 
different  patients’  prostates.  This  indicates  that  characterization  of  the  tissue  MLu  concentration 
for  each  patient  is  necessary,  but  will  only  result  in  optimal  treatment  of  the  point  of  measurement. 
To  optimize  treatment  for  the  entire  prostate,  it  will  be  necessary  to  build  up  a  map  of  the  MLu 
concentration  in  three  dimensions.  This,  coupled  with  a  similar  map  of  tissue  optical  properties, 
will  allow  the  optimization  of  the  light  source  distribution  to  provide  uniform  photodynamic  dose. 

As  a  first  step  to  creating  a  map  of  MLu  distribution,  we  have  acquired  fluorescence  profiles  in  all 
four  quadrants  of  the  prostate  of  one  patient,  shown  in  figure  6.  As  expected,  we  see  significant 
variation  among  the  four  quadrants.  However,  the  general  features,  namely  a  peak  in  MLu 


concentration  around  5  to  15  mm  and  a  shallow  minimum  around  30  mm,  are  reproduced  in  three 
out  of  the  four  quadrants.  This  indicates  that  the  variations  in  MLu  concentration  occur  on  a  scale 
of  approximately  1.5  to  2  cm  (the  spacing  of  the  catheters  used  to  make  these  measurements)  in 
all  three  dimensions.  We  therefore  expect  a  set  of  measurements  with  a  catheter  spacing  of  1 
cm  to  be  sufficient  to  characterize  the  MLu  distribution  in  a  typical  prostate.  This  can  be  achieved 
in  our  current  clinical  protocol  by  using  the  catheters  currently  reserved  for  treatment,  in  addition 
to  the  four  dedicated  detector  catheters,  for  fluorescence  measurement.  Work  in  this  area  is 
ongoing. 

The  observation  of  MLu  photobleaching  may  have  significant  implications  for  photodynamic 
dosimetry.  It  has  long  been  appreciated  that  photobleaching  may  reduce  the  photodynamic  dose 
by  reducing  the  sensitizer  concentration  as  time  progresses.19  More  recently,  it  has  been 
suggested  that  photobleaching  itself  can  be  used  to  monitor  the  progress  of  photodynamic 
therapy.20  For  the  case  of  singlet-oxygen-mediated  photobleaching,  photodynamic  destruction  of 
target  tissues  and  photobleaching  of  the  drug  are  inseparably  related.21  If  this  is  the  case  for  MLu, 
it  is  possible  that  the  quadrants  where  we  observe  significant  photobleaching  are  those  where  the 
greatest  dose  of  singlet  oxygen  is  delivered.  An  alternative  explanation  is  that  those  quadrants 
where  no  photobleaching  or  a  slight  increase  in  MLu  concentration  is  observed  are  better 
vascularized,  and  have  the  bleached  sensitizer  replaced  from  the  circulating  pool  or  from 
adjacent  tissue.  These  two  explanations  make  opposite  predictions  concerning  the  effectiveness 
of  treatment  in  those  quadrants  that  exhibit  photobleaching  versus  those  that  do  not.  Further 
research  into  the  photobleaching  behavior  and  in  vivo  vascular  effects  of  MLu  is  needed  to 
resolve  the  underlying  dynamics  responsible  for  the  changes  in  concentration,  or  lack  thereof, 
that  we  observe. 
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APPENDIX  5 


Optimization  of  light  sources  for  prostate  photodynamic  therapy 
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Department  of  Radiation  Oncology,  University  of  Pennsylvania,  Philadelphia,  PA 

ABSTRACT 

To  deliver  uniform  photodynamic  dose  to  the  prostate  gland,  it  is  necessary  to  develop  algorithms  that  optimize 
the  location  and  strength  (emitted  power  x  illumination  time)  of  each  light  source.  Since  tissue  optical 
properties  may  change  with  time,  rapid  (almost  real-time)  optimization  is  desirable.  We  use  the  Cimmino 
algorithm  because  it  is  fast,  linear,  and  always  converges  reliably.  A  phase  I  motexafin  lutetium  (MLu)-mediated 
photodynamic  therapy  (PDT)  protocol  is  on-going  at  the  University  of  Pennsylvania.  The  standard  plan  for  the 
protocol  uses  equal  source  strength  and  equal  spaced  loading  (1-cm).  PDT  for  the  prostate  is  performed  with 
cylindrical  diffusing  fibers  (CDF)  of  various  lengths  inserted  to  longitudinal  coverage  within  the  matrix  of 
parallel  catheters  perpendicular  to  a  base  plate.  We  developed  several  search  procedures  to  aid  the  user  in 
choosing  the  positions,  lengths,  and  intensities  of  the  CDFs.  The  Cimmino  algorithm  is  used  in  these  procedures 
to  optimize  the  strengths  of  the  light  catheters  at  each  step  of  the  iterative  selection  process.  Maximum  and 
minimum  bounds  on  allowed  doses  to  points  in  four  volumes  (prostate,  urethra,  rectum,  and  background) 
constrain  the  solutions  for  the  strengths  of  the  linear  light  sources.  Uniform  optical  properties  are  assumed.  To 
study  how  different  opacities  of  the  prostate  would  affect  optimization,  optical  kernels  of  different  light 
penetration  were  used.  Another  goal  is  to  see  whether  the  urethra  and  rectum  can  be  spared,  with  minimal 
effect  on  PTV  treatment  delivery,  by  manipulating  light  illumination  times  of  the  sources.  Importance  weights 
are  chosen  beforehand  for  organ  volumes,  and  normalized.  Compared  with  the  standard  plan,  our  algorithm  is 
shown  to  produce  a  plan  that  better  spares  the  urethra  and  rectum  and  is  very  fast.  Thus  the  combined 
selection  of  positions,  lengths,  and  strengths  of  interstitial  light  sources  improves  outcome. 

Keywords:  photodynamic  therapy,  prostate,  Cimmino  Optimization,  light  dosimetry,  tissue  optical  properties, 
diffusion  theory. 


I.  INTRODUCTION 

Photodynamic  therapy  (PDT)  is  a  treatment  modality  employing  light  of  an  appropriate  wavelength  in  the  presence 
of  oxygen  to  activate  a  photosensitizing  drug  which  then  causes  localized  cell  death  or  tissue  necrosis.  Using  a  surface 
illumination  technique,  PDT  has  been  used  to  treat  many  superficial  tumors  including  skin,  lung,  esophagus,  and 
bladder  (1).  This  technique  is,  however,  inadequate  when  applied  to  large  bulky  tumors  or  solid  organs  due  to  limited 
light  penetration  into  tissue.  A  more  efficient  illumination  scheme  would  be  interstitial  light  delivery  whereby  optical 
fibers  are  placed  directly  into  the  bulky  tumors  or  organs. 

The  prostate  gland  is  an  organ  that  appears  to  be  a  good  target  for  interstitial  PDT.  Tumors  of  the  prostate  are  often 
confined  to  the  prostate  itself  and  brachytherapy  techniques  used  for  the  placement  of  radioactive  seed  implants  can  be 
adapted  for  the  placement  of  interstitial  optical  fibers  (2).  Several  preclinical  studies  have  evaluated  the  feasibility  of 
delivering  PDT  to  the  prostate  via  this  interstitial  approach  (3-7).  The  development  of  this  light  delivery  technique  has 
necessitated  an  improved  understanding  of  light  dosimetry,  critical  in  planning  the  configuration  of  multiple  fibers 
within  the  organ  or  tumor.  Based  on  results  of  a  preclinical  study  in  canine  (8),  we  have  initiated  a  motexafin  lutetium 
(MLu)-mediated  PDT  of  the  prostate  in  human  at  University  of  Pennsylvania^*))  MLu  is  a  second  generation  synthetic 
photoactive  drug  that  has  a  Q-band  absorption  peak  at  732  nm.  (10-11)  Ideal  optimization  of  the  photodynamic  linear 
light  sources  depends  on  knowledge  of  the  spatial  distributions  of  (1)  tissue  light  opacity  within  the  prostate,  (2) 
photosensitizing  drug,  and  (3)  oxygen  within  tissue.  Since  these  spatial  distributions  can  vary  in  time,  measurements 
must  be  done  just  prior  to  the  clinical  procedure.  Moreover,  the  opacity  distribution  may  be  affected  by  bleeding 
associated  with  insertion  of  the  light  sources  and  ideally  should  be  monitored  at  a  significant  number  of  points  during 
the  entire  procedure.  At  present,  measurements  at  more  than  a  few  sample  points  within  the  patient  during  the  clinical 
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procedure  are  difficult  to  make.  The  state  of  the  art  is  to  obtain  measurements  before  the  actual  clinical  procedure,  and 
to  assume  that  during  the  procedure  all  these  distributions  are  uniform  (of  constant  value)  in  the  prostate  region  and 
static  in  time.  A  comprehensive  study  of  the  inhomogeneous  light-opacity  distribution  in  vivo  is  performed  and 
reported  elsewhere.  (12-13)  For  the  present  study,  we  will  concentrate  on  optimization  of  light  fluence  distribution. 

A  number  of  optimization  algorithms  used  in  brachytherapy  are  of  interest  for  prostate  photodynamic  therapy.  The 
most  common  ones  are  simulated  annealing  algorithms  (14-16)  and  genetic  algorithms  (17-19).  Gradient  algorithms 
also  have  been  applied  (20).  In  general,  gradient  algorithms  give  reproducible  solutions  but  may  be  trapped  in  local 
minima  far  from  the  global  minimum  (21).  Simulated  annealing  and  genetic  algorithms  avoid  getting  trapped  in  local 
minima,  but  are  relatively  slow  because  they  are  stochastic  algorithms. 

To  the  best  of  our  knowledge,  optimization  algorithms  for  photodynamic  therapy  are  not  yet  in  the  literature.  In 
this  study,  we  describe  a  systematic  search  procedure  based  on  the  Cimmino  feasibility  algorithm  (22,23)  to  obtain  the 
locations  and  strengths  of  light  sources  for  photodynamic  treatment.  The  Cimmino  algorithm  is  an  iterative  linear 
algorithm  which  was  first  applied  to  radiotherapy  inverse  problems  by  Censor  et  al.  (24-26).  The  algorithm  is  safer  than 
most  common  optimization  algorithms  outlined  above  since  it  always  converges  and,  if  no  solution  exists  for  the 
inequalities  (i.e.  the  prescribed  dose  constraints  are  not  all  satisfied),  the  Cimmino  algorithm  reverts  to  a  least-square 
solution  (26). 

II.  METHODS  AND  MATERIALS 
1.  Diffusion  theory  and  Determination  of  optical  properties 


The  transport  scattering  (p’s)  and  absorption  (pa)  coefficients  characterize  the  scattering  and  absorption  properties 
of  tissue.  With  the  diffusion  approximation,  the  light  fluence  rate  tf>  at  a  distance  r  from  a  point  source  of  power,  S,  can 
be  expressed  as  (27) 
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where  S  is  the  power  of  the  point  source  (in  mW);  </>(r)  is  the  the  fluence  rate  in  mW/cm2;  juejff  =  ^3  •  jua  ■  ju's  (27)  is 
the  effective  attenuation  coefficient  in  tissues  and  is  applicable  for  a  wider  range  of  pa  and  ps’  combinations  than  the 
traditional  definition  of  =  V3  /V0“'.s+/G)  (28). 

For  a  cylindrical  diffusing  fiber  (CDF)  with  length  /,  the  light  fluence  rate  can  be  calculated  by  a  superposition  of 
Eq.(l): 
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-where  s  is  the  energy  released  by  light  per  unit  time  per  unit  length,  also  called  unit  length  source  strength  (mW/cm). 
rt  =  -Jx,2  +h 2  ,  where  x,  =  (/  - 1  -  ( Ar  - 1)  /  2)  ■  Ax  and  h  is  the  perpendicular  distance  from  the  center  of  the  linear  fiber. 


Ax  is  the  length  step  of  point  sources  and  N  (odd  integer)  is  the  number  of  equal  spaced  point  sources  used  in  the 
summation  (parenthesis  in  Eq.  (2)).  The  numerical  value  of  the  summation  should  be  independent  of  N  (or  Ax)  if  N  is 
large  enough.  We  found  that  accurate  results  of  the  summation  can  be  obtained  if  Ax  $0. 1  cm  or  A  >  25  for  /  =  2.5  cm. 
In  all  our  calculations  A  =  201  was  used.  The  two  free  parameters  (pa  and  ps’)  are  inherently  separable  in  that  for  a 
CDF  with  a  given  length:  the  magnitude  of  the  fluence  rate  near  the  light  source  ( h  =  0)  is  determined  by  ps’  only  and 
the  slope  of  the  spatial  decay  of  the  light  fluence  rate  is  determined  by  pefronly. 

In  theory,  measurements  of  tj>  at  two  different  distances  r  from  the  point  source  with  known  unit  length  source 
strength  5  and  length  /  are  sufficient  to  determine  both  pa  and  ps\  Measurements  at  multiple  sites  allow  evaluating  the 
variation  of  these  optical  characteristics  within  the  prostate  volume.  Since  Eq.  (1)  is  a  non-linear  equation  of  two  free 
parameters  pa  and  ps’,  we  used  a  differential  evolution  algorithm  developed  by  Storn  et  al  (29).  This  algorithm  is 
simple  and  robust,  which  converges  faster  and  with  more  certainty  than  adaptive  simulated  annealing  as  well  as  the 
annealed  Nelder  &  Mead  approach  (29).  We  modified  the  algorithm  to  require  that  all  parameters  (pa  and  ps’)  are 
positive  (30). 


Table  1  In-vivo  Optical  properties  measured  at  732  nm  in  human  prostate.  The  values  in  the  parentheses  are  the 
standard  deviation  (s.d.)  of  the  mean  values  measured  from  different  locations  of  the  same  prostate  (’x’  in  Fig.  lb). 
No  s.d.  is  listed  if  only  one  data  point  is  available. _ 


Patient 

number 

pa  (cm'1) 

Before  PDT 
Ps’  (cm'1) 

5  (cm) 

pa  (cm'1) 

After  PDT 

Ps’  (cm'1) 

5  (cm) 

2 

0.09 

29.8 

0.34 

0.09 

43.7 

0.29 

3 

0.15 

22.0 

0.31 

0.07 

33.4 

0.37 

4 

0.43  (0.28) 

7.69  (4.76) 

0.41  (0.14) 

0.51 

1.67 

0.63 

5 

0.21 

11.8 

0.37 

0.13 

7.18 

0.60 

6 

0.27  (0.27) 

10.5(11.2) 

0.50  (0.05) 

0.19(0.20) 

18.9(18.4) 

0.45  (0.06) 

7 

— 

— 

— 

0.30  (0.08) 

23.7(13.9) 

0.24  (0.11) 

9 

0.53  (0.36) 

6.61  (4.51) 

0.41  (0.09) 

0.64  (0.25) 

7.00  (5.59) 

0.33  (0.10) 

10 

0.63  (0.32) 

4.62  (2.87) 

0.42(0.10) 

0.19(0.05) 

9.27  (4.47) 

0.54(0.31) 

11 

0.67  (0.17) 

6.39(3.18) 

0.32  (0.10) 

0.83  (0.45) 

5.45  (3.89) 

0.38(0.16) 

12 

0.71  (0.43) 

8.99(6.51) 

0.32  (0.12) 

0.30  (0.06) 

20.2  (4.8) 

0.28  (0.08) 

13 

0.27(0.14) 

18.5  (11.6) 

0.30  (0.07) 

0.26  (0.09) 

17.0  (8.8) 

0.31  (0.07) 

14 

0.72(0.11) 

3.37(1.37) 

0.39  (0.11) 

— 

— 

— 

Optical  properties  measured  in  13  patients  have  been  summarized  in  Table  1  (12).  The  heterogeneity  of  optical 
properties  in  human  prostates  is  somewhat  smaller  than  that  observed  in  canine  prostates  at  732  nm  (7).  Overall  pa 
varied  between  0.07  -  1.62  cm'1  (mean  0.3±0.2  cm'1)  and  ps’  varied  between  1.1  -  44  cm'1  (mean  14±1 1  cm"1).  The 
effective  attenuation  coefficient  peff  varied  between  0.91  —  6.7  cm'1,  corresponding  to  an  optical  penetration  depth  (5  = 
1/peff)  of  0.2  -  1.1  cm.  The  mean  values  of  pcff  and  5  were  2.9±0.8  cm'1  and  0.4±0.1  cm,  respectively.  This  penetration 
depth  is  substantially  larger  than  that  of  0.1  -  0.25  cm  predicted  for  630  nm  (6)  but  is  smaller  than  0.5  -  3  cm  observed 
in  normal  canine  prostate  at  732  nm  (7).  The  most  probably  explanation  is  that  canine  prostate  has  different 
grandular/structure  content  than  that  of  human  prostate.  While  the  mean  reduced  scattering  coefficient  in  canine  was 
3.6  ±  4.8  cm'1  (7),  it  was  15  ±  1 1  cm'1  in  human  at  the  same  wavelength  (732nm).  The  increased  reduced  scattering 
coefficient  resulted  in  increased  effective  attenuation  coefficient,  or  a  reduction  of  optical  penetration  depth,  assuming 
the  absorption  coefficient  remains  the  same. 

2.  Patient  Selection,  Surgical  and  PDT  Procedure 

A  Phase  I  clinical  trial  of  motexafin  lutetium  (MLu)-mediated  PDT  in  patients  with  locally  recurrent  prostate 
carcinoma  was  initiated  at  the  University  of  Pennsylvania.  The  protocol  was  approved  by  the  Institutional  Review 
board  of  the  University  of  Pennsylvania,  the  Clinical  Trials  and  Scientific  Monitoring  Committee  (CTSRMC)  of  the 
University  of  Pennsylvania  Cancer  Center,  and  the  Cancer  Therapy  Evaluation  Program  (CTEP)  of  the  National  Cancer 
Institute.  A  total  of  15  patients  were  treated.  Each  patient  who  signed  the  informed  consent  document  underwent  an 
evaluation,  which  included  an  MRI  of  the  prostate,  bone  scan,  laboratory  studies  including  PSA  (prostatic  specific 
antigen),  and  a  urological  evaluation.  Approximately  two  weeks  prior  to  the  scheduled  treatment  a  transrectal 
ultrasound  (TRUS)  was  performed  for  treatment  planning.  An  urologist  drew  the  target  volume  (the  prostate)  on  each 
slice  of  the  ultrasound  images.  These  images  were  spaced  0.5  cm  apart  and  were  scanned  with  the  same  ultrasound  unit 
used  for  treatment. 

A  built-in  template  with  a  0.5-cm  grid  projected  the  locations  of  possible  light  sources  relative  to  the  prostate.  A 
treatment  plan  was  then  prepared  to  determine  the  location  and  length  of  light  sources.  Cylindrical  diffusing  fibers 
(CDF)  with  active  lengths  1-5  cm  were  used  as  light  sources.  The  sources  were  spaced  one  centimeter  apart  and  the 
light  power  per  unit  length  was  less  than  or  equal  to  150  mW/cm2  for  all  optical  fibers.  The  length  of  the  CDF  at  a 
particular  position  within  the  prostate  was  selected  to  cover  the  full  length  of  the  prostate  (see  Fig.  la).  The  final  plan 
often  required  that  the  prostate  be  divided  into  four  quadrants.  Four  isotropic  detectors  were  used,  each  placed  in  the 
center  of  one  quadrant.  A  fifth  isotropic  detector  was  placed  in  a  urethral  catheter  to  monitor  the  light  fluence  in  the 
urethra.  (Fig.  lb). 

The  patients  were  anesthetized  in  the  operating  room  with  general  anesthesia  to  minimize  patient  movement  during 
the  procedure.  Transrectal  ultrasound  unit  was  used  to  guide  needle  placement  in  the  operating  room.  A  template  was 
attached  to  the  ultrasound  unit  and  was  matched  to  the  same  5-ram  grid  used  for  treatment  planning.  Four  detector 
catheters  (one  for  each  quadrant)  were  inserted  into  the  prostate.  These  detectors  were  kept  in  place  during  the  entire 


procedure  of  PDT  treatment.  Four  additional  pre-planned  treatment  catheters  for  light  sources  were  then  inserted  0.5  or 
0.7  cm  away  from  the  detector  catheters  (Fig.  lb).  These  source  catheters  were  used  for  light  delivery  and  optical 
properties  measurements.  A  15-W  diode  laser,  model  730  (Diomed,  Ltd.,  Cambridge,  United  Kingdom)  was  used  as 
the  732  nm  light  source. 


CDF 


Template 


Prostate 


(A) 


Figure  1.  (a)  Experimental  setup  for 
measuring  the  in-vivo  optical  properties 
of  canine  prostate.  The  prostate  template 
was  drilled  with  a  5-mm  equal  spaced 
grid.  Cylindrical  diffusing  fibers  (CDF) 
were  inserted  into  the  catheters  to 
illuminate  the  entire  prostate  gland.  An 
isotropic  detector  (not  shown)  is  placed 
in  one  of  the  catheters,  which  is  moved 
to  locations  at  a  distance  of  3,  6,  9,  12, 
15  mm  from  the  light  source.  The 
detector  reading  at  each  location  is 
peaked  to  ensure  that  it  is  at  the  middle 
of  the  linear  fiber. 


3.  Search  procedure  with  the  Cimmino  optimization  algorithm 

At  present  a  medical  physicist  chooses  the  number  of  light  sources,  the  particular  template  holes  (or  “slots”)  for 
source  insertion,  the  length  of  the  linear  light  source  to  insert  in  each  slot,  the  position  of  the  source  within  the  slot,  and 
a  single  duration  of  illumination  for  the  entire  set  of  sources.  This  is  a  tedious  and  time-consuming  job  that  requires 
contours  from  ultrasonic  tomographic  images,  visualization  of  3D  volumes  (prostate,  urethra,  and  rectum)  and  their 
intersections  with  linear  light  sources,  estimation  of  the  mean  opacity  of  the  prostate,  and  visualization  of  the  scattered 
light  distribution  within  the  prostate,  for  different  choices  of  source  parameters. 

The  algorithms  discussed  below  try  to  achieve  (1)  a  prescribed  minimum  dose  within  the  prostate,  and  (2)  doses 
not  exceeding  the  maximum  doses  specified  separately  for  the  prostate,  urethra,  rectum,  and  background  tissues.  The 


contours  of  the  prostate,  urethra,  and  rectum  in  each  transverse  slice  (parallel  to  the  template  plane  and  perpendicular  to 
the  linear  light  sources)  are  assumed  available  in  computer  memory.  At  present,  up  to  13  transverse  slices  spaced  5 
mm  apart  are  allowed. 

The  template  currently  being  used  for  source  insertion  is  a  plate  with  a  square  array  of  13  x  13  =  169  holes  (slots 
for  linear  light  sources)  spaced  5  mm  apart  (Fig.  lb).  For  the  particular  patient  data  being  used  as  a  benchmark,  only  51 
slots  are  situated  to  allow  the  light  source  either  to  penetrate  the  prostate  or  approach  within  a  1-mm  margin.  In  present 
clinical  practice,  sources  are  separated  by  10  mm,  that  is,  every  other  slot.  Thus  for  the  benchmark  patient,  12  template 
slots  are  used  for  sources. 

The  light  source  is  a  tube  of  illumination,  1  mm  in  diameter  and  at  most  50  mm  in  length  from  template  base  to 
maximum  penetration.  Creating  the  illumination  within  the  tube  are  “light  seeds”  of  5  mm  length  (thus  5  mm  between 
the  centers  of  adjacent  light  seeds).  The  algorithms  of  this  paper  enforce  two  clinical  requirements:  (1)  the  seeds  within 
a  light  source  are  contiguous  (i.e.,  no  gaps  occur  between  light  seeds),  and  (2)  each  light  source  has  at  least  two  seeds. 
Although  violation  of  these  restrictions  might  yield  mathematically  improved  light  distributions,  the  clinical  use  of 
short  discontinuous  light  sources  requires  greater  precision  and  increased  time  for  in  vivo  placement,  thus  an  increased 
risk  to  the  patient. 

There  are  four  mathematical  problems.  (1)  Given  linear  light  sources  with  specified  slots  and  specified  parameters 
(source  length,  and  retraction  of  the  first  seed  into  the  slot),  find  the  source  strengths  (or  product  of  power  and  durations 
of  illumination  in  units  of  energy)  to  satisfy  the  prescribed  dose  constraints.  (2)  Given  the  source  slots,  find  the  source 
parameters  and  source  strengths.  (3)  Given  only  the  number  of  sources  and  the  allowed  set  of  slots,  find  the  source 
slots,  source  parameters,  and  source  strengths.  (4)  Guarantee  the  convergence  of  each  mathematical  procedure  to  an 
acceptable  solution;  if  no  solution  exists  (that  is,  not  every  constraint  can  be  satisfied),  then  guarantee  convergence  to  a 
least-squares  (best  compromise)  solution.  These  four  problems  must  be  solved  for  two  cases:  individual  source 
strengths  (sources  may  have  different  strengths),  and  uniform  source  strength  (all  sources  have  the  same  strength).  The 
case  of  uniform  source  strength  with  all  source  slots  and  parameters  specified  (problem  1)  is  the  present  practice  in  the 
clinic  and  is  our  baseline  for  improvement. 

Our  search  algorithm  systematically  changes  the  position  and  retraction  of  each  light  source  within  the  prostate, 
calculates  the  best  source  strengths  with  the  Cimmino  algorithm  (after  each  change),  and  then  checks  whether  the 
discrepancy  from  the  dose  prescription  has  decreased.  The  discrepancy  function  we  use  is  the  sum  of  the  absolute 
values  of  the  overdose  or  underdose  at  each  sample  point  with  respect  to  the  prescribed  constraint.  An  importance 
weight  is  pre-chosen  for  each  named  volume  (prostate,  urethra,  rectum,  and  background)  and  normalized  so  that  the 
sum  of  the  importance  weights  is  unity.  The  importance  weight  of  a  sample  point  is  the  importance  weight  of  its 
volume  divided  by  the  number  of  sample  points  in  that  volume.  Thus  the  discrepancy  at  each  sample  point  is  given  a 
normalized  importance  weight.  If  there  is  no  solution  with  zero  discrepancy,  the  normalized  weighting  will  put  most  of 
the  error  where  it  is  least  important. 

A  feasibility  algorithm  rigorously  solves  the  first  problem.  The  second  and  third  problems  require  search 
procedures  over  allowed  slots,  source  lengths,  and  retractions  (as  appropriate).  For  each  choice  of  source  slots  and 
source  parameters  that  is  assayed  in  a  search,  the  discrepancy  between  the  solution  of  the  feasibility  algorithm  and  the 
prescribed  dose  at  a  point,  summed  over  the  constraint  points,  is  calculated.  If  this  total  discrepancy  is  less  than  any 
previously  discovered,  then  that  parameter  set  becomes  the  current  best  solution.  The  second  and  third  problems  are 
“combinatoric”,  where  exhaustive  searches  are  not  possible  in  reasonable  time.  Non-exhaustive  searches  risk 
encountering  local  minima,  so  that  finding  the  solution  for  the  “absolute-minimum”  discrepancy  cannot  be  guaranteed 
or  even  recognized.  Thus  the  key  to  the  second  and  third  problems  is  a  good  search  strategy. 

Solution  of  the  first  problem  allows  comparison  of  individual  strengths  vs.  a  uniform  strength  for  user-chosen 
sources.  Solution  of  the  second  problem  allows  a  similar  comparison  but  with  source  lengths,  retractions,  and  strengths 
chosen  by  prostate  geometry  and  minimum  dose  discrepancy.  A  solution  of  the  third  problem  allows  the  treatment 
planner  to  find  automatically  the  best  (or  almost-best)  light  sources  and  source  strengths  from  just  the  specified  dose 
constraints  to  the  prostate  and  organs,  and  the  choice  for  the  number  of  sources.  Such  a  solution  would  be  virtually 
impossible  to  find  by  human  visualization  and  estimation.  Automatic  source  selection  and  weighting  becomes  even 
more  important  if  it  is  necessary  to  limit  dose  to  the  urethra  and  rectum  without  compromising  treatment  effectiveness 
to  the  prostate. 

The  selection  of  the  number  of  light  sources  is  as  much  a  clinical  as  a  mathematical  problem.  To  find 
mathematically  the  best  number  of  linear  light  sources  to  insert,  one  can  rerun  repeatedly  any  procedure  used  to  solve 
problem  3  with  different  numbers  of  sources.  However  changing  the  number  of  sources  differs  from  rearranging 
source  positions  and  solving  for  source  strengths;  it  involves  a  tradeoff  between  (a)  fewer  sources  —  less  homogeneous 


dose  coverage,  higher  source  strengths,  but  fewer  surgical  complications,  and  (b)  more  sources —  better  dose  coverage, 
lower  source  strengths,  but  more  surgical  complications).  Thus  choosing  the  number  of  sources  requires  clinician  input 
based  on  medical  experience  and  judgment. 

To  focus  this  paper  on  algorithmic  procedures,  we  omit  any  discussion  on  choosing  the  number  of  light  sources. 
The  number  of  sources  is  always  assumed  to  be  given.  Several  elements  are  common  to  the  last  three  algorithms. 
These  are  (1)  a  kernel  (Eq.  2)  to  calculate  dose;  (2)  a  set  of  constraint  points  in  the  volume  of  interest  to  guide  the 
optimization  procedure;  these  points  derive  from  both  a  sampling  grid  and  extra  points  chosen  on  the  periphery  of  the 
prostate  (target);  (3)  the  Cimmino  feasibility  algorithm. 

Each  search  procedure  allows  the  number  of  sources  to  be  decreased  by  one.  A  new  optimization  calculation  is 
then  performed  and  the  resulting  dose  distribution  compared  with  the  previous.  The  decrease  in  the  number  of  sources 
can  be  continued  iteratively.  If  the  number  of  sources  can  be  decreased  without  increasing  discrepancy  between  the 
prescribed  and  optimized  dose  distribution,  fewer  sources  need  be  used,  thereby  reducing  the  complication  of  the 
procedure  and  discomfort  of  the  patient. 

Each  algorithm  also  allows  the  template  to  be  shifted  horizontally  or  vertically  (anterior  to  supine  patient).  Up  to 
seven  shifts  are  allowed  in  each  direction  (e.g.  from  -3  mm  to  +3  mm  in  steps  of  1  mm  for  the  5-mm  grid),  thus  49 
separate  optimizations.  These  can  be  used  to  check  the  sensitivity  of  the  dose  distribution  to  small  shifts  in  template 
position. 

To  check  the  effect  of  Cimmino  optimizations,  two  different  optical  properties  were  chosen:  (a)  the  average  optical 
properties  of  all  prostate  patients,  pa  =  0.3  cm'1  and  ps’  =  14  cm'1;  and  (b)  an  extreme  case  when  the  optical  penetration 
depth  is  the  longest  observed  in  human  prostate  (8  =  1.3  cm),  pa  =  0. 1  cm'1  and  ps’  =  2  cm'1 . 

Table  2.  Results  of  source  strengths  (J)  obtained  using  various  Cimmino  search  algorithms  for  1 2  sources.  Standard 
plan:  user  selected  source  parameters  and  uniform  strengths;  Cimmino  1:  optimize  source  strengths  only; 
Cimmino  2:  optimize  source  lengths,  strengths,  and  template  slots;  and  Cimmino  3:  the  same  as  Cimmino  2  with 
constraint  for  rectum.  Opt  1:  pa  =  0.3  cm'1,  ps’  =  14  cm'1;  Opt  2:  pa  =  0.1  cm'1,  ps’  =  2  cm'1.  The  Calculation  time  is 
obtained  from  a  Dell  PC  with  a  2.8  GHz  Pentinum  IV  processor. 


Source  No. 

Standard  plan 
Opt  1  Opt  2 

Cimmino  1 

Opt  1  Opt  2 

Cimmino  2 

Opt  1  Opt  2 

Cimmino  3 

Opt  1  Opt  2 

1 

410.83 

69.56 

286.22 

43.36 

267.66 

57.48 

296.82 

57.48 

2 

410.83 

69.56 

232.24 

63.72 

271.56 

51.48 

292.12 

51.48 

3 

410.83 

69.56 

264.92 

67.04 

222.82 

49.06 

312.78 

49.06 

4 

410.83 

69.56 

242.40 

40.74 

213.14 

78.18 

240.62 

78.18 

5 

410.83 

69.56 

350.54 

64.38 

349.38 

50.20 

344.94 

50.20 

6 

410.83 

69.56 

293.08 

83.90 

358.48 

55.54 

244.36 

55.54 

7 

410.83 

69.56 

595.94 

75.60 

219.50 

56.84 

239.14 

56.84 

8 

410.83 

69.56 

429.26 

80.60 

293.76 

81.30 

159.44 

81.30 

9 

410.83 

69.56 

306.20 

71.52 

219.16 

44.42 

171.24 

44.42 

10 

410.83 

69.56 

289.82 

47.46 

351.70 

56.56 

229.26 

56.56 

11 

410.83 

69.56 

290.28 

51.14 

258.80 

56.44 

241.88 

56.44 

12 

410.83 

69.56 

391.96 

77.00 

209.24 

69.92 

183.74 

69.92 

Calculation  time  (s) 

1 

1 

9 

9 

260 

173 

309 

173 

III.  RESULTS  AND  DISCUSSIONS 

We  used  our  standard  plan,  based  on  geometrical  coverage,  1-cm  spaced  loading,  and  uniform  source  strength  as 
our  default  plan  to  judge  the  improvement  made  by  the  Cimmino  algorithms. 

Figure  2  compares  computer  runs  of  optimized  100%  isodose  distributions  to  a  single  3D  prostate  volume  using 
the  Cimmino  algorithm  (solid  lines)  compared  with  the  standard  plan  (dashed  line).  The  source  locations  for  the 
Cimmino  run  were  chosen  manually  and  kept  the  same  as  the  uniform  loading  plan  (dashed  line).  There  is  no 
substantial  difference  between  the  equal-source-strength  plan  and  the  Cimmino  individual-source-strength  plan  for 
either  the  mean  or  most  penetrating  optical  properties  found  in  the  human  prostate.  This  is  not  surprising  since  when 


the  light  penetration  depth  is  short  enough,  the  light  coverage  is  determined  more  by  the  light  source  location  than  the 
source  strength,  provided  we  have  renormalized  the  plan  to  guarantee  that  the  entire  prostate  is  covered. 


z  =  0  cm  z  =  0.5  cm  z  =  1  cm  z  =  1 .5  cm 


Figure  2.  Comparison  of  100%  isodose  lines  for  uniform-strength  (solid  line)  vs.  Cimmino  optimized  loading  (dashed 
line)  —  for  fixed  template  and  CDF  source  lengths —  for  two  optical  properties:  left  (a)  pa  =  0.3  cm4,  ps’  =  14  cm'1  and 
right  (b)  pa  =  0.1,  ps’  =  2  cm'1.  Contours  for  prostate,  urethra,  and  rectum  are  also  shown  in  the  figure. 

To  further  explore  the  importance  of  the  source  location,  we  let  the  Cimmino  algorithm  determine  the  source 
locations,  lengths,  as  well  as  the  source  strengths.  We  kept  the  total  number  of  sources  to  be  the  same  as  the  standard 
plan,  12.  The  results  are  shown  in  Figure  3.  The  source  positions  (o)  are  chosen  by  the  computer  with  the 
corresponding  isodose  line  plotted  as  solid  line.  The  source  positions  for  the  standard  plan  are  the  same  as  those  shown 
in  Fig.  2.  For  both  optical  properties,  the  Cimmino  algorithm  provides  a  better  plan,  sparing  the  urethra,  rectum,  and 
background  while  covering  the  entire  prostate.  The  light  source  positions  are  now  chosen  more  peripherally  compared 
to  the  equally  spaced  coverage,  i.e.,  fewer  sources  are  placed  near  the  urethra  compared  with  the  uniform  loading.  The 


improvement  is  more  significant  for  shorter  light  penetration  depth  (Fig.  3a)  than  that  for  longer  light  penetration  depth 
(Fig.  3b). 


z  =  0  cm  z  =  0.5  cm  z=  1  cm  z  =  1 .5  cm 


Figure  3  Comparison  of  100%  isodose  lines  of  Cimmino  optimized  results  (optimized  source  lengths,  loading,  and 
template)  for  12  linear  sources  (solid  lines)  and  Cimmino  optimized  results  (optimized  loading  only,  user  selected 
source  lengths  and  template)  (dashed  lines)  for  two  optical  properties:  (a)  pa  =  0.3  cm'1,  ps’  =  14  cm"1  and  (b)  pa  = 
0.1,  ps’  =  2  cm'1.  The  source  location  for  the  dashed  lines  is  the  same  as  Fig.  2. 

We  then  examined  whether  maximum  dose  constraint  for  the  critical  structure  makes  any  difference.  We 
compared  two  plans  made  with  the  Cimmino  algorithm  (computer  choosing  positions,  lengths,  and  source  strengths) 
but  one  had  the  maximum  dose  to  the  rectum  reduced  from  300%  (dashed  line)  to  200%  (solid  line).  Figure  4a  shows 
that  the  rectum  sparing  is  improved  for  the  average  optical  properties  in  prostate,  which  is  less  than  5  mm.  No 
significant  improvement  is  observed  if  the  optical  penetration  depth  is  more  than  1  cm  (Fig.  4b). 

Finally  we  examined  the  best  plan  we  have  obtained  so  far  using  12  CDF  sources  (Fig.  4,  solid  line)  and  compared 
it  with  the  Cimmino  algorithm  using  51  CDF  sources.  The  latter  represents  the  best  mathematical  plan  for  linear 


sources  and  the  given  template,  given  no  discontinuity  along  any  catheter  for  active  light  sources.  The  result  shows  that 
the  rectal  sparing  of  the  12-source  plan  is  not  much  worse  than  that  of  the  plan  using  51  sources  for  low  opacities  (Fig. 
5b).  For  the  average  optical  parameters  (Fig.  5a),  some  improvement  of  rectum  sparing  can  be  achieved  by  using  more 
CDF  sources. 


z  =  0 cm  z=  0.5  cm  2=1  cm  z  =  1 .5 cm 


Figure  4  Comparison  of  100%  isodose  lines  of  Cimmino  optimized  results  for  12  linear  sources  with  0.5  cm  step 
size  with  the  maximum  dose  constraint  of  200%  (solid  lines)  and  300%  (dashed  lines)  for  rectum  for  two  optical 
properties:  (a)  pa  =  0.3  cm'1,  ps’  =  14  cm'1  and  (b)  pa  =  0.1,  ps’  =  2  cm"1.  In  this  example,  the  source  locations, 
lengths,  retractions,  and  strengths  are  chosen  by  the  Cimmino  algorithm  (see  Table  2). 

Figure  6  shows  the  effect  of  optical  properties  on  the  prostate  light  dose  coverage  when  incorrect  optical  properties 
are  used.  In  this  figure,  same  source  strengths  and  source  parameters  are  used  for  two  different  optical  properties.  The 
100%  isodose  line  provides  adequate  prostate  coverage  for  one  optical  properties  (solid  line)  but  is  very  bad  if  the 
optical  properties  is  in  fact  shorter  than  what  is  assumed  (dashed  line).  This  indicates  that  the  values  of  the  optical 
parameters  are  very  important  in  determining  the  actual  light  fluence  distribution. 


The  present  paper  assumed  uniform  (homogeneous)  optical  properties.  An  open  question  is  whether  optimized 
solutions  for  inhomogeneous  media  will  further  improve  over  the  present  uniform-medium  calculations.  It  is  also  an 
open  question  on  the  number  of  control/constraint/sample  points  required  to  guarantee  the  optimization  outcome. 
Further  studies  are  needed  to  determine  the  minimum  number  of  CDF  needed  to  achieve  complete  coverage. 


z  =  0  cm  z=0.5cm  z  =  1  cm  z=1.5cm 


Figure  5  Comparison  of  100%  isodose  lines  of  Cimmino  optimized  results  (solid  lines)  for  51  linear  sources  with 
0.5  cm  spacing  and  the  best  Cimmino  plan  using  12  sources  (dashed  lines,  loading  pattern  in  Fig.  4)  for  two  optical 
properties:  (a)  pa  =  0.3  cm'1,  ps’  =  14  cm'1  and  (b)  pa  =  0.1,  ps’  =  2  cm'1. 

IV.  CONCLUSION 

The  question  addressed  is  whether  any  significant  advantage  may  derive  from  methods  that  weight  each  CDF 
source  separately  and/or  choose  the  geometry  of  the  light  sources  as  well.  In  summary,  our  comparison  shows  that:  (1) 
It  is  important  to  measure  the  optical  properties  of  a  patient  because  it  determines  the  light  fluence  distribution.  This 


4 


effect  is  more  predominant  than  optimizing  the  source  positions,  lengths,  and  strengths.  (2)  For  the  range  of  the  optical 
properties  in  the  human  prostate,  the  individual-source  strength  optimization  does  not  significantly  improve  over  the 
equal-source-strength  optimization  when  given  the  positions  of  the  light  sources.  (3)  Improved  results  with  Cimmino 
optimization  compared  to  user-specified  optimization  indicate  that  computer  optimization  saves  the  user  time  and  setup 
and  reduces  human  stress.  The  Cimmino  optimization  is  fast  enough  for  this  problem  to  obtain  clinical  real-time 
optimization  (less  than  300  s,  Table  2). 


Figure  6.  A  comparison  of  prostate  coverage  for  pa  =  0.1,  ps’  =  2  cm'1  (solid  line)  and  pa  =  0.3  cm'1,  ps’  =  14  cm'1 
(dashed  line)  using  the  source  strengths  listed  in  Table  2  for  Cimmino  3,  Opt  2. 
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